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Motion of Individual Domain Walls in a 
Nickel-Iron Ferrite 


By J. K. GALT 
(Manuscript received May 11, 1954) 


Samples have been cut from single crystals of the nickel-iron ferrite 
CNiO)o.75 (FeO)o.2F e203 in such a way that they contain one and only one 
movable ferromagnetic domain wall. The viscous damping coefficient for 
this wall, which is a measure of the losses associated with domain wall 
motion in this material, has been measured as a function of temperature. 
This damping shows a very large increase as the temperature goes down to 
the region of 77°K. The value of this damping is correlated with the Landau- 
Inifshitz equation for the rotational motion of magnetization by means of 
previously available theoretical analysis. In addition, it 1s suggested that 
the sharp increase in damping at low temperatures is due to a relaxation 
associated with a rearrangement of the valence electrons on the divalent and 
trivalent tron tons in the ferrite. A tentative phenomenological theory of the 
losses based on this mechanism ts presented. 


INTRODUCTION 


The mechanism which contributes most to the permeability of high- 
permeability magnetic materials is the motion of ferromagnetic domain 
walls. These walls are thin lamellae in which the direction of the spon- 
taneous magnetization of the material changes from one domain to 
another. As a result, this mechanism contributes a major part of the 
energy losses which accompany rapid changes in the direction of the 
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magnetization in such materials. In the ferromagnetic metals, it is well 
known that these losses ordinarily arise largely from the eddy currents 
which are induced by the motion of the domain walls. In the ferrites, 
however, the conductivity is so low that the contribution of eddy cur- 
rents to the losses is never overwhelming and is often negligible; the 
losses must therefore in large part arise from other sources not yet under- 
stood. It is the purpose of this paper to present some recent studies of 
these losses and to discuss their relevance to the losses in ferrites gen- 
erally. 

In any ordinary sample of a ferromagnetic material, a study of domain 
wall motion and the associated energy losses is complicated by the fact 
that the domain pattern is very complex. Any attempt to provide a 
theoretical explanation of data taken on such samples must involve an 
averaging process over many domain walls of varying area, crystal 
orientation, etc. This makes it extremely difficult to describe the be- 
havior of such patterns uniquely and quantitatively, although some 
progress has been made.’” A method of avoiding this difficulty has been 
developed by Williams, Bozorth, Shockley, Kittel? and Stewart** in 
working on silicon iron. This method consists in cutting a polygonal 
ring from a single crystal in such a way that each leg of the ring hes 
along one of the easy directions of magnetization in the crystal. In 
silicon iron this leads to a rectangular ring with each leg along a [100] 
crystal direction. In the ferrite which we use this technique to study 
here, the easy directions are [111] directions, and we use a diamond 
shaped sample as shown by the solid lines in Fig. 1. Each leg is along a 
[111] direction, and the major face is a (110) plane. If the sample is good 
enough, the domain pattern is that indicated by the dotted lines in Fig. 
1. This pattern consists of four stationary walls, one at each corner, and 
one movable wall which goes all the way around the sample. The mag- 
netization thus travels around the sample in two paths, one clockwise 
and the other counter-clockwise, and the position of the movable wall 
therefore determines the net circumferential magnetization. In such 
samples we study quantitatively and in some detail the motion of an 
individual movable wall. 

Williams, Shockley and Kittel’ studied the motion of the movable 
wall on one of the rectangles cut from a single crystal of silicon iron. 
They found the motion to be viscously damped, as Sixtus and Tonks* 
had in earlier experiments with more complicated domain walls. Be- 
cause of the simplicity of their domain pattern, Williams, Shockley and 
Kittel were able to calculate the eddy current losses in their experi- 
ments, and to show that they accounted for most of the observed damp- 
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ing as expected. There was an additional contribution, however, which 
Kittel suggested was due to mechanisms of the sort which give rise to 
the width of ferromagnetic resonance lines.’ These mechanisms of course 
are the controlling ones in the ferrites, where eddy current losses are 
small. The motion of a domain wall damped by such effects and un- 
affected by eddy currents was first discussed in a classic paper by Landau 
and Lifshitz.” 

The experiments reported in the present paper consist of measure- 
ments of the velocity of a movable wall as a function of applied magnetic 
field in a sample like that shown in Fig. 1. The measurements are made 
by observing the voltage induced in a secondary winding on such a 
sample when a known field is applied by means of a pulse of current 
in a primary winding. The composition of the ferrites used in these 
studies is given by the approximate chemical formula (NiO)9.75(FeO)o.25- 
Fe,0;. Data have been taken as a function of temperature on several 
samples. The large, perfect crystals of the ferrites which are essential 
to the success of these experiments have been obtained through Dr. 
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Fig. 1 — Sample of ferrite used in this study. 
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G. W. Clark from the Linde Air Products Company. Similar studies 
have been performed previously at room temperature on FeO,’ by 
the author, and in a preliminary way on (NiO)o.1 (FeO)o.09 Fe2O3 ° by 
the author in collaboration with J. A. Andrus and H. G. Hopper. 


THE EXPERIMENTS 
Preparation of Samples 


The key to the success of experiments of this sort is of course obtain- 
ing the domain pattern shown in Fig. 1, so that we have only one movable 
wall in the sample. The achievement of this pattern, and the observa- 
tion of it when achieved, depend in turn on success in producing a 
perfect or almost perfect sample. It therefore seems worth while to 
describe the process which has finally emerged as a satisfactory way of 
producing these samples. 

The rough crystal is first oriented by means of X-rays (Laue and X-ray 
goniometer techniques) to an accuracy of a few minutes while it is 
mounted in such a position that afterwards we can grind a flat on it 
which is coincident with the (110) plane. This flat is ground simply with 
a belt grinder. A cut is then made with a diamond saw parallel to this 
flat, so that we have a disc whose faces are (110) crystal planes. This 
disc is usually made about 24 to 3 mm thick. The second face is ground 
accurately parallel to the first in a paralleling block. A flat coincident 
with the (100) plane is ground on the edge of this disc for later use in 
orienting the sample in the plane of the disc. This flat is also ground with 
a belt grinder after orienting the disc with Laue and X-ray goniometer 
techniques. 

The disc is ground as smooth as possible with 30314 emery, and then 
very carefully polished. The polishing is done first on a lap surfaced 
with No. 0000 french emery paper, with Linde A abrasive loose on top 
of the emery paper. After this the lap is surfaced with a sheet of very 
smooth paper and then Linde B abrasive is used. The polishing process 
takes four to eight hours per disc, and removes all pits visible under 
50X magnification, except those inherent in the crystal. Only on such a 
smooth surface can the very small holes which sometimes occur in these 
crystals be seen. Sometimes, however, the polishing process conceals 
fine cracks. It also cannot reveal variations in chemical composition 
which sometimes occur from point to point in the crystals. Such com- 
position variations are presumably variations in the concentration of 
divalent iron from point to point in the crystal. In order to reveal these 
latter imperfections, the disc is etched for three hours by boiling it in 
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50 per cent H.SO, under a reflux condenser. An asbestos pad is placed 
between the flame and the bottom of the flask containing the H.SO, 
in order to prevent sharp temperature fluctuations in the bath. The rate 
of the etching attack, and the quality of the surface it leaves are sharply 
dependent upon temperature. It should be mentioned that if this etch 
is used on the discs before polishing, the surface remains rough or may 
even be made rougher, so that it is impossible to detect the imperfec- 
tions in the disc. The etch must start on a smooth surface. 

Once the disc is cut, polished, and etched, if it is found to be sufficiently 
free of imperfections, a sample is cut from it in such a way as to avoid 
those which there are, as they are revealed by the polishing and etching 
processes. First the diamond-shaped hole is cut by means of a jig whose 
rotational position with respect to the disc is determined from the (100) 
flat on the edge of the disc. This jig is a piece of steel which is driven in 
vibration vertically with a magnetostrictive drive.’ The surface of the 
disc is covered with a slurry of carborundum or diamond dust, and this 
abrasive is made to cut a hole in the disc as the vibrating jig is slowly 
lowered. With the hole cut in the proper orientation, the outer parts of 
the disc are ground down to form the legs of the sample. Another jig of 
the proper shape 1s used to hold the sample in position during this process. 

A hysteresis loop is taken as soon as the sample is cut. A relatively . 
good loop taken on our best sample is shown in Fig. 2. All such loops 
on these samples are taken on the Cioffi recording fluxmeter.” This 
loop is obviously not yet in the form which we finally need. In order to 
square the hysteresis loop, we anneal the sample for approximately an 
hour at 600°C in a magnetic field of 10 to 20 oersteds. The field is pro- 
duced by running a current through a few turns of glass insulated wire 
wound on the sample. After such a heat treatment the hysteresis loop 
of this sample assumed the form shown in Fig. 3. 

Once the sample is prepared, the next problem is to observe the domain 
pattern and find if any important deviations from the pattern shown 
in Fig. 1 occur. The heat-treatment we give them corrodes the polished 
surfaces of the sample, and of course the faces exposed when the sample 
is cut from the disc have not yet been polished. Consequently both the 
major (110) faces of the sample and the outer faces of the legs are pol- 
ished, and the sample is then etched in the same way as before. Usually 
a hysteresis loop is again taken at this point as a check. If the sample 
is good, it is not significantly different from the loop taken immediately 
after heat-treatment. The sample is brought to a demagnetized con- 
dition at this point so that the movable wall will be near the center of 
the sample where it can be observed. This completes the process of 
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Tig. 2 — Hysteresis loop taken on sample before heat treatment in a magnetic 
field. 


preparing the samples for observing their domain pattern and for per- 
forming our experiments on them. 

Four samples have been prepared in this way for our studies on 
(NiO)o.75 (FeO)o.25 FeO; - 


Domain Pattern Observations 


The method used to observe the domain walls on these surfaces is 
the same as that used by Williams and his collaborators.’ We will there- 
fore not describe it in detail. It consists essentially of observing through 
a microscope the pattern formed by a magnetic colloid on the surface. 

The observation of domain patterns, even on these carefully pre- 
pared samples, is difficult. There is still some pitting on the surfaces. 
Also, many of the surfaces become rounded in the process of polishing. 
This produces surface spikes of the sort discussed by Williams, Bozorth 
and Shockley.” The result is that on many surfaces the domain pattern 
of the sample as a whole has to be discerned in a substantial amount 
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of extraneous structure such as surface spikes and pits. It is therefore 
impossible to show in one picture the whole pattern as diagrammed 
in Fig. 1. However, more detailed pictures of parts of the pattern do 
show that it is there. The essential features of the pattern on our best 
sample are shown in Figs. 4 and 5. Fig. 4 shows the stationary wall at 
one corner, and Fig. 5 shows a section of the movable wall on one leg. 
Differentiation of the domain walls in Figs. 4 and 5 from the many 
scratches is not very difficult after one has some experience in such 
observations. 

The variation in wall position along the leg shown in Fig. 5 is due to 
the effects of strains and other imperfections, present even in this care- 
fully prepared sample, in determining the position of the wall at rest. 
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Fig. 8 — Hysteresis loop taken on same sample as loop in Fig. 2 after annealing 
for one hour at 580°C in a magnetic field of approximately 20 oersteds. 
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Fig. 4 — Picture of domain pattern at one corner of sample. The stationary 
wall is indicated by arrows. 


MOTION OF INDIVIDUAL DOMAIN WALLS 1031 





Fig. 5 — Picture of a section of the movable wall along one leg of sample. The 
line has been slightly emphasized by retouching after the picture was taken. The 
edge of the leg can be seen at the bottom of the figure. 


It is unlikely that the wall is so distorted from a plane when in rapid 
motion, however, since then the driving force and the viscous damping 
resistance to motion are both much larger than the effects of these 
imperfections. The imperfections, of course, are primarily effective in 
determining the coercive force, as read from the hysteresis loop. 

The domain pattern as traced out on each of four samples is discussed 
below: 

Sample 1. Although they had spikes associated with them, the sta- 
tionary walls expected at the corners could be seen, at least in part. In 
addition, at one of the acute angle corners there was some rather ex- 
tensive domain wall structure. This structure had one form when the 
sample was magnetized in one direction, another when it was magnetized 
in the other. It was due to the presence of a small void at this corner 
whose magnetic energy was reduced by having domains of reversed 
magnetization around it. The existence of this void was established from 
structure which appeared in the spots on an X-ray Laue photograph 
taken at this point. The process of magnetization in this sample con- 
sisted in (a) the growth of the wall from the nucleus around this void 
until it existed all around the ring, and (b) the motion of this wall to 
the other side of the sample, where it again shrank to a configuration 
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which minimized the magnetic energy of the void. As a result of this 
state of affairs, the wall was not in equilibrium in the middle of the 
sample, but always shrank around the void so as to magnetize the sample 
in one direction or the other. It was impossible therefore to demagnetize 
the sample so that the wall was at the center of the legs and then have 
the wall stay there to be observed. The wall could be brought to the 
center, however, and held there using the technique mentioned by 
Williams and Shockley® (see Fig. 8 in their paper). The legs of the sample 
were so small that it was difficult to do this, but the wall was found on 
three legs of the sample at different times in this way. The wall curved 
a good deal in traveling along the legs. This sample was etched repeatedly 
so that data could be obtained as a fuction of sample dimensions. The 
variations observed led to a viscous domain wall damping independent 
of dimensions if it was assumed that the domain wall was perpendicular 
to the major (110) face. Therefore it was to this (112) plane that the 
wall was brought for observation. 

Sample 2. Stationary walls at the corners were rather patchy but 
visible. The movable wall was traced along the outside faces of two legs, 
which indicates that it lay in the (110) plane. The wall curved a good 
deal. There were also other walls which enclosed patches of surface. It 
is suspected that these patches were the bases of spike domains extending 
into the sample from strain patterns on the surface. The sample was 
therefore etched again. Unfortunately, the bath apparently became 
locally overheated, and this etch took off rather more material than 
expected. It also left a matte surface on which domain walls could not 
be observed. The data taken on this sample, however, check those on 
other samples if we assume a pattern in which the movable wall is in 
the (110) plane, as our observations lead us to suspect. 

Sample 8. The stationary walls at the corners were seen, but only 
with difficulty. They were patchy. There was a good deal of structure 
all along the legs on the major (110) face of this sample, but no wall 
which ran around the sample could be seen on this face. 

On the outside of two of the legs, which are (112) faces, pitting and 
extraneous walls were so bad that the main wall could not be discerned. 
On a third, the wall could be traced most of the way. On the fourth, 
however, there were two walls, one of which could be traced along the 
whole leg, the other of which went only three-fourths of the way along 
the leg. Both walls on this leg showed a good deal of curvature. It there- 
fore appears that the movable wall lies in the (110) plane, but that 
there is another wall big enough so that it may move and affect our data. 
This picture of a domain pattern with two movable walls was confirmed 
by checking the data obtained on this sample with those from others. 
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Sample 4. The stationary walls at the corners, although patchy, were 
seen. There was some extraneous domain wall structure on the major 
(110) faces, but nothing which looked at all like the main wall. 

On the outside (112) faces of the legs, however, the wall was traced 
almost all the way around. Only short sections were impossible to trace. 
The wall curved as usual, and there was some extraneous domain wall 
structure on these faces, but substantially the whole of the ideal pattern 
shown in Fig. 1 was seen on this sample. The hysteresis loops shown in 
Figs. 2 and 3 and the domain pattern pictures shown in Figs. 4 and 5 
were taken on this sample. 

Not only was the domain pattern on sample 4 the best and most com- 
plete, but the data taken on this sample was much the most reproducible. 
We shall therefore report the data taken on this sample in detail, and 
simply refer to the results on other samples as a check and to indicate 
the sort of variations which occurred from sample to sample. 


Afeasurements 


Our procedure in making the measurements is as follows. The sample 
is wound with a primary and a secondary winding. A square pulse of 
positive voltage is applied to the primary winding in series with a re- 
sistor which is large enough to keep the pulse rise time short. The rise 
time must be short compared to the time required for the field produced 
by the pulse to reverse the magnetization of the sample. On the other 
hand, since the pulse is applied for the purpose of reversing the magneti- 
zation of the sample, the length of the pulse must be at least comparable 
with the time required for the reversal to occur; if possible, it should be 
longer than this. The reversal time, of course, is the time required for 
the mobile domain wall to move from one side of the sample to the other 
under the field produced by the applied pulse. A second pulse, of nega- 
tive voltage, is applied to the primary during each cycle of the pulser 
in order to bring the wall back to its original position so that the phe- 
nomenon may be observed repetitively. 

By synchronizing an oscilloscope sweep with the pulser, the signal 
induced in the secondary winding is observed while the applied pulse 
is on the primary. Since this signal is proportional to the velocity of the 
wall, it is constant to a first approximation during the application of a 
constant field. Irregularities in the crystal may cause the velocity of the 
wall to vary somewhat as it moves across the sample, however, and this 
will cause the signal to vary too. In this case, the observer reads the 
average value. Fig. 6 shows an example of the signal induced in the 
secondary winding as seen on an oscilloscope. Sample 4 was used to 
obtain this picture. 
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WALL VELOCITY IN CM PER SEC 





Fig. 6 — Oscilloscope trace of the induced voltage in the secondary winding 
while a square current pulse 8u sec long was applied to the primary. The voltage 
spike at the beginning is not completely understood, but is thought to be con- 
nected with the fact that the domain wall spikes associated with imperfections 
do not pull back on the domain wall until it has moved a little distance. In this 
sample, once this initial peak is over, the wall velocity comes to its steady state 
value and is not much disturbed by imperfections, so that the observer need do 
no averaging. In this case, the wall was moving with a velocity of 3500 cm/sec. 
At this velocity, the wall was unable to reverse the magnetization in &u sec, so the 
signal ends as the magnetic field goes to zero at the end of the pulse. After the 
pressure on the wall due to the magnetic field stops, the domain wall spikes asso- 
ciated with imperfections pull the wall back slightly, giving rise to the voltage 
spike in the opposite direction. 


The applied field due to the primary pulse is deduced from the cur- 
rent in the primary winding (measured by observing the voltage across 
the series resistor) using the solenoid formula H = 41NI. To obtain the 
relation between wall velocity and induced voltage per secondary turn 
we have: 


Volts/turn = (d®/dt) X 10°° = 8rM,(Az/At)wvon X 10°, (1) 


where (Az/At) is equal to the domain wall velocity v, and wwan is the 
width of the wall between the boundaries of the sample in the direction 
perpendicular to the direction of magnetization. It is in deriving (1), 
of course, that we use our detailed knowledge of the domain pattern in 
the sample. 

We are thus able to obtain a value of the domain wall velocity v for 
each value of the applied field H. These data are the results of the 
experiment. 

The value of M, used in (1) is the measured value (822.5 cgs sania) cc) 
at room temperature. The values used at other temperatures have been 
deduced on the assumption that M, varies the same way with tempera- 
ture in this material as it does in magnetite as measured by Weiss and 
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Forrer.” We have extrapolated their data to get the variation up to our 
highest temperatures. 

In general, a plot of the data turns out to have the form shown in Fig. 
7. This is the data taken on Sample 4 at 201°K. The wall does not move 
until the field exceeds the coercive force required to get it past various 
imperfections in the crystal. Its motion in fields higher than this is 
viscously damped. The wall velocity, v, therefore follows the relation: 


v= GH — Z.,), | (2) 


where H, is the coercive field and G is the slope of the line drawn through 
the data. The value of G is high if the losses are low, and vice versa, of 
course. 


RESULTS 


Data on Sample 4 of the sort shown in Fig. 7 have been taken at 
various temperatures. We show in Fig. 8 a plot of v/(H — H,) as a 
function of temperature for this sample. Clearly, the outstanding fea- 
ture of the data is a tremendous increase in the viscous damping of the 
domain wall at low temperatures. . 

Since the other samples were not as satisfactory, for reasons given 
above, we do not reproduce the data on them explicitly. Similar data 
have been taken on Samples 1 and 2, however, and they show the same 
behavior within their accuracy except that the very sharp decrease in 
v/(H — H,) seemed to occur at a somewhat higher temperature. This 
difference may be due to slight variations in composition among the 
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Fig. 7 — Typical plot of actual data for domain wall velocity as a function of 
applied field. 
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samples (Chemical analysis indicated that Sample 1 was 
| (NiO)o.m(FeO)o.2sFe20s , 


as distinct from 
(NiO)o.75(FeO)o.25!'e203 


for Sample 4) but it is possible that other more subtle differences such as 
the arrangement of the divalent nickel and iron ions are involved. Volt- 
ages induced in the secondary winding on Sample 3 for various applied 
fields were much higher than for the other samples at room temperature. 
This confirms the presence of more than one wall as indicated by the 
domain pattern. Therefore no further data were taken on this sample. 

Each sample, after it had been cooled to the temperature of liquid 
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Fig. 8 — Plot of v/(H — H.) for Sample 4 as a function of temperature. 
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Fig. 9 — Plot of wall velocity as a function of applied field at 77°K. Note that 
at the higher fields the velocity is no longer a linear function of the applied field. 


air once, gave the same value of v/(H — H,) as before. Samples 1 and 2, 
however, were cooled several times, and after the later runs they no 
longer did this. In both cases the value of v/(H — H.) as deduced from 
the ideal domain pattern and (1) was lower by about one-half; we inter- 
pret this to mean that the domain pattern in these samples was changed 
by the repeated thermal shock. Direct confirmation of this interpreta- 
tion by observation was not possible, however, for reasons indicated 
above in connection with domain pattern observations on these samples. 

In view of the fact that only one point is plotted beyond the knee of 
the curve in Fig. 8, it should be emphasized that continuous qualitative 
observations made while the sample was cooling showed that the change 
was continuous and monotonic. On Sample 1, furthermore, the data at 
201°K was somewhat down from the knee of the curve [v/(H — H,) = 
18000 cm/sec/oe] because of the fact that the knee occurred at higher 
temperature as mentioned above. 

Another feature of the data is indicated by Figs. 9 and 10. Data dis- 
cussed thus far have been taken at the lowest convenient velocities in 
order to minimize the possibility of wall distortion. However, when 
data were taken at higher fields, a non-linearity of the sort shown in 
Fig. 9 appeared. The average velocity increases more rapidly at the 
higher fields than our viscous damping coefficient would lead us to ex- 
pect. This effect was observed at all temperatures, but as comparison 
of Figs. 7 and 9 will show, it set in at lower velocities at low temperatures 
where the enlarged viscous damping appeared. Simultaneously with the 
appearance of this non-linearity in the apparent average velocity of the 
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wall, the voltage induced in the secondary winding during the pulse of 
constant applied field is seen from the oscilloscope trace to distort 
with time. Fig. 10 shows a series of traces observed at room temperature 
on Sample 4 which show this distortion increasing from (a) to (d) as the 
applied field is increased. At the highest fields the trace forms a peak 
which is almost triangular in shape. 

We shall discuss the theoretical implications of these data in the next 
two sections. 


3500 — 
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Fig. 10 — A series of oscilloscope traces showing the deviation of the wall ve- 
locity from a constant with time at the higher applied magnetic fields and there- 
fore at the higher velocities. These pictures were taken at room temperature, but 
similar phenomena occur at all other temperatures in the range of applied fields 
where the nonlinearity in velocity shown in Fig. 9 becomes apparent. The ve- 
locities, as shown, increase from a to d. 
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THEORY 


A theoretical analysis of the experimental results given in the last 
section divides itself rather naturally into three parts. First we char- 
acterize the data in terms of an equation of motion for unit area of 
domain wall. This means we determine the constants of motion (viscous 
resistance and coercive force in our case) of unit area of wall. Secondly, 
we show the relation between the viscous resistance of unit area of 
domain wall, and the constants which characterize the ferromagnetic 
material in general (saturation magnetization, crystal anisotropy, etc.). 
This is essentially an application of the recent work of Becker” and 
Kittel. Lastly, we calculate the magnitude of the damping from a 
relaxation mechanism which accounts for the low temperature effect 
shown in Fig. 8 at least qualitatively. 

Consider unit area of a 180° domain wall between two regions of 
saturated material. Such a system has an equation of motion for small 
amplitudes of the applied magnetic field H which may be written: 


mé + B2 + az = 2M.H, (3) 


where z is the displacement of the domain wall along its normal, m is 
its mass per unit area, 8 isa parameter measuring viscous resistance, and 
a is a stiffness parameter which has meaning only for small fields such 
as those used in initial permeability measurements. When fields larger 
than the coercive force are applied, as in our experiments, the term con- 
taining a disappears and the field effective in moving the wall is less 
than the applied field by an amount equal to the coercive force; this is 
shown by the data given previously in the section on results. This re- 
formulation of (3) is quite reasonable when one remembers the spikes 
which pull back on the wall, in the experiments of Williams and Shockley, 
for small wall motions and snap off entirely if the wall moves a large 
distance. Furthermore, since the velocity of the domain wall rises to 
its steady value in a negligible time in these experiments (Fig. 6) the 
initial term in (8) is also negligible. As these remarks indicate, under the 
conditions of the experiment in wall velocity, (2) takes the form: 


6z = 2M.(Happ — H-). (4) 


This relation obviously fits the data given previously. 
The second step in the theoretical analysis starts from an equation of 


_—_> 
motion for the magnetization M in a small volume which was first used 
by Landau and Lifshitz,° and takes advantage of more recent work of 
Becker” and Kittel.” If we consider a volume small enough so that the 
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magnetization in it is everywhere uniform even if we are inside the 
domain wall, our equation of motion is: 
dM _ 2 

+ yiM x H] — /M[if x (Mt HL (5) 


Here y is the gyromagnetic ratio (ge/2 mc), and X is a parameter, as- 
sumed to be characteristic of a given ferromagnetic material, which is 


—>> 
determined by the magnitude of the damping effects in the motion of M. 
The magnitude of the last term on the right in (5) is thus determined 
by the amount of the damping losses. 


> > 
The rate of dissipation of energy in the small volume is H-(dM/dt), 
— — 
where M is the magnetic moment of the volume, and H is the total 
> 
magnetic field in the volume. The value of H requires some discussion. 


Outside the wall H = Hy where i is the applied field, which is parallel 
to the wall. When the wall is moving, however, there is an additional 
field H, inside it. This field, which is normal to the wall, is a demag- 
netizing field which arises from the tendency of WM in the moving wall to 
have a component normal to the wall. This field has just such a value 
that the magnetization in the moving wall prowess about it with the 
Larmor frequency. Its value is: 


H, = —(v/y)(00/a2), (6) 


as Becker” has shown. Here » is the velocity of the wall, z is a distance 
coordinate normal to the wall, and 6 is the rotational angle of the mag- 
netization as we pass through the wall along z. Inside the wall, H, is 
much larger than ih , but in any case H =H. + H). From (5) we find: 


H-dM/dt 2 HZ, (7) 


as Kittel” first showed. In the theory of the domain wall it is shown 
that 00/dz in the wall is equal to the square root of the ratio of the in- 
crease in anisotropy energy as the magnetization turns away from the 
easy direction of magnetization, to the exchange energy constant. That 
is: 


= ([9(0) — g(00)]/A)*”. (8) 


The exchange energy constant A is defined by the following expression 
for the exchange energy per unit volume due to gradients in the direction 
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of magnetization: 
Exchange energy/unit vol = A[(Va)” -+ (Vaz)? + (Vas)"], (9) 


where a;, a2, and a; are the direction cosines of the magnetization. 
g(0) is the anisotropy energy: 


g(@) = Ky (ayas + as ats” -+ a3 01), (10) 


expressed in terms of 0, and g(6)) is the anisotropy energy along the 
direction of easy magnetization. Note that [g(@) — g(@)| is always posi- 
tive. K, is the first order anisotropy constant. 

If we use (6) and (8) in (7), and integrate over z along a cylinder of 
unit cross-section normal to the wall to get the rate of energy dissipa- 
tion for unit area of moving wall, we have: 

“ae = tPA" [ g@ — glo’? da = 2HoMw, (11 

[B Seae = ow'/'a") [lo — glo" ao = 2HeAt2, (11) 
where 6, and 2 are the angular positions of IZ on the two sides of the 
wall. 6. — 6, = a, of course, since we are considering a 180° domain 
wall. In order to obtain (11), we have used (8) to transform from in- 
tegration over z to integration over 6 as well as to evaluate (7). We set 
our result equal to 247,H (pressure on the wall) times v since this is the 
rate at which the wall, considered phenomenologically, does work. We 
may now write: 


2 My? A 1/2 


rT 

12 

rN | [9(0) — g(6)"" de 

This is the desired relation between wall velocity and applied field which 
is to be compared with (4). In this way we find: 


G2 
B= Q/7'A™) | [9(0) — g(oo)}? aa. (13) 


We have thus shown the relation between the wall parameter 8 and 
the parameter \ which measures in general the losses associated with 
motions of the magnetization. 

The third part of our theoretical analysis is concerned with a cal- 
culation of the damping parameter \, or rather the relation between 
v and Z itself, from an explicit physical mechanism. Such a calculation 
has not been made in the past since the appropriate mechanism on 
which to base it has remained obscure. Verwey and his co-workers” 
have explained the well known transition at about 115°K in Fe,O, as 
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an order-disorder transition in the arrangement of the divalent and 
trivalent iron ions. M. Fine of Bell Telephone Laboratories has found a 
remnant of this transition in crystals of the same composition as those 
studied in the present research, by means of ultrasonic measurements 
of elastic constants.’° We propose that the mechanism which causes the 
sharp rise in the damping of domain wall motion at low temperatures is a 
relaxation associated with this transition. Wijn and van der Heide” 
have explained in this way observations of theirs of losses associated 
with initial permeability at low frequencies in certain other polycrystal- 
line ferrites. The time associated with this relaxation should be short 
because this rearrangement of ions involves only the motion of electrons 
from one site to another.” It should be of the order of the relaxation 
time associated with the electrical conductivity of Fe;0,. Snoek” has 
suggested some time ago that losses in the ferrites were due to an after- 
effect (relaxation) which, because of the short time constant involved, 
must be associated with electron migrations. 

It is extremely useful to compare our data with a theory of the damp- 
ing based on the above relaxation mechanism no matter what assump- 
tions we make in detail about what it is that relaxes. We shall see that 
we are led quite generally to the result that v/H ~ 1/7 where 7 is the 
relaxation time for the process. However, in order to perform this cal- 
culation explicitly we must make more detailed assumptions about 
exactly what quantity relaxes with the relaxation time 7. Changes in the 
direction of the magnetization cause changes in stress in the sample 
because of magnetostriction. One possible assumption is that the re- 
sulting strain would lag behind this stress and mechanical energy would 
be dissipated in the crystal. This mechanism, however, cannot act in 
our case. The magnetization in the two domains on each side of the wall 
points in opposite directions, but causes the same strain in both, and 
they have such a large stiffness that the thin region occupied by the 
domain wall assumes this same strain even though the direction of the 
magnetization is different there. Thus regardless of the stresses produced 
in the wall, the strains remain the same, inside and outside the wall, 
whether it is moving or not; under these conditions no work is done on 
the lattice, and no energy can be lost in this way by the moving wall. A 
calculation has been made by the author on the assumption that it is 
the magnetization itself which relaxes with relaxation time 7. This 
assumption leads to the result that wall velocity is not linearly dependent 
upon (H — H,); it is therefore not correct, since the data shows such a 
linear dependence. A similar result is to be expected if we assume that 
the dielectric polarization relaxes. 
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What seems at present likely to be the approximate nature of the 
mechanism, and what we will assume is the nature of the mechanism 
for the purposes of an illustrative calculation is as follows. As the domain 


=> 
wall passes a point in space, and the direction of Af changes, the electrons 
on the divalent and trivalent iron ions tend to rearrange themselves so 


as to minimize the magnetocrystalline anistropy energy.” If i changes 
slowly this anistropy energy is near the minimum possible value (the 
reversible value) at all times, and the process is.almost isothermal. As a 
result of the fact that the process deviates irreversibly from equilibrium, 
however, net work is done in bringing about the change. If, on the other 


—> 
hand, the direction of Af changes so suddenly that the electrons have no 
time to rearrange, the process is adiabatic, and the magnetocrystalline 


anistropy energy varies more widely with the angular position of M1 : 

Since our data is taken at low velocities and extrapolated to zero 
velocity, it seems most appropriate for us to make a calculation of the 
losses on the assumption that as we increase the velocity of the wall 
we are deviating from the isothermal condition. Let us define as a ther- 
modynamical system the part of the magnetic lattice which lies in a 
small volume fixed in space. This volume is a sheet of unit cross-section 
in which the magnetization is uniform and which is part of the cylinder 
of unit cross-section normal to the wall mentioned in connection with 
(11). From the first law of thermodynamics, as the wall passes the small 
volume, we have: 


dw = dU — dQ = dg, (14) 


where dQ is heat added to the system, dU is a change in internal energy, 
dw is work done on the system, and g is the anisotropy energy as- 
sociated with our rearranging electrons. Note that g, is a term in the 
free energy of the magnetic lattice. The free energy includes other 
terms in addition to g, and indeed there is in general another term in 
the magnetocrystalline anistropy energy; we will not consider any of 
these, however, since they also integrate to zero as the domain wall 
passes our small volume. Similarly, there are many contributions to 
dw, -but we will concern ourselves only with the term which does net 
work on our system, the pressure on the wall times its velocity. If we 
now consider changes in our system with time we have from (14): 


Cea 
di dt’ 


The rate of doing work on unit area of the wall moving along our 


(15) 
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cylinder, which is the integral of dw/dt over the cylinder, is 2/M,Hw 
as in (11). The rate of dissipation of energy in this section of the wall 
is somewhat more difficult to calculate. It is the sum of contributions 
from all the small volumes along the cylinder. In order to calculate the 
contribution from the small volume we are considering, we first note 
that if the process is reversible, 


_ dg 

dg = 6 dé, (16) 
where @ is the angle of rotation of M. . Physically, the factor dg/d@ repre- 
sents a torque on the magnetization. We will assume that (16) holds 
even when the domain wall passes our system at a finite velocity and 
we have departed slightly from isothermal equilibrium. The field H, 
defined in (6) transmits this torque to make the magnetization rotate 
through the angle 6, but we will not go further into the details of this 
process. 

As the domain wall goes by, dg/dé, which we will abbreviate as g’, 
changes. If the process were reversible, g’ would be zero at all times and 
the torque on the magnetization would always have its equilibrium 
value. Actually, as we deviate more and more from the reversible process 
by moving the wall faster, the electrons are no longer able to rearrange 
fast enough, and g’ deviates from zero while continually relaxing toward 
it. We must form our analysis in such a way that a maximum is estab- 
lished for g’, since even if the magnetization moves infinitely rapidly, 
g’ does not become infinite. Since the electrons minimize their free energy, 
it must be positive, and we write g = g1,, — gi so that g’ = gi, — 91, 
where g; relaxes toward g;,,. Note that the torque relaxes downward 
from a value, gj, , associated with the adiabatic anisotropy energy (fast 
motions of the magnetization) to a value zero, associated with the iso- 
thermal anisotropy energy (slow motions of the magnetization). We 
may now write: 

dg: _ Jie — 1 (17) 

dé oor - 
Physically, this relation assumes that the torque on the magnetization 
relaxes in the same way if gio is a continuous function of time as it 
would if gj.. and g; differed but gj. was constant. Here 7 is the relaxation 
time associated with the rearrangement of divalent and trivalent iron 
ions. If we write (17) in the form: 

dg 


fi _ Sis 


/ 
1 
T T 
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we see at once that the solution for g; as a function of ¢ as the domain 
wall passes the position of our small volume is: 


ds ayy : 
gl) ==" | gig(et dt (19) 
We have ignored the solution of the homogeneous equation, as it con- 
tains a factor ¢””. 

In order to evaluate (15) for our small volume, we need a more ex- 
plicit value for g;. From the Fourier transformation we may write: 


giolt) = | Gralwe" de (20) 
Consequently: 


gp er i! | GCaolewett des dt. (21) 
T — oo 


This can be integrated over ¢. If we do this, then bring the factor e‘” 
out from under the integral sign, we find: 


, = 4 G0 (wo) twt 
gx(t) — i, ear ee du. (22) 


Now the frequencies for which we get a contribution to the integral in 
(22) have an upper bound determined by the velocity and the thickness 
of the domain wall, of course. The faster the domain wall moves, the 
higher the upper bound on the frequency. We have been careful to make 
our measurements near zero wall velocity, and the data clearly extra- 
polates linearly to zero. This was done in order to minimize the possibility 
of wall distortion, but it also is consistent with our assumption of al- 
most isothermal conditions. We are therefore justified in assuming that 
wr <1 for all the frequencies in Gj.,.(@), and we can expand the factor 
1/(1 + dwr) in (22). We obtain: 


gi(t) = [. [1 — (wr) + (iwr)” + -:- 1Giuloje du, (23) 


and if we form the derivatives of (20) and compare them with the terms 
in (23), we see that: 
di (t) oh = @ Ji0(t) 


dt apt “ep 





gilt) = Giolt) — 7 
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Since the value of 7 is very much less than one, the series in (24) con- 
verges quite rapidly at low domain wall velocities, and we only keep 
the first two terms. 

Until now, we have been concerned with a small volume of the mag- 
netic lattice fixed in space. As in deriving (11), if we wish to calculate the 
rate of loss of energy for unit area of the domain wall moving with con- 
stant velocity, we must integrate dgi/dt over a cylinder of unit cross- 
section normal to the wall. g:,, is now a function of (¢ — z/v) where v is 
the velocity of the domain wall, and z is a coordinate normal to the wall. 
If we form this integral, and use (16) and (24), we find: 


[Mae = [| aha — 27) — Be @ - af) [MET ae (05) 
We note that: 


dgialt— 2/0) _ _, dnlt — 2/2) 
dt dz j 
do(t — z/v) =e dot — z/v) 
dt dz ; 


For the first term in (25) we find, using (26): 


(26) 


[ Jig(t — 2/v) oe el) dz = —v [ " Jialt — 2/v)d0=0, (27) 


since 91, is the same on both sides of the domain wall. Now if we use (26) 
in the second term of (25), and remember the result in (27), we have: 


Ot ay = — gy? [Wiel — 2/0) dO(t — 2/2) 2 


3d a dz dz (28) 


We may, without loss of generality, evaluate the intergral in (28) 
at t = 0. It is therefore the integral of a function of z over z. In order to 
evaluate it, we wish to express g;,, aS a function of 0. We have: 





dg. _ a f” dia, (2) 
[% dz = TU | rT dz dz, (29) 
°2 a1, d0 
sp | Sa Yio WW 
= —7V » de? de dé, (30) 


where we have replaced g;,, by dg:,,/d0. We note now, from the relation 
between g and gi, that dgi/dt = —dg/dt. We therefore set the right hand 
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side of (80) with reversed sign equal to: 


* dw 
¥ or dz = 2M.Hw, 
as mentioned after (15) to get the relation between the velocity of the 
domain wall and the applied field. As we shall see later, the right hand 
side of (80) is positive if we use a g;,, associated with a positive contribu- 
tion to the anisotropy energy. We find: 





1 2M, 
ae ie Bin dB ag” (31) 
6, de dz 


The derivatives of @ with respect to z in (29), (80), and (31) are to be 
evaluated by means of (8) and (10). In using these equations, of course, 
we are assuming that the wall is moving slowly enough so that its shape 
remains that of the wall at rest. 

Equation (81) shows that v ~ H/r as we mentioned earlier. Inspec- 
tion of Fig. 8 shows that this relationship explains very satisfactorily 
the sharp drop in v/(H — H,) at low temperatures if we remember that 
7 depends on temperature as follows: 


T = To eo. (32) 


where ¢ is an activation energy. 
Finally, the right hand side of (80) may be set equal to: 


(i -dit /at) de, 
as calculated from the Landau-Lifshitz equation, see (11), to obtain a 
value for \. We find: 


2 41/2 2 gy, 0 
x 41 de dz 


Spies eb 
I [9() — g(0)]"" do 


dé 





(33) 


DISCUSSION 


It is clear that (4) fits our experimental data at each temperature. 
By fitting our data to this expression we obtain values for 6, the para- 
meter characteristic of the material which measures the damping of the 
wall. Values of 6 obtained in this way are given in Table I. A more 
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correct value of M, for Fe;0, at room temperature (475 c.g.s. units) has 
been used in calculating the 6 given in Table I than has been used in 
previous calculations. 

In Table I we give in addition to values of 8, values of \ at room tem- 
perature obtained from (13). The value for Fe30, is calculated from the 
new value of 6 and also corrected for the error mentioned in Reference 
14. Values given in References 7 and 8 are somewhat in error. The values 
of 8 and 2 given for Fe30, in Table I are the remainders after the contri- 
bution due to eddy currents has been subtracted out by means of the 
low field calculation of Williams, Shockley, and Kaittel,? [see their Equa- 
tion (11)]. 


TaBLE I— Room TEMPERATURE Data on Fe;04 AND 
(NiO)o.75 (FeO)o.25 Fe.03 








(c.g.s. units) 











B (corrected for eddy domain d ferromagnetic 
currents) wall resonance 
BiG gO fic beecy hc eee ease atta ee 0.44 5.5 X 108 9 X 108 
(NiO) 0.75 (FeO)o.25 FesOs....... 0.023 4.6 X 107 10 X 107 


TaBie I[— Dara For (NiO)o.75 (FeO)o.25 Fe2O3* 











T(°R) Ms (c.g.s. units) Ki (ergs/cc) ie a aber ay 
77 341 —8.1 X 101 158 6.3 X 109 
201 335 —5.4 X 103 26150 4.4 & 107 
300 322.5 —3.8 X 10! 28500 4.6 X 10° 
363 309 —3.1 X 101 32500-85300 4.0-4.3 X 10’ 
400 298 —2.8 X 104 31750 4.5 X 107 
445 281 —2.4 X 104 36850 3.9 X 107 














* The value of M, at 300°K is measured from the hysteresis loop of Fig. 3. 
Other values were obtained by assuming that Af, varied with 7 in the same way 
as observed by Weiss and Forrer! in Fe;Q, . 


‘The evaluation of \ from (13) is done as follows. Since the wall is in a 
(110) plane, we find from (10): 





0(0) — g(a) = #2! (//3 — V8 sin’ 0) (34) 


med 
where @ is the angle between M and the [100] direction which lies in the 
plane of the domain wall. @ is cos "(1/+/3) on one side of the wall, and 
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a + cos” 1(/+/3) on the other. Then, 
— 
cos—1(1/+/3) 


In performing this integration, care must be taken to use the positive 
value of the square root over the whole interval. It should be noted that 
in using (8) and (10) to evaluate (13) we are assuming that the wall is 
moving slowly enough so that its shape is the same as that of the wall 
at rest. A is best evaluated from a fundamental relation derived by 
Herring and Kittel” between A and the Bloch constant: 


[9(0) — g(60)}'"? da = 0.92 | Ki]. (35) 


A = [S/o]"[k/13.3C7"), (36) 


where & is Boltzmann’s constant, C is Bloch’s constant as used in the 
relation M, = M)(1 — CT*”), So is the atomic spin, and Q is the atomic 
volume. (So/Q) is equal to the saturation magnetization at O°K divided 
by twice the Bohr magneton. 

For Fe;0,, we find C = 4 X 10° by fitting the Bloch T*” law to 
the saturation magnetization measurements of Weiss and Forrer.” 
From (86), assuming J/7, at O°K is 505 ¢.g.s. units, we then find A = 
1.24 X 10°°. Furthermore, K; = —1.1 X 10° as given by Bickford,” and 


y = (1.76 X 10")g/2 = 1.865 X 10’, 


where we have used Bickford’s” value (2.12) of g.” Now from (13) we 
find 


\ = 5.5 xX 10° 


in Fe;0, at room temperature. 

For (NiO)o.75(FeO)o.2sFe2O3:, we assume that M/,, while different 
from that for Fe;O, , varies in the same way with temperature so that 
C = 4 X 10°. From our measurement of I, at room temperature 
(322.5 c.g.s. units) and this assumption about the variation of M/, with 
T, we find that M, at O°K is 342 c.g.s. units. This leads by (86) to A = 
1.09 X 107°. Ferromagnetic resonance experiments” done by W. A. Yager 
and F. R. Merritt in collaboration with the author on spherical single 
crystals of the same ferrite material as that used in the present research 
give a g value of 2.14 at all the temperatures mentioned in Table IT 
except 77°KK, where g = 2.19. These values are used in determining the 
value of y[= 1.76 X 10’ g/2] in (13). The anisotropy values in Table IT 
are also taken from the results of these ferromagnetic resonance experi- 
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ments. The Ki at room temperature is —3.8 X 10°. The room tempera- 
ture value of \ domain wall given in Table I and the values at various 
temperatures given in Table II are obtained from (13) using these data. 
An independent value of \ may be obtained from the ferromagnetic 
resonance line width. The relation between the observed line width 
2AH and X has been given elsewhere.” Sample shape enters this relation, 
but not in a critical way, and we therefore ignore it except as it affects 
the value of the dc magnetic field at resonance, H,.;. The relation is: 


\ = AHyYM./Hees - (37) 


From Bickford’s” data on line width and (37) we find \ = 9 X 10° for 
Fe,;0, at room temperature. From the ferromagnetic resonance data 
reported elsewhere” on (NiO)o.75(FeO)o.25Fe203, the material used in 
the present research, we find \ = 10 X 10’ at room temperature. Table I 
compares the room temperature values of \ obtained in the two ways on 
the two materials. 

The differences between the domain wall experiments and the fer- 
romagnetic resonance experiments lead to quite different behavior of 
» in the two cases at low temperatures. These differences can be under- 
stood in terms of the frequency dependence of \ as given by an extension 
of the relaxation theory given in the third part of the theoretical dis- 
cussion. A discussion of these relationships must await the detailed 
report on the ferromagnetic resonance results which is now in prepara- 
tion, where such an extension will be given. As Table I shows, the room 
temperature values obtained in the two ways are of the same order of 
magnitude. 

Let us now turn to a discussion of the relation between (31) and (83) 
and the data shown in Fig. 8. Qualitatively, of course, the 1/7 factor in 
v/Hy which (81) reveals, taken together with (82) explains most satis- 
factorily the sharp increase in viscous damping of the domain wall at 
low temperatures. Furthermore, it seems quite possible, although the 
author has not investigated it, that the higher order terms in (24) account 
for the nonlinearities in Figs. 9 and 10. 

It should be mentioned that an increase in relaxation time at low 
temperatures which is consistent with (32) has been deduced by Bloem- 
bergen and Wang” and Healy” from ferromagnetic resonance data 
taken by them. 

Quantitatively, we have inadequate data for a satisfactory compari- 
son between Fig. 8 and (31), and the assumptions of the theory should 
perhaps be investigated further before any such comparison is taken 
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seriously. Nevertheless, plausible assumptions can be made which make 
an instructive comparison possible. g1,, in (31) is equal to the difference 
in the variation of the anisotropy energy when measured adiabatically 
and when measured isothermally. At this stage of our knowledge, many 
assumptions which still retain the symmetry of the crystal are possible 
concerning the form of q,,(@). For our present purposes we will assume 
that it is given to within a constant by (34), but we must introduce a 
minus sign to take account of the fact that the rearranging electrons 
must have a positive anisotropy energy associated with them. Since we 
differentiate before substituting in (81), we do not need to subtract out 
the constant. The amplitude of qi, is not given by | Ki| of course; we 
will call this amplitude | K, |. When we introduce the minus sign into 
(34), and calculate the integral in (31), we find: 


1 40M, q 
v=- Cll : 
ane res (38) 
. A 


This relation points up the fact that the mechanism we are discussing 
here is characterized by two parameters, an amplitude factor | K, | and 
a time 7. Our experiment tells us nothing about | K,|, so we will ar- 
bitrarily assume it is about 14 of | Ki| at room temperature, or 20000 
ergs/cc. | K, | can be measured by comparing values of | Ki | determined 
isothermally, say by measuring the torque on a disc, and values of 
| K,| ‘determined adiabatically, say by ferromagnetic resonance ex- 
periments, but this has not been done as yet. To determine 7 at 77°K, 
assume that the maximum frequency of rotation of dipoles in the wall is 
such that wr 2 0.1 when the non-linearity shown in Fig. 9 first occurs 
(at v = 150 cm/sec). We calculate the thickness of the wall to be 4 X 10~° 
cm from standard formulae (see Kittel’s review article, Reference 14) 
and find 7 20.5 X 10°°. These values of 7 and | K,|, together with 
| | from Table II at 77°K and A as calculated above, give v/Hy = 40 
when inserted into (88). This is to be compared with a measured 
v/(H — H.) of 160 at 77°K. In view of the preliminary nature of (81), 
and the arbitrariness of some of our assumptions, this check is quite sat- 
isfactory. It would be naive to expect the theoretical result to be closer 
than an order of magnitude to the experimental one. Inspection of Fig. 
8 suggests that the mechanism which gives the sharp increase in damp- 
ing at 77°K is submerged at room temperature in the effects of other 
mechanisms, perhaps other electronic rearrangements, perhaps exchange 
effects of the sort recently suggested in metals by Rado.” Equation (39) 
confirms this expectation when we use 7 & 10°”, as determined either 


1052. THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1954 


from ferromagnetic resonance data” or from the conductivity of Fe;O,. 
The above method of calculating 7 from domain wall data is less sat- 
isfactory at room temperature. 

The idea of associating losses in ferrites such as this one with 
changes of order in the divalent and trivalent iron ions explains the fact 
that the damping observed in domain walls at room temperatuie in 
Ie,O; ,’ where none of the divalent iron is replaced by nickel, is larger 
than that observed in (NiO)o.75(FeO)o.25F e203. Finally, as Wijn and van 
der Heide” have pointed out, and as (31) and (33) show more analyti- 
cally, this mechanism is a very satisfactory explanation of the sharp 
change in domain wall relaxation frequency observed by Galt, Matthias, 
and Remceika.” 
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Negative Impedance Telephone Repeaters 


By J. L. MERRILL, JR., A. F. ROSE, and J. O. SMETHURST 
(Manuscript received June 8, 1954) 


In the exchange telephone plant, speech is transmitted largely at voice 
frequencies over a single pair of wires that carries conversation in both direc- 
tions. Until recently, adequate transmission was assured through a suitable 
choice of coil loading and conductor size. The need for voice frequency amplt- 
fiers had long been recognized, but wide use of the conventional hybrid type 
of repeater with separate amplifiers for conversation in the two directions 
had not been economic. In 1948, however, a new type of repeater was intro- 
duced that used the same source of gain for transmission in both directions, 
and did not requtre costly filters. This repeater is known as the 1 repeater. 
It operates as a negative impedance in the line. 

This paper describes the wide field of usefulness of two new types of 
negative umpedance repeaters, the H2 and ES. Sufficient theory is given to 
show the advantages of using two types of negative impedance in a line, the 
series type and the shunt type, and how, by the addition of. the shunt type 
to the series type, it 1s possible to insert gain without serious reactions on 
line wmpedance. The new equipment is described and maintenance fea- 
tures, together with methods of testing this relatively new type of repeater, 
are discussed. 


APPLICATION IN THE BELL SYSTEM 
INTRODUCTION 


The application of the telephone repeater, the development of which 
made countrywide telephone service practicable, had been confined 
largely to toll plant from the year 1915 when the transcontinental line 
was first established, until a few years ago. About 1948 the negative 
impedance repeater’ was developed and placed in production. This 
repeater operates on the principle of inserting negative resistance (and, 
if desired, negative inductance or capacitance) in series with the line, 
thus reducing the overall impedance and increasing the current in the 
line. This results in transmission gain in the same sense as that resulting 
from a repeater of the conventional type. This principle and the package 
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nature of the assembly resulted in a telephone repeater so low in cost 
and so simple in application and installation that it has found extensive 
use primarily in the local telephone plant. 

In the period from 1948 to the present time, over 50,000 series-type 
negative impedance repeaters were manufactured and incorporated in 
the Bell System telephone plant. These repeaters have been used largely 
on intraexchange trunks and on trunks extending from the exchange 
areas to near-by smaller towns. Such installations have been very effec- 
tive in improving the transmission on short haul calls and in many cases 
have also reduced trunk costs by permitting the use of smaller and 
cheaper conductors. They are usually operated at gains that reduce the 
nonrepeatered trunk loss by more than half. 

Negative impedance repeaters are especially suited to exchange trunk 
use, not only because of their simplicity and ease of maintenance, but 
also because unlike earlier types of repeaters, they preserve the de 
continuity of the circuits on which they are installed. This latter fea- 
ture means that they do not interfere appreciably with the signaling 
methods ordinarily used on such trunks. Also, they have the added 
advantage that, in the event of tube failure, the circuit still functions 
but with its loss substantially raised until the defective tube has been 
replaced. 


APPLICATION OF SERIES-TYPE NEGATIVE IMPEDANCE REPEATERS 


In a multioffice exchange area, trunks can generally be grouped in 
three classes. Largest in number are those known as interoffice trunks 
which extend directly from one local operating center to another local 
center in the same exchange operating area. In some cases direct trunks 
between two centers can not be justified and each office has trunks to a 
central location known as a tandem center where a through connection 
is made when required. This second class of trunks is known as ‘‘tandem”’ 
trunks. Frequently trunks of this type are provided to supplement the 
direct trunks and thus give the added advantage of alternate routing. 
A third type of trunk extends between the local office and the toll office 
and therefore always forms part of a toll connection to the local office. 
These are known as toll connecting trunks. In single office areas the toll 
connecting trunks to neighboring small towns are known as ‘‘tributary”’ 
trunks. 

In some large multioffice exchange areas there are also trunks between 
the toll office and tandem centers which may be some distance from the 
toll office and through which local offices reach the toll office for long 
distance calls. These trunks are also called tandem trunks but are more 
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in the nature of intertoll links and are required to operate at very low 
loss. similar to an intertoll circuit.’ 

The most extensive use of negative impedance repeaters has been in 
connection with the interoffice trunks. In this role they have been ex- 
tremely useful and have contributed greatly to transmission improve- 
ment in the last few years. It would also be desirable to utilize this new 
and effective tool in reducing the losses of toll connecting and tandem 
trunks. There have been many cases of this type where improvement 
has been effected, but a very wide use has been prevented by their 
effect on return loss. 

The insertion of a series negative impedance in an otherwise uniform 
loaded circuit introduces an impedance discontinuity which means that 
a substantial amount of energy is returned to the sending end of the _ 
circuit as “echo.” As the introduction of negative impedance is the 
method of obtaining gain from the series repeater it follows that the 
magnitude of the impedance discontinuity and therefore the “echo,” 
is proportional to the gain of the repeater. The effect of echo on ease of 
conversation depends both on its magnitude and on the amount of delay® 
before the echo reaches the talker’s ear. When the delay is small as in 
the case of an interoffice trunk a large amount can be tolerated, but in 
the case of a toll connecting trunk which becomes a part of a toll con- 
nection the delay may be large and only a small amount of echo can be 
permitted. This limits the amount of gain from a series repeater when 
it is used in toll connecting trunks. 

As a result of this feature of the series negative impedance repeater, 
most cases requiring substantial gains in toll connecting or tributary 
trunks usually had to be handled by the older and more expensive hybrid- 
coil type repeater. With this latter repeater, it is practicable to introduce 
gain without serious reflection effects if proper attention is given to 
network design and to uniformity of construction of the line itself. 
However, the cost of the more complicated repeater was relatively 
high so that the new development described herein was undertaken to 
meet the indicated need, i.e., a repeater embodying the same desirable 
features of simplicity of design but still approaching in performance 
that of the hybrid-coil repeater in its effect on return loss. 

The new repeater (E23) consists essentially of the earlier series re- 
peater with the addition of a shunt negative impedance element. This 
combination has approximately as small an effect on return loss as the 
hybrid-coil repeater and will give about the same gain when used under 
similar line conditions. As a result, the field of use of the negative im- 
pedance type repeater will be greatly extended, especially in the toll 
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connecting plant where, from a practical standpoint, certain features of 
the hybrid-coil repeater are not needed. For example, more than 10 
db gain is rarely required, and unequal gains in opposite directions would 
not ordinarily be useful. Besides, methods of signaling requiring a 4-wire 
split of the repeater are not utilized on exchange area trunks. 

In addition to its important new characteristics, the series-shunt 
combination retains the desirable features of the earlier series repeater, 
i.e., it preserves de continuity, it is simple in design and it is easy to 
install and maintain. If properly engineered and installed on circuits of 
reasonably uniform impedance, the new repeaters can be operated in 
tandem, do not have serious reactions on return loss, and are capable of 
reducing the losses of the trunks to about the same values as hybrid-type 
repeaters. 


SCOPE OF APPLICATION OF THE NEW REPEATER 


With the reduction of the return loss reaction of the earlier form of 
negative impedance repeater, the extent of application of the new E23 
repeater becomes largely a matter of economics. In the telephone trunk 
plant, there are generally three gauges of conductors available for trunk 
use, namely, 19, 22 and 24. The larger gauges, utilizing more copper and - 
requiring more conduit space, naturally cost more than the smaller 
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Fig. 1 — Illustrative field of use; interoffice trunks. 
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Fig. 2 — Illustrative field of use; toll connecting or tandem trunks. 


gauge. In general, it is the practice to select the smallest gauge which 
will give the required transmission loss. For short distances, the dif- 
ferential in cost between trunks of different gauges is not large but as 
trunk length increases it will be found that a smaller gauge trunk with 
a repeater will cost less than one of larger gauge without a repeater. 
For example, at 10 miles a 22-gauge circuit with a repeater will cost 
substantially less than a 19-gauge circuit with no repeater. Also, as 
transmission objectives are improved, there may be cases where a smaller 
gauge trunk with two repeaters will be cheaper than a larger gauge with 
one repeater. Furthermore, in cases where the costs are about equal or 
even where the larger gauge is slightly cheaper, it will be found that 
lower losses can be obtained with repeatered circuits and the engineer- 
ing choice would accordingly favor the smaller gauge. 

Construction costs differ considerably depending on local conditions; 
hence the differential cost between the various gauges of conductors will 
differ correspondingly. However, for illustrative purposes Figs. 1 and 
2 have been prepared which utilize average Bell System costs for both 
outside plant and telephone repeaters. In the case of the E23 repeaters, 
the experience so far has been somewhat limited so there is a greater 
degree of uncertainty as to the actual costs than for the older type re- 
peater or the outside plant construction. 

Fig. 1 shows the field of use for various conductor gauges when used 
in interoffice trunks where an over-all transmission objective of 6 db for 
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the trunks has been assumed. In considering this chart it will be seen 
that non-loaded 24-gauge conductors can be utilized for the shorter 
distances up to about two miles. The next step is to apply loading, which 
is indicated in the chart by the symbol H88 meaning 88 millihenry coils 
at 6,000 foot spacing. After reaching the limit of about 7 miles without 
repeaters on 22-gauge, the differential between 22- and 24-gauge is 
more than enough to pay for a repeater, the simple scries repeater (22) 
being used for distances up to 10 miles and then the improved E23 re- 
peater extending the use of 24-gauge to about 14 miles. Beyond this 
point the most economical combination is indicated. 

The general effect of the introduction of the negative impedance re- 
peater is to shift the average gauge distribution so that more small gauge 
cable can be economically utilized, accompanied in most cases by im- 
proved transmission. It will be noted from the figure that no 19-gauge 
cable will need to be added in the future to care for interoffice trunks 
up to distances as great as 25 miles between offices. Of course, in cases 
where 19-gauge is already available in plant it can be used to advantage 
despite the fact that for new construction a smaller gauge with the nega- 
tive impedance repeaters would be cheaper. 

It has also been assumed in making up Fig. 1 that supervision or 
pulsing requirements will not limit the use of the smaller gauges. In a 
specific case where some of the older type central offices are involved, 
signaling may have an important bearing and substantially distort the 
economic ranges indicated on the chart. 

Fig. 2 illustrates the field of use of toll connecting or tandem trunks. 
As a toll connecting trunk forms part of a multilink connection, and 
since there are always at least two (and often more) trunks in series on 
connections involving these trunks, the Bell System companies find it 
economical to plan for a maximum loss of 3 db for this type of trunk, 
Likewise for the tandem trunks, where two in series may be used in 
place of an interoffice trunk, 3 db is considered a reasonable objective. 

It will be noted again by referring to this chart that 19-gauge has very 
little future field of use and in all cases the new ‘“‘series-shunt”’ repeater 
will be utilized rather than the earlier series type because of return loss 
considerations and also because of the greater gain required to reduce 
the trunk losses to the desired values. In the tandem and toll connecting 
case, signaling may again be a distorting factor though not to so great 
an extent as in the direct interoffice trunk case. To this extent, however, 
the curves are theoretical, as they have been made up without regard 
to this limitation which may apply in a few practical cases. 
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SPECIFIC APPLICATIONS 


Many installations of the new repeaters have been engineered for 
completion in 1954. In all cases economic studies were made and re- 
sults of these studies broadly confirmed the indications of the two charts 
discussed above. To bring out more clearly the effectiveness of the 23 
repeater, it may be worth while to consider a few specific cases involving 
installations of the new repeaters being made this year. 

The first case shown on Fig. 3 is in the area of the Pacific Telephone 
and Telegraph Company. The figure shows the two alternatives that 
were considered for providing additional tandem trunks between San 
Francisco proper and the East Bay Area, needed this year because of 
an extension of customer toll dialing arrangements. 

The engineering study for this project was somewhat more com- 
plicated than would ordinarily be the case because two possible routes 
of unequal length and with slightly different amounts of submarine 
cable were involved. The present route which touches Yerba Buena 
Island is subject to some hazard from dragging ship anchors but the 
circuit relief would have been cheaper here than on the other route 
shown, were it not for the new repeaters. The second route is consider- 
ably less hazardous and, in addition, is sufficiently removed from the 
present one to provide increased reliability under disaster conditions. 
Without the repeaters, however, the second route would have been 
somewhat impracticable, since it does not allow easy installation and 
access to required loading points between the two shore lines. However, 
with two of the new repeaters on each pair of conductors, nonloaded 
22-gauge cable gives substantially the same effective transmission loss 
as loaded 19-gauge conductors on the shorter route and over a future 
period will involve less annual cost per circuit. Furthermore, the use of 
22-gauge cable permits more than twice as many pairs to be included in 
the same size sheath in the expensive submarine section. 

The initial installation of the new type repeaters on this cable will 
total about 1,300, divided between the Main Office in San Francisco 
and the Main Office on the East Bay side. Individual repeater gains 
are 5 to 7 db at 1,000 cycles but the repeater gain characteristic is shaped 
to offset the increasing cable losses at the higher frequencies so that the 
over-all circuit has relatively uniform transmission in the voice range. 

The second case is one where the need for the repeaters results from 
the complex toll switching system at Chicago, Illinois. Here it has been 
found desirable, because of the great volume of toll traffic, to have several 
toll offices scattered throughout the city and suburbs. In general, toll 
calls coming into these offices from other cities are completed over toll 
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Fig. 3 — Typical E23 repeater application, San Franciso — East Bay tandem trunks. 
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connecting trunks direct to the called subscribers’ central offices. It is 
not economical, however, to provide enough such trunks to handle peak 
loads. When all of these direct trunks are busy, an incoming toll call is 
switched to the subscriber’s central office via a tandem office serving 
his general area. Since the losses of trunks between the tandem and 
local offices are of the same order as those of the direct trunks from toll 
to local offices, it would be desirable to operate the trunks between the 
toll and tandem offices at losses close to 0 db if this were practicable. In 
this way the same transmission objectives would be met on both rout- 
ings and no contrast in transmission would be evident to the same 
subscribers on calls completed over the different routes at different times. 
Vig. 4 illustrates a specific case of 36 trunks between toll office No. 3 
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Fig. 4 — Typical E23 repeater application; Chicago toll tandem trunks. 
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and Belle Plaine tandem, a distance of about 16 miles. Without repeaters 
these trunks would have been operated at a loss of about 10 db but with 
the repeaters they are operated at about 2 db. The only alternative to 
using | repeaters in this specific case would have been to use the more 
expensive hybrid-type repeaters. 

In some cases similar to this one, it was found practicable to temporize 
by installing the series repeater first and operating it at limited gain 
until the shunt element became available. Later, when the production of 
the shunt element was started, the additional, units were added and 
full advantage of the new series-shunt design utilized in reducing the 
equivalent of the trunks to the lowest value permitted by echo return 
loss considerations. 

The third example shown on Fig. 5 is in the city of Pittsburgh between 
Churchill tandem and the town of New Kensington, approximately 17 
miles. Here the transmission on existing 19-gauge loaded cable without 
a repeater would have been about 8 db. The repeaters reduce this figure 
to between 2 and 3 db which should be satisfactory in this case. It should 
be noted, however, that 22-gauge with two repeaters would also provide 
a low enough equivalent, and should it become necessary to supple- 
ment the present cable at some future date, there will be an opportunity 
for further savings from the E23 repeaters. 
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Fig. 5 — Typical E23 repeater application; Pittsburgh toll connecting trunks. 
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THEORY 
NEGATIVE IMPEDANCE CONCEPT 


Both the E2 and E38 repeater elements contain an amplifier having 
multiple feedback paths. The operation of an amplifier circuit of this 
type can be explained by classical feedback theory. However, experience 
with the E1 repeater over the past four years has shown the value of 
using a negative impedance concept in engineering such a device. Hence, 
in the explanation of operation given here, the repeater units will be 
treated simply as two-terminal networks which have negative impedance 
‘inputs over the frequency band of interest. The effect of introducing 
these impedances into telephone circuits can then be computed by the 
same simple network theory used to determine the effects of passive 
impedances. 


THE E2 REPEATER UNIT 


The E2 repeater is essentially a two-terminal network the impedance 
of which has a magnitude | Z | and a negative phase angle that can vary 
with increasing frequency from minus 90 degrees, or less, through minus 
180 degrees to at least minus 270 degrees. This type of negative im- 
pedance is shown in the diagram of Fig. 6(a). It has been known for 
many years as the series type because it could be produced by con- 
necting the output of an amplifier back in series with its input. More 
recently it has come to be known as an open circuit stable, negative 


jx 





(a) 


— The two types of negative impedance: (a) Open Circuit Stable and 
(by § Fort Circuit Stable. 
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impedance, because it will not oscillate when its two terminals are 
open circuited. 


THE E3 REPEATER UNIT 


The E3 repeater is essentially a two-terminal network the impedance 
of which has a magnitude | Z | and a positive phase angle that can vary 
with increasing frequency from plus 90 degrees, or less, through plus 
180 degrees to at least plus 270 degrees. This type of negative impedance 
is shown in Fig. 6(b). It has been known as the shunt type because it 
could be produced by connecting the output of an amplifier back in 
shunt with its input. In more recent years it has become known as the ~ 
short circuit stable type because it will not oscillate when its two termi- 
nals are short circuited. 


THE NEGATIVE IMPEDANCE CONVERTER 


The amplifier circuits of both of these repeater units perform the 
same function: that of a negative impedance converter. The operation 
of such converters is illustrated in Fig. 7. 

Fig. 7(a) shows the converter as a four-terminal network having a 
ratio of transformation k and a shift of phase through a negative angle 
of approximately 180 degrees over the operating band of frequencies. 
If, as shown, an impedance Zy is connected to terminals 3 and 4 then the 
impedance seen at terminals 1 and 2 will be the impedance Zy multi- 
plied by the ratio k and shifted in phase through a negative angle of 180 
degrees. This impedance will (over the frequency range of zero to in- 
finity) fulfill the definition given for the impedance presented by the E2 
repeater. Hence, Fig. 7(a) can represent the operation of the E2 repeater. 
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7 — The Negative Impedance Converter: (a) E2, open circuit stable and 
(b} 3, short circuit stable. 
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Fig. 7(b) shows the same converter, but here the impedance Zy is 
connected to terminals 1 and 2. The impedance seen at terminals 3 and 
4 (at least over the frequency band of interest) will be Zy divided by k 
and shifted in phase through a positive angle of approximately 180 de- 
grees. This impedance will (if frequencies from zero to infinity are con- 
sidered) fulfill the definition given above for the E38 repeater impedance. 
Thus Fig. 7(b) can represent the operation of the E8 repeater. 

From Fig. 7, it is apparent that the same converter circuit could have 
been used for both the E2 and the E3 repeaters. For practical reasons 
it was not. However, the ratio k and the phase shift in both the con- 
verter of the E2 and that of the E3 were made approximately the same. 


OPERATION IN TRANSMISSION LINES 


Within limitations, the E2 repeater can be represented by a negative 
impedance, —Z, and the E8 repeater can be represented by a negative 
admittance, — Y. With a negative impedance and a negative admittance 
available, losses of transmission lines can be reduced in the manner 
illustrated in Fig. 8. The transmission line is represented by two net- 
works as shown in Fig. 8(a). One of these (Network A) is in the form 
of a T network the series arms of which are represented by impedances 
Z; and the shunt arm, by an admittance Y. This network has a propaga- 
tion constant a; + j6,. The attenuation a; represents the major portion 
of the line attenuation, and the phase shift 6; is that just sufficient to 
make Network A realizable physically. This representation is necessary 
because Network A has image impedances each equal to the character- 
istic impedance (Zo) of the line. If the characteristic impedance of the 
line were a pure resistance, then the phase shift through this network 
could be zero and 6, could be zero. But the characteristic impedances of 
actual lines are not pure resistance; thus the phase shift 8; must be in- 
cluded in Network A. The other network (Network B) is shown as a 
box. It has a propagation constant az + j@.. Here ® represents the 
remaining phase shift in the transmission line and a: is an attenuation 
just sufficient to make Network B physically realizable in view of the 
image impedances which are both equal to Zp, the characteristic im- 
pedance of the line. Fig. 8(b) shows the addition to this line of a repeater 
consisting of a T network made up of negative impedances —Z in the 
series arms and a negative admittance — Y in the shunt arm. The arm 
—Z of the repeater adjacent to the line cancels Z of the line. The two 
admittances —Y and Y cancel and the other series arms —Z and Z 
also cancel. The result, as shown in Fig. 8(c), is that only the attenua- 
tion and phase shift of Network B remain. 
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In practice the amount of attenuation (a1) which can be canceled by 
the repeater depends on the uniformity, the loss and the type of line. 
The permissible magnitude is computed by conventional methods which 
are beyond the scope of this paper. 

Fig. 8 has shown how the combination of a series and a shunt repeater 
can annul a large part of the attenuation of a telephone line. Much may 
be accomplished also by a series negative impedance alone as illustrated 
in Fig. 9. In Fig. 9(a) a transmission line is again represented by two 
networks A and B. However, Network A now is shown as an ZL con- 
figuration having a series arm Z and a shunt admittance Y together with 
an ideal transformer of ratio 1: where K exceeds unity. Network A 
of Fig. 9 is equivalent to Network A of Fig. 8, and Network B of Fig. 9 
is the sameas Network B of Fig. 8. In Fig. 9(b) the addition to this trans- 
mission line of a single —Z such as the E:2 repeater is shown. This nega- 
tive impedance —Z cancels the series impedance Z of Network A. The 
result is shown in Fig. 9(c). 


NETWORK A NETWORK B ss 








XH +jAr 
pee REPEATER AND TRANSMISSION LINE------~-— >| 





(Cc). 


Fig. 8 — Operation of the ‘‘T”’ in a transmission line: (a) transmission line, 
(b) repeater and transmission line, and (c) result of addition of repeater. 
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Thus Fig. 9 shows that when a single —Z is added in series with the 
conductors of a transmission line the attenuation is reduced, but the 
equivalent of a shunt conductance Y together with an impedance trans- 
formation is left. The impedance transformation 1:K could be cor- 
rected by means of a transformer if it were not for the shunt element Y. 
Thus an impedance irregularity is introduced, and this irregularity re- 
flects power back toward the source. When the reflected power becomes 
a significant part of the total power passing through the line, transmis- 
sion is unsatisfactory. Echo is excessive and therefore the use of the 
single series repeater is limited. 


THE BRIDGED T STRUCTURE 


The discussion of the T repeater, illustrated in Fig. 8, was based on the 
use of two series negative impedances and a shunt negative admittance. 
It is perfectly possible, and more economical, to obtain the same effect 
by using a single series impedance in a bridged I’ structure and this 
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Fig. 9 — Operation of —Z in a transmission line: (a) transmission line, (b) 
repeater and transmission line, and (c) result of adding repeater. 
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is, in fact, what is done with the E23 combination repeater. Fig. 10 shows 
how this is accomplished by using a center tapped line coil for the 2 
repeater with the E3 connected as a shunt element to the midpoints of 
the coil. If the coil is considered to be ideal this arrangement is equiva- 
lent to the configuration shown in Fig. 11 in which the coil provides the 
basic bridge structure: the E2 repeater is the series arm; and the E38, 
the shunt arm. Incidentally this arrangement of E2 and E33 repeaters 
is similar to G. Crisson’s twin 21-type repeater.* 

For a bridged-T structure, the image impedance equals the square 
root of the product of the series and shunt arms and the attenuation 
(in db) is as indicated on Fig. 11. 

If a network is to be inserted in an electrically long line without in- 
troducing an irregularity, its image impedance must match the char- 
acteristic impedance of the line. This would be the case for the bridged T 
network if Z4 were set equal to NZ and Zz were set equal to Zp divided 
by N. Then the square root of the product of Z, and Zs would be Zo. 

A network make up of negative impedances is designed to match a 
line in the same way and Fig. 12 is a representation of such a structure. 
Here, as a matter of convenience, the shunt arm is shown as an impedance 





Fig. 10 — The bridged T repeater E23. 
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with a +180° phase shift instead of an admittance as used previously. 
The letter N designates a numeric or proportionality constant. It will 
be observed that the product of the shunt and series impedances is a 
real or positive impedance and hence the image impedance is a positive 
impedance, Zo. The gain is determined entirely by the value of N. 
Thus if the characteristic impedance of a transmission line is known, 
together with the gain that the line can support without risk of oscil- 
lation, then N is known and the repeater network can be adjusted to 
give the required gain. 

The advantage of the bridged T as compared to a single series nega- 
tive impedance such as the E2 can be demonstrated by comparing the 
relative transmission gains obtainable from the two arrangements. Fig. 
13(b) shows the insertion loss of a single impedance Z4 connected in 
series with a transmission line having a characteristic impedance Zo . 
If Z,4 is a negative impedance such as that produced by the E2 repeater 
then the repeater gain becomes a function of N as shown in Fig. 18(c). 
If N equals 2 the gain is infinite and the system will oscillate. Thus NV 
must always be less than 2 where Z, is a negative impedance of the 
series or open circuit stable type. Practically, the impedance of the 
transmission line is not a constant Za but varies with termination, line 
construction and temperature. Thus N should be decreased until the 
negative impedance is always less than the sum of the two line im- 
pedances in series with it taking into account all possible variations in 
these impedances. 

The same limitations on N apply to the bridged T repeater of Fig. 12 
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Fig. 11 — Schematic of the bridged T network. 
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Fig. 12 — Schematic of the bridged T repeater. 
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as apply to the single series repeater. If N is 2 the gain is infinite and the 
circuit will sing. This similarity goes further. Assume that a single 
negative impedance Z4 equal to NZ_/180° is inserted in series in an 
electrically long line and N is adjusted for stability. If this series element 
is removed and the bridged T of Fig. 12 is inserted in the same place 
and adjusted by changing N until the system is stable it will be found 
that N will have the same value in the bridged T structure as it had 
for the single series negative impedance. 

Thus, if N is the same in the case of the bridged T as in the single 
series impedance, the gain advantage can be obtained by comparing 
formulas on Figs. 12 and 13 from which it can be seen that the gain 
advantage of the bridged T is equal in db to 20 logiy [1 + (N/2)]. If a 
single series repeater can be used in a line to give an insertion gain of 
6 db (N = 1) then a bridged T can be used to provide 20 logiy (1 + 0.5) 
or 3.5 db additional. Thus, in this case the series repeater gives 6 db 
gain as compared to 6 + 3.5 or 9.5 db for the bridged T. These gains are 
theoretical; in actual lines with simply constructed repeaters the com- 
parison may not be quite so favorable to the bridged T. 





THE NEGATIVE IMPEDANCE CONVERTER 


So far the discussion of the E2 and E8 repeaters has been in terms of 
a “black box’ which translates a positive impedance into a negative 
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Fig. 13 — Insertion gain of the E2 repeater. 
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impedance through a multiplying and phase shift operation. It will be 
interesting to examine these boxes to see what factors determine their 
characteristics. 


THE E2 CONVERTER 


The E2 negative impedance converter is the same as the E11. As dis- 
cussed elsewhere® it can be represented schematically as in Fig. 14(a) 
and also in terms of the equivalent circuit of Fig. 14(b) if the coils are 
assumed to be ideal. The converter performs much like a transformer. 
An impedance seen through it is not only transformed in magnitude by 
the ratio of |1 — y8|to|1-+ u| it is also modified by the phase shift 
of this factor which over the operating band of frequencies approximates 
180 degrees. The symbol y» stands for the voltage gain of the electron 
tube and 8 is the ratio of | to 1 + (1/joCR). If both C and R are large 
8 approaches unity in magnitude and the ratio of conversion approaches 
1 —yztol+ yu. If wis large compared to unity then this conversion ratio 
approaches —1. This ratio of —1 is approximately realized in the E2 
converter, and therefore the conversion ratio is not changed appreciably 
by small variations in p. 
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Fig. 14 — E2 Converter; (a) schematic and (b) equivalent circuit. 
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In addition to the transformation term there is also, as shown in Fig. 
14(b), a series term, 2, divided by 1 + yu. Here R, is the plate resistance 
of the tube, and u is the voltage gain as mentioned before. The factor 2 
results from the use of two tubes in push pull. If » is large compared to 
unity then this series term becomes approximately 2f,/u. It is entirely 
dependent upon the characteristics of the electron tube. As the char- 
acteristics change from tube to tube with manufacturing variation or 
in the same tube over a period of time or with variation in battery supply 
potential, the term 2R,/(1 + x) will change accordingly. Percentage- 
wise this change may be large. This is the largest source of variation in 
the E2 converter. It can be minimized by operating the converter be- 
tween impedances much larger in magnitude than 2R,/(1 + yw) so that 
variations in this term have relatively small effect. This has been done 
in the E2 repeater by stepping up the impedance of the transmission 
line by about 1:9 by means of the transformer shown in Fig. 14. 


THE E3 CONVERTER 


Theoretically, the same converter used for the E2 and shown in Fig. 
14 could have been used for the E38. Instead of connecting the line to 
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Fig. 15 — E38 Converter; (a) schematic and (b) equivalent circuit. 
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the converter through terminals 1 and 2, terminals 3 and 4 would be 
used. However, because the E3 must be designed for connection across 
a transmission line a coil or transformer input is not practical since the 
coil would shunt the line at low frequencies and introduce excessive loss 
to dial pulsing and 20 eps ringing. Without a coil to step up the im- 
pedance of the line, variations in 2R,/(1 + yu) with standard triodes are 
too large to be neglected. For this reason another converter circuit was 
designed for the E3 repeater. 

This circuit is shown in schematic form in Fig. 15(a). It consists of 
tio resistances, a and c, respectively, and an amplifier poled according 
to the plus and minus designation on Fig. 15(a). The output impedance 
of the amplifier has been designated as b. If the input impedance of 
this amplifier is high compared to other circuit impedances, Fig. 15(a) 
can also be represented by the equivalent circuit of Fig. 15(b). Here is a 
conversion factor similar to that in the E2 converter and also a series 
impedance. The factor u is the ratio of the open circuit output voltage 
to the input voltage of the amplifier. In the E3 converter this voltage 
ratio uw is quite high because the amplifier is a two-stage arrangement. 
In the design of the E3 both ; and #2 are approximately one half. Thus 
uG, and pP: are both large compared to unity so that the conversion ratio 
(1 — 28:)/C + ue) is approximately unity and relatively independent 
of variations in uy. Furthermore, because uf, is large compared to b, the 
series term in the converter circuit is relatively small and variations in 
this term have little effect on the operation of the converter. 


Circuit DESCRIPTION 


THE E2 TELEPHONE REPEATER 


The circuit function of the E2 telephone repeater can be divided into 
two parts: the electron tube (negative impedance) converter; and the 
adjustable two-terminal network associated with the converter. 

In order to reduce the effect of variations in the electron tubes to 
negligible proportions, and at the same time to operate the tubes with 
load impedances that will permit optimum energy transfer from tube to 
connected circuit, the impedance of the telephone line is stepped up 
by means of the input transformer. To insure adequate balance for use 
in the telephone lines, the low voltage side of the transformer is divided 
into two equal, balanced windings. Each winding is center-tapped and 
connected in series with a line conductor. The circuit of the E2 repeater 
is shown in Fig. 16. 

In practice, it is advantageous to limit the conversion bandwidth so 
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Fig. 16 — Circuit of E2 repeater. 
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that the line and network impedances do not have to be controlled over 
an extended frequency band. The conversion gain is limited by reducing 
8 to a small value at low frequencies by capacitors Cs and C, of Fig. 16 
in the plate-to-grid coupling network, and at high frequencies by the 
small capacitors C; and C2 across the grid resistors Rs and FR. respectively. 

The conversion ratio is affected by small losses inherent in the elec- 
tron tube and transformer. These are balanced out by a fixed resistor 
R9 connected in series with the gain adjusting network to increase the 
amount of positive feedback. 

The final negative series impedance presented to the line is equal to 
approximately —0.1Zy over the frequency band of 300 to 3,500 cycles. 
The impedance Zy is determined by the configuration of the gain-ad- 
justing network comprising several inductors, capacitors, and resistors. 
These components may be arranged in any form to obtain the desired 
negative impedance, which in turn, introduces the gain and frequency 
shaping characteristic desired for each type of line facility. 

The E2 repeater employs a Western Electric 407A twin-triode elec- 
tron tube of the 9-pin miniature type. The tube heater circuit can be 
operated from 24- or 48-volt office battery. The heater current is 100 
milliamperes for 20-volt operation, 50 milliamperes for 40 volt opera- 
tion and the plate current is 11 milliamperes. 
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Fig. 17 — Schematic circuit of E3 repeater. 
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Fig. 18 — Circuit of E3 repeater. 
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THE E3 TELEPHONE REPEATER 


The E38 repeater is a two-terminal high impedance device designed to be 
bridged across the line in contrast with the E2 repeater which is a low 
impedance series connected device. Schematically shown in Fig. 17, 
its circuit is considerably different from that of the E2. The bridged 
telephone line and an adjustable network form two arms of a bridge 
connected to output and input circuits of a two-stage amplifier. 

The bridge is so connected that a fairly large proportion of the am- 
plifier output current flows into the telephone line and a fairly large 
proportion of the original line voltage is developed between the grid 
and cathode of the first tube since these circuits are connected to non- 
conjugate points of the bridge. The output and input terminals of the 
amplifier are connected to conjugate points of the bridge in such a 
manner than when the bridge circuit is balanced no feedback is effective 
at the input. 

It has been shown, in the preceding Section, that the negative im- 
pedance generated by this form of converter is equal to the network 
impedance Zy divided by a conversion factor. To obtain a practical 
design of the E3 repeater for the faithful conversion of the network 
impedances, with a minimum of spurious components, it is necessary 
to balance out, as nearly as possible, all converted circuit elements 
associated with the output bridge and connections to the line. Accord- 
ingly, the two line capacitors are balanced out by a network capacitor. 
The battery supply resistors, the resistor and capacitor of the plate- 
battery filter, and the plate-load resistor of the first-amplifier stage are 
all balanced out in the network by placing suitable values of resistors 
and capacitors in the cathode circuit of the second-amplifier stage. All 
of these elements combine to form an equivalent two-terminal network. 

An ideal negative impedance device would convert any impedance 
in the network over a wide frequency band but it is advantageous to 
limit the negative impedance, in so far as practicable, to the frequency 
bandwidth required by the particular application. This is accomplished 
primarily in the network associated with the converter. The conversion 
bandwidth of the E3 repeater is restricted at the low frequencies by the 
design of the resistive-capacitive feedback coupling network, between 
the output bridge and the input grid, and at the high frequencies by 
the shunt capacitance connected across one resistive arm of the bridge 
circuit. The circuit is shown in Fig. 18. 

The final negative impedance shunted across the line is equal to 
—Zy/0.94, within +2.5 per cent, over the frequency range of 200 to 
5,000 cycles. The magnitude and phase of the negative impedances are 


1080 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1954 


controlled by the configuration of the gain-adjusting network, consist- 
ing of several inductors, capacitors, and resistors. These components 
may be arranged in a variety of ways to obtain the gain and frequency 
shaping characteristics desired for each type of line facility. 

The E3 repeater employs a Western Electric 407A twin-triode elec- 
tron tube of the 9-pin miniature type. The circuit is arranged so that the 
current for the heater of the tube can be obtained from 24- or 48-volt 
office battery. The heater power is 2 watts and the plate current is 
nominally 5 milliamperes. 


EQUIPMENT 


The objective of the equipment design of the E2 and E3 repeaters was 
to produce a repeater that would be simpler to manufacture, easier to 
engineer, install, and maintain, and make more efficient use of the 
mounting space, particularly on 23-inch bays, than the present E1 
repeater. Because of the large demand, savings in manufacturing costs 
were realized by using a compact aluminum die-cast shell to house the 
repeater components. To facilitate manufacture, parallel thermoplastic 
strips were used for the mounting of “pigtail” type of components.° 
Further savings in shop costs were obtained by coordinating the designs 
of the E2 and E38 repeaters for maximum interchangeability of parts. 

Engineering and installation effort has been reduced considerably by 
avoiding engineered options and by arranging the equipment so that 
the maximum portion of the assembly and wiring work is performed in 
the shop. Testing and maintenance routines are simplified by arranging 
the repeaters as plug-in units that can be removed from their bay posi- 
tions and plugged into a portable test set located in a more convenient 
working space. In this way, the network strapping and any repair work 
which may be required on a repeater need not be performed while the 
repeater is in place, or is in a congested area. Maintenance and service 
interruptions are reduced to a minimum length of time by replacing a 
defective plug-in unit with a spare for immediate restoration of service, 
and the faulty repeater repaired at a more convenient time. 


REPEATER UNITS 


Both E2 and E3 repeater units use the same rectangular die-cast 
chassis, shown in Fig. 19. The front section of each unit carries the elec- 
tron tube and test pin jacks which must be accessible for testing and 
routine maintenance. The gain-adjusting network strapping terminals 
are arranged, in three rows, along the left side of the repeater, and are 
accessible only after the unit is removed from its mounting shelf. 
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Fig. 20 — Assembly and wiring of E2 repeater. 
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Fig. 21 — Repeater mounting shelves. Upper for 19-inch bays, lower for 23-inch 
bays. 


All external connections, to either type of repeater, are made through 
a male connector rigidly mounted into the repeater chassis, as shown in 
Fig. 20. The matching female socket is suspended in the mounting shelf, 
on a floating assembly, to relieve the strain on contacts and wiring when 
the repeaters are plugged into a shelf. Both connector units consist of 
molded rectangular phenolic blocks equipped with 15 gold plated con- 
tacts. Proper alignment between the male and female connectors is 
maintained by using a positioning key and guide pin on the repeater 
chassis, shown in Fig. 20, and a track on the mounting shelf. After a 
repeater is seated into its shelf position it is made secure by means of a 
screwdriver operated quick-acting fastener. 
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MOUNTING SHELF 


An aluminum die-cast shelf, shown in Fig. 21, is used for mounting 
the repeaters on the relay rack bays. The shelf comes in two widths; 
one holding 8 repeaters is used for 19-inch bays and the second, con- 
taining 10 repeaters, is used on 23-inch bays. Grooves cut into the base 
of the shelf match a projection on the repeater to act as a track system 
for positioning the repeater into its connector socket. The shelf connector 
is mounted in a small die-cast block which, in turn, is fastened to the 
shelf by shoulder screws to provide a floating assembly. A tapered key, 
molded into the repeater casting, engages a slot in the connector block 
to secure the horizontal alignment, and a tapered pin, at the bottom of 
the repeater casting, raises the connector block for the vertical align- 
ment. 

Fifteen mounting shelves can be arranged on standard 11-foot 6-inch 
_ relay rack bays. The maximum complement of repeaters will be 120 

for 19-inch bays and 150 for 23-inch bays. 


POWER DISTRIBUTION PANELS 


Fabricated power distribution panels which occupy 134-inches of 
mounting space are used for supplying plate and filament power to the 
repeater shelves. The primary unit, shown in Fig. 22, is required for 
furnishing power to the first three repeater shelves and one test-set power 
outlet. Four sets of plate and filament alarm-type fuses are provided, 
one set for each shelf and one set for the test power outlet. Four filament- 
adjusting rheostats are furnished and the panel is equipped with an 
alarm relay and lamp. Pin jacks are available on the front of the panel 
to measure the filament voltage. 

A supplementary power distribution panel is available, equipped with 
plate and filament fuses, rheostats and pin jacks sufficient for six ad- 
ditional repeater mounting shelves. One primary and two supplementary 
power panels are required to fully equip an 11-foot, 6-inch bay with 15 
repeater shelves. Both panels are completely wired in the shop to sim- 
plify installation. Fig. 23 shows a 19-inch relay bay mounting shelf 
complete with E2 and E38 repeaters and the two types of power distri- 
bution panels. 


REPEATER TEST SET 


A new test set has been developed, which will simplify the installa- 
tion and maintenance of these devices. This test set has been designed 
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Fig. 22 — Power distribution panels. Upper is primary unit, lower is supple- 
mentary unit. 
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Fig, 23 — Shelf of 2 and E38 repeaters, with power distribution panels, for 
19-inch bays. 
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not only to test the new E2 and E3 repeaters but also the older type of 
E1 repeater, which has been in use for several years. All test connections 
which are required to measure individual repeaters or combinations of 
El o- E2 and E38 repeaters are made by the operation of a single rotary 
control function switch. Such a system will avoid the unintentional 
errors and time consuming operations which attend the setting-up of 
complicated patch cord connections. 

The first five positions of the function switch are arranged to make 
transmission measurements of an individual repeater or any combina- 
tion of E1 or E2 and E38 repeaters. Insertion gain and loss measurements 
are performed with the repeaters working between their normally con- 
nected line impedances. 

In making transmission and stability measurements, it is sometimes 
necessary to set up trial connections of the repeater networks during 
the initial test period or for unusual line conditions. If the actual re- 
peater networks were used, it would require a large number of soldered 
connections to be made for each condition of the tests. Equivalent jack- 
ended E2 and E38 networks have been included in the test set, so that 
with small patch cords, a rapid interconnection of the network compo- 
nents can be made. 

As gain and stability of the E-type repeaters are both directly de- 
pendent upon impedance variations of the lines, it is sometimes neces- 
sary to measure line impedance to determine the causes for low return 
loss and instability. The last four positions of the function switch are 
_arranged to make positive impedance measurements of lines, negative 
impedance measurements of the repeaters and return loss measurements 
between any two lines or impedances. A decade resistance standard has 
been built into the test set to furnish a direct indication of the magni- 
tude of the unknown impedance. The phase angle of the impedance is 
determined by a return loss measurement and the angular degrees are 
read from a chart supplied with the equipment. 

Although this test set was designed primarily for tests on E type 
repeaters, it is possible to connect other types of four-terminal networks 
into the test set and measure impedances, return loss against known 
impedances, and insertion gains or losses. 


THEORY OF TRANSMISSION MEASUREMENTS 


One method of measuring the insertion gain of a negative impedance 
repeater, without affecting the transmission or stability of the circuit, 
is to introduce the test voltage, through a low impedance source, in 
series with the input line and measure the resultant current picked off 
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through a similar low impedance connection in series with the output 
line. A comparison of the received currents, with and without a repeater, 
will give an indication of the insertion gain or loss. 

The principles involved in making this type of transmission measure- 
ment are shown in Fig, 24. The voltage source and the current measuring 
detector are assumed to have negligible impedances compared to the 
impedances of the lines Z; and Z:. 

In the reference condition of Fig. 24, the driving voltage 4» and the 
current measuring device are in series with the two line impedances 2, 
and Ze’. The resultant current will be determined by the source voltage 
FE) divided by the vector sum of impedances Z; and Z2 : 


Zi + Ze 


where Z; + Zs indicates a vector addition of the two impedances. 

In the “measure” condition of Fig. 24 the repeater, represented by the 
network JN, is inserted into the circuit between the voltage source and 
the measuring device, and, although the voltage Ho is assumed to re- 
main constant, the addition of network N into the circuit will change 
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Fig. 24 — Principles of insertion measurements. 
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Fig. 25 — Test set showing repeaters and test cords. 


the current to a new value J; . The ratio of the currents [,/Io is a direct 
measure of insertion gain or loss between Z, and Z, due to inserting the 
network. The db change in transmission caused by the insertion of N is 
; eer I 
Insertion transmission in db = 20 logio 7 : 
0 
When J, is less than J) , the addition of N has caused an insertion loss, 
and when J; is greater than I» the addition of N has caused an insertion 
gain. In the test set the current indicating device has a db scale so that 
the change in transmission can be read directly in db. 


TRANSMISSION MEASUREMENTS 


In setting up to make repeater tests, the E2 and E3 units are removed 
from their shelf positions and plugged into adapters in the test set, as 
shown in Fig. 25. A test plug, fashioned to simulate the male connector 
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Fig. 26 — Sequence of transmission measurements. 
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of the repeater, is inserted into the vacant E2 repeater position for 
picking up the connections to the incoming and outgoing lines. A jack 
ended cord, connected to the test plug is patched into the test set for 
applying the line connections to the repeaters under test. 

The first position of the rotary function switch arranges the test set 
for the reference condition. The incoming and outgoing E2 repeater 
lines provide the terminations, as shown in Fig. 26(a). The low impedance 
voltage source is obtained by means of an oscillator working through 
a step-down transformer having an impedance ratio of 600:2 ohms. The 
low impedance side consists of two equal well balanced windings, one 
winding being inserted into each side of the line, to form a balanced-to- 
ground voltage source. The received current is measured by means of a 
detector working through an identical transformer having one of the 
2-ohm windings connected into each side of the line for maintaining 
the balanced-to-ground circuit. 

Switch position 2 connects the El or E2 repeater into the test cir- 
cuit between the sending and receiving impedances as shown in Fig. 
26(b). The change in detector reading will indicate the insertion gain 
or loss of the connected repeater. 

The third and fourth switch positions change the arrangements of 
E1, E2 and E8 repeaters for the insertion measurements of the various 
repeater combinations, as shown in Figs. 26(c) and (d). Fig. 26(e) in- 
dicates the connections for testing the E3 repeater alone on switch posi- 
tion 5. 


IMPEDANCE MEASUREMENTS 


The second section of the type-E repeater test set has been arranged 
to measure the impedance of two-terminal networks and the return 
loss between any two impedances or between an unknown impedance 
and a specified network. The impedance measuring circuit consists of a 
hybrid coil arranged in the form of a balanced bridge for measuring 
return loss. The driving voltage and detecting device are connected 
across conjugate arms of the bridge and the unknown impedance and a 
resistance standard are applied across the opposite conjugate arms of 
the bridge circuit. 

Simple transmission measurements made across the hybrid coil, 
between the oscillator and detector sides of the bridge, are the only 
determinations required to find the return loss between any two im- 
pedances or the magnitude and phase angle of an unknown impedance. 
Such a device has several advantages over other types of impedance 
bridges in that no critical balances are required, no calibrations are 
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necessary, the critical elements are stable passive networks and changes 
in the sensitivity of the measuring device, which includes variations in 
oscillator output and sensitivity of the detector circuit, will not affect 
the accuracy of measurement. 

An impedance reference is established, as shown in Fig. 27(a), when 
the transmission through the hybrid network is measured with con- 
nection to the resistance standard. A second transmission measurement 
is made through the hybrid coil network with the standard impedance 
connected instead of the unknown. After replacing the unknown with 
the resistance standard, the transmission through the hybrid coil is 
adjusted by varying the resistance standard to obtain the same value 
of loss, as measured with the unknown impedance. When the two trans- 
missions are equal, the resistance in ohms, as read on the resistance 
standard, will be the same as the magnitude of the unknown impedance. 
With the magnitude of the resistance standard and unknown impedance 
the same, the phase angle of the unknown is readily determined by com- 
paring the transmission through the hybrid coil for two polings of one 
of the hybrid windings. The first poling, Fig. 27(c), which is used as a 
reference, provides a measure of the transmission through the hybrid 
network when the current in the unknown branch and the current in 
the resistance branch are added vectorially in the detector winding. The 
reverse poling of one of the coil windings, Fig. 27(d), provides a measure 
of the transmission through the hybrid when the currents in the re- 
sistance and unknown branches are subtracted vectorially in the de- 
tector circuit. In addition to being a method of measuring the phase 
angle of an impedance, this comparison of two measurements with 
different poling is a return loss measurement of an unknown against a 
known impedance of equal magnitude. A curve of phase versus return 
loss may be used for convenient interpretation of the return loss measure- 
ment in terms of phase angle. After the phase angle has been ascertained 
it is necessary to determine the sense of the angle, that is, whether the 
unknown impedance contains positive or negative reactance. A reference 
condition is established, Fig. 27(e), by shunting a small value of capaci- 
tance first across the unknown impedance and then across the resistance 
standard and noting the change in transmission through the hybrid 
network. An increase in the transmission, in going from the reference 
to measure conditions, indicates a positive reactance and a decrease in 
transmission indicates a negative phase angle. The reference condition 
is established to minimize indicational errors at small phase angles of 
the unknown impedance. 

Two additional pieces of information are revealed in the return loss 
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measurement. When the return loss is positive, the unknown impedance 


contains positive resistance and the phase angle must lie between 0 and 
+90 degrees. For negative return losses the unknown impedance con- 
tains negative resistance and the phase angle must lie between +90 
degrees and 180 degrees. 


ACKNOWLEDGMENTS 


The design of the E3 repeater was based upon theoretical work by 
S. T. Meyers; the design of the test set was based upon theoretical work 
by H. Kahl and S. T. Meyers. In regard to the theory of negative im- 
pedances the contributions of F. B. Llewellyn are gratefully acknowl- 
edged. 


REFERENCES 


1. J. L. Merrill, Jr., A Negative Impedance Repeater, A. I. I. E. Trans.,69, Part 
2, pp. 1461-1466; 1950. 

2. J. J. Pilliod, Fundamental Plans for Toll Telephone Plant, A. I. E. E. Trans. 
71, pp. 248-256. 

3. O. B. Blackwell, Time Factor in Telephone Transmission, A. I. E. E. Trans., 
51, pp. 141-7. 

4. G. Crisson, Negative Impedance and The Twin 21-Type Repeater, B.S. T. J., 
31, pp. 485-5138, July, 1931. 

5. J. L. Merrill, Jr., Theory of The Negative Impedance Converter, B. 8. T. J., 
51, pp. 88-109, Jan., 1951. 

6. W. Ik. Kahl and L. Pedersen, Some Design Features of the N-1 Carrier Tele- 
phone System, B. 8. T. J., 51, pp. 418-446, April, 1951. 


Stress Systems in the Solderless Wrapped 


Connection and Their Permanence 


By W. P. MASON and O. L. ANDERSON 


(Manuscript received December 16, 1953) 


The solderless wrapped connection is initially held together by the hoop 
stress in the wire which enters the connection as a result of the tension put 
on the wire by the wrapping tool. Measurements made out to a time of 1.5 
years at room temperature show that the tension has decreased to 70 per cent 
of the one day value (8000 lbs per square inch) in this period. Two methods 
of extrapolation are discussed, both of which indicate that at least half of the 
initial one day value will remain at the end of forty years at room tempera- 
ture. 

Another set of stresses enters the connection as a function of time, namely 
the diffusion forces produced by diffusion of the tin plating into the brass 
terminal and copper wire. A number of experiments are discussed which 
show that the activation energy of diffusion is materially reduced by the 
‘shearing stresses tn the connection. Measurements at two temperatures, 
which allow extrapolation to room temperature, indicat. that at the end of 
two years the force required to strip the wire from the terminal has increased 
by & per cent over the initial value and that at the end of forty years the 
increase will be 20 per cent. Support for these conclusions is furnished by 
tests on actual connections that have been in the field for one year and ten 
months, which show an increase of 5 per cent in the stripping force even 
though the relaxed hoop stress is only 68 per cent of the initial value. The 
increase, which is due to diffusion forces, can be made higher by using zinc, 
cadmium or aluminum plating, and the fusion occurs in a shorter time. 


INTRODUCTION 


As discussed in a series of papers,’ the solderless wrapped connection 
is an efficient and inexpensive method of connecting a wire to a terminal. 
All the tests made so far indicate that it is mechanically sound and 
sufficiently free from the effects of corrosion to have a trouble free life 
of at least forty years. Photoelastic and stress studies show that the 


~ 1'The Solderless Wrapped Connection, B. 8. T. J., 32, May 1953. 
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connection is initially held together by the hoop stress in the wire which 
enters the joint as a result of the tension put on the wire by the wrap- 
ping tool. As time goes on, another system of stresses are gencrated, 
namely, the diffusion stresses caused by the diffusion of one part of the 
connection into the other which eventually eliminate the surface between 
the wire and the terminal and in effect join the two together. 

It is the principal purpose of this paper to describe a number of ex- 
periments which show how these stresses develop, how large they are 
and how they can be increased by substitution of a different type of 
plating between the terminal and the wire. A comparison of the strength 
of a joint formed by a tin plated copper wire and a bare copper wire 
both wound on nickel silver or brass terminals shows that the diffusion 
forces develop more quickly when the wire is tin plated than when it is 
bare. Measurements made at different temperatures show that the 
activation energy of diffusion is decreased in proportion to the hoop 
stress in the wire indicating that the shearing stress at the contact sur- 
faces aids diffusion. This activation energy is considerably less for tin 
than for copper. The diffusion joint has a strength per unit area equal 
to the limiting shearing stress of the tin plating. This is in the order of 
3,000 pounds per square inch for tin but is considerably higher for other 
types of plating such as zinc, aluminum or cadmium. Measurements of 
the stripping force of connections made with bare and tinned copper 
wire on zinc, aluminum or cadmium plated terminals show that the 
stripping force increases by a factor of two as a function of time, and 
the time required for the diffusion forces to operate is considerably less 
with these types of plating. . 

The combination of relaxation and diffusion stresses that are discussed 
later show that as the mechanical strength due to the hoop stress de- 
creases, the strength due to diffusion increases and at the end of forty 
years the standard tin plated wire on nickel silver. or brass terminals 
will be at least 20 per cent stronger than it is initially. The extensive 
corrosion tests described in Footnote 1, taken together with the mechan- 
ical strength tests described here, show that the standard connection 
should not fail in the forty year period under consideration. 


RELAXATION OF HOOP STRESS AS A FUNCTION OF TIME 


The rate of relaxation of the hoop stress in the wrapping wire is an 
important quantity for the stability of the connections. This has been 
studied by wrapping 24 gauge wire with a constant tension of three 
pounds around a spring steel terminal 0.0124 inches thick and 0.062 
inches wide. As shown by Figure 11 of Mallina’s paper,” this causes the 
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terminal to twist through an angle of about 25° when 100 turns are 
wrapped around the terminal. This twist is the result of a torque which 
is caused by the fact that the wire has a hoop stress and the wire does 
not come back on itself but advances by the thickness of the wire for 
each turn. As shown previously,” the total torque is equal to 





2ab 
T = (WLI. L 
orque = (WF) (22) (1) 
where (W.F.) is the wrapping force, 2a and 2b the thickness and width 
of the terminal and L the total length of the wrapped section. Hence, by 
calibrating the spring constant, the average hoop stress can be evaluated. 
In this manner, Mallina” has found that after transient creep (defined 
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Fig. 1 — Relaxation of stress in copper solderless wrapped connection at room 
temperature. 


in this paper as the loss in stress during one day) has occurred, the stress 
level for the standard 24 gauge connection is about 8,000 pounds per 
square inch. 

This same technique can be used to measure the loss in hoop stress as 
a function of time for all one has to do is to put a pointer on the spring 
and observe the rate the angle decreases. It will be noted that this 
system reproduces all of the stress systems and strain hardening that 
occur in the wrapped solderless connection, and hence a measurement 
of the relaxed angle is the most representative measurement of the 
relaxation of the solderless wrapped connection. Fig. 1 shows an average 
of the results of four such springs from the time of their intial formation 
out to a period of one year and one month. The rate of relaxation is 
quite rapid during the first day (transient creep) but at the end of one 


2R. F. Mallina, Solderless Wrapped Connections, Part I — Structure and 
Tools, B. 8. T. J., pp. 525-555, 32, May 1953. 
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year and one month the average hoop stress is 5,900 pounds per square 
inch and it is decreasing at the rate of 240 pounds per square inch be- 
tween the first and second year. In the previous paper,’ since the tran- 
sient creep occurring during the first day is somewhat variable, the angle 
of twist reached 24 hours after the connection had been made, was taken 
as the initial 100 per cent value and the relaxation has been plotted 
with respect to this. Fig. 2 shows this plot with data from Mason and 
Osmer and a set of eight relaxing springs measured by Mallina and 
McKettrick all plotted on the curve. At the end of 1.5 years, the ob- 
served stress is 70 per cent of 8,000 pounds per square inch or 5,600 
pounds per square inch on the average. 

The measurements out to 1.5 years by this method appear to be the 
longest measurements of stress relaxation in annealed copper wire in the 
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Fig. 2 — Relaxation of stress in tinned copper wire as a function of time and 
temperature. 


literature. Since the wrapped solderless connection is expected to last 
at least forty years, it is desirable to be able to extrapolate the stress 
value out to forty years time. Several methods are possible. According 
to theoretical treatment,’ stress relaxation proceeds as follows: 


o—-c= = log. (1 + t/t) (2) 


3 W. P. Mason and T. F. Osmer, Solderless Wrapped Connection, Part II — 
Necessary Conditions for Obtaining a Permanent Connection, B. 8. T. J., 32, 
pp. 557-590, May 1953. 

4See Progress in Metal Physics, Volume IV, Pergamon Press, Limited, 1953, 
Chapter V, Theory of Dislocations by A. H. Cottrell, pp. 233 and 251 to 260; 
Kuhlmann, D. and Z. Physik, 24, p. 48, 1952; and Proc. Phys. Soc., 64, p. 64, 1951; 
and A. H. Cottrell and V. J. Ayetekin, Inst. Metals, 77, p. 389, 1952. The last 
reference gives data for stress relaxation in zinc which shows that it obeys equa- 
tions (2), (3) and (4). 
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where oo is the initial stress, o the final stress, k Boltzmann’s constant, 
T the absolute temperature, fo a time determined from the equation 


= _kT BAty 
d09 = 3 log () (3) 
and A and @ are constants in the fundamental equation 

do _ __ A gH Ba) kT 

ie Ae (4) 


where H is the activation energy in the absence of any stress. 

The result of (2) is that for times large compared to fy , the stress dif- 
ference oo — o, plotted on logarithmic paper will be a straight line 
relation. Hence, in Fig. 2, for times in the order of forty years, the in- 
dicated stress level is 60 per cent of 8,000 pounds per square inch or 
4,800 pounds per square inch. 

Another method of extrapolation which essentially makes use of (4) 
is to measure the rates of relaxation at different temperatures and de- 
termine an activation energy for each stress level by comparing the 
logarithms of the times for the same stress levels as a function of tem- 
perature. The constant A in (4) is usually considered to have a tempera- 
ture factor T in it and is usually written as A’7’. With this relation, (4) 
can be written 

de sg OE Heer 


~ A Pe a Balke AT 


Integrating this equation between the limits oo and o, we have 


do = dt (5) 


oe a gegen = e Paolkry = a Pei oe presean ener (6) 





Hence, if the quantity 


gro euee 


is small compared to unity, we have 


(H — Bo) = k loge (t1/te) J, E — i (7) 


As shown by Fig. 3, the value of 8 is about 6 calories per pound per 
square inch and hence if op — a is at least 800 pounds/square inch (2.e., 
for all values of relaxation of 0.9 or less) this factor will be less than 1 
per cent for all temperatures used and hence we can use (7) to evaluate 
values of the activation energy H-Bc. 
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Fig. 3 — Activation energy for stress relaxation as a function of hoop stress. 


Measurements of the rate of stress relaxation for tinned copper wire 
at 200°C, 175°C, 150°C, and 100°C are shown by Fig. 2, and the activa- 
tion energies plotted against average hoop stress are shown by Tig. 3. 
Values are given using 150°C to 200°C as the temperature range and 
100° to 175° also. Both ranges give the same activation energies within 
the experimental errors and show that down to about 0.4 relaxation the 
activation energies satisfy an equation of the type 


H'= H — Bo (8) 


The curvature exhibited by the relaxation versus log ¢ shown for all 
temperatures indicates that the activation energy must increase faster 
than (8) for low values of relaxation and the dotted line shows a hypo- 
thetical curve ending up at the self diffusion activation energy for cop- 
per, 57 kilocalories per mole. 

To apply this method in general, one has to take account of any trans- 
formation such as recrystallization in the temperature range of measure- 
ment. For example, Fig. 4, shows similar curves for aluminum wire. 
Recrystallization in aluminum is known to occur at temperatures above 
150°C and this change is shown in the relaxation measurements by the 
lower values of relaxation that occur for long times. If one takes values 
of time above and below the recrystallization temperature, the activa- 
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tion energy will appear higher for this range than for a temperature 
range below the recrystallization temperature. The agreement of the 
activation energy for copper for the two temperature ranges shown by 
Tig. 3, shows that no transformation occurs from 25°C to 200°C and 
hence we can extrapolate the relaxation to room temperature taken as 
25°C, with the result shown by the solid line labelled 25°C. This agrees 
well with the measured values and indicates a stress at the end of forty 
years equal to 5,200 pounds per square in. 


DIFFUSION STRESSES IN SOLDERLESS WRAPPED CONNECTIONS 


In addition to the hoop stresses, another set of stresses develops as a 
function of time, namely, the diffusion stresses caused by the diffusion 
of one part of the connection into the other. The first experiment that 
showed the presence of these stresses was the stripping force tests of 
Fig. 23 of the previous paper referred to in Footnote 3. These measure- 
ments were carried out on connections which had been held at 175°C 
for lengths of time up to ten days and it was found that the stripping 
forces did not decrease with time. A more careful set with twenty con- 
nections for each point have recently been run with the results shown 
by Fig. 5, solid line labelled 175°C. From this curve, it is seen that the 
stripping force decreases to 88 per cent of its initial value of 15.5 pounds 
average for six turns of 24 gauge tinned wire and then increases to 120 
per cent of the initial value at the end of ten days. Similar increases are 
shown at 100°C over a longer period of time and recent tests of solder- 
less wrapped connections that have been in the field for one year and 
ten months show that the stripping force is about 5 per cent higher on 
the average than it was when the connection was formed. 
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Fig. 4 — Relaxation curves for aluminum wire in solderless wrapped connec- 
tions. 
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Fig. 5 — Ratio of measured stripping force to initial stripping force as a func- 
tion of time and temperature for tinned copper wire on nickel silver terminals. 


The dashed lines show the corresponding hoop stresses as a fraction 
of the initial hoop stress and hence it is evident that as the hoop stresses 
go down the stripping forces first decrease and then increase to higher 
values than were effective originally. We shall presently show that the 
difference between the measured stripping force and the proportionally 
relaxed hoop stress is due to a stress caused by diffusion of the tin in the 
plating into the wire and terminal of the solderless wrapped connection. 

One experiment which shows that the initial stripping force is caused 
by the hoop stress in the wire is the experiment shown by Fig. 6. Here a 
terminal is made which has a slightly tapered tin plated pin in the middle 
and is cut back for some distance beyond the pin. The terminal is wound 
with five turns of No. 14 gauge (0.065 inch diameter) tinned copper wire. 
The initial stripping force to pull off the winding was determined and 
it was found to take 72 pounds force on the average to strip the wire 
off the terminal. A similar set of measurements was made on the force 
required to pull the pin out of the terminal and this averaged about 50 
pounds or 70 per cent of the stripping force for the wire. Since the pin 
and terminal had the same coefficient of friction as the wire and ter- 
minal, and since the pin is required to support all the compressional 
stress in the terminal due to the hoop stress in the wire, it is evident 
that at least 70 per cent of the force required to strip the wire off the 
terminal is plain frictional force between the wire and the terminal. 
It is thought that the remainder of the force is due to the gouges cut in 
the terminal by the winding process. When the wire is stripped off the 
terminal, these cuts gouge out parts of the wire and hence require a 
higher force. This effect is equivalent to friction for a rough surface, 
which is higher than that for a smooth surface. 
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Next the terminal was held at 175°C for various times as shown by 
Tig. 6, where the stripping forces and the forces to pull the pins are 
plotted. It is seen that the force required to pull the pin decreases with 
time while the force required to strip the terminal increases with time. 
The pin force duplicates the stress relaxation curve initially but departs 
from this curve more and more as time progresses. This gradual depar- 
ture appears to be due to some diffusion in the tin which makes partial 
contact with the terminal. The amount of diffusion between pin and 
terminal is less than that between wire and terminal for two reasons: 
(a) the contact area between pin and terminal is not as intimate as the 
contact area between wire and terminal due to excessive plastic flow in 
the latter but not the former case, and (b) the contact area between 
pin and terminal is much greater than between wire and terminal. 
Nevertheless, the force required to pull the pin is substantially the same 
as the frictional force holding the wire on the terminal due to the hoop 
stress in the wire. Hence, we can conclude from this experiment that the 
initial stripping force is due to the frictional force resulting from the 
hoop stress plus shearing forces required to gouge the wire. 

Several experiments have been undertaken to measure the effect of 
diffusion separate from friction. The most successful of these was the 
arrangement shown by Fig. 7. Here a wire was pressed against a double- 
toothed sharp edged block, and a constant pressure was maintained for 
various times at various temperatures. An attempt was made to produce 
an indentation of the same magnitude as that developed in the wrapped 
solderless connection, although, of course, the area of contact increased 
slightly with time. 

It was found that for aluminum or copper on nickel silver at room 
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Fig. 6 — Experiment showing difference between stripping force and frictional 
force rdile to hoop stress — temperature 175°C. 
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temperature, no adhesion occurred even when the time of loading was 
very long. It was also found that if the load was removed, even mo- 
mentarily, before the sample was placed in the oven, no adhesion oc- 
curred. For a given temperature there is an induction period before the 
wire adheres to the block, and this induction period increased, in general, 
as the temperature decreased. This induction period was found to be a 
time of nucleation. This is shown by the fact that the period increases 
as the square of the contact dimensions. Since diffusion is a function of 
a//Dt where «x is the distance, D the diffusion constant and ¢ the time, 
this observation shows that nucleation starts at a given point and pro- 
ceeds for a certain fraction of the contacting surface before fusion 
strengths are observed. After the induction period, the fusion force 
—the force required to pull wire and block apart — increased at a 
rapid rate for the cases of copper, tinned copper, and zinc wire on a 
nickel silver block. 

In order to account for the effect of fusion due to the increase of con- 
tact area, the ratio of the fusion force to the contact area was determined, 
yielding a shearing stress. The area of contact could be easily measured 
with a microscope since a bright surface was produced by the shearing 
process. The shearing stress, Fig. 8, for the case of tinned copper wire 
on nickel silver is shown to approach a limit of about 3,000 psi which is 
approximately the limiting shearing stress for tin. Hence, it appears 
that tin diffuses into the copper wire and the nickel silver base. Further- 
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Fig. 7 — Diffusion forces for tinned and bare copper and zinc wire on nickel 
silver base. 
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Fig. 8 — Diffusion shearing strength for tinned copper wire on nickel silver at 
175°C as a function of time. 


more, it appears that the fused bond will just withstand a shearing stress 
equal to the limiting shearing stress of tin. Metallurgical evidence for 
this diffusion was presented in Fig. 19 of the paper referred to in Foot- 
note 4. In this figure, a tin plate layer between the copper and nickel 
silver was diffused in the two metals in a time of about 400 hours at 
180°C. The new evidence indicates that this layer of tin has the shear 
strength of bulk tin. 

The previously described data constitutes measurements of fusion 
forces separate from the complications of friction in the solderless 
wrapped connection. However, the same measurements of the effect 
of fusion are obtained if one subtracts from the stripping force the 
product of the relaxed force times the initial stripping force (about 15.5 
pounds) of the actual wrapped connections. As shown by Fig. 6, the 
fusion force so computed divided by the area of contact of a six turn 
connection (about 0.0045 sq in) yields the shearing stress, previously 
defined. In Fig. 9, the shearing stress is given for a tinned copper wire 
on a nickel silver terminal as a function of time for 175°C and 100°C. 
By way of comparison, the constant stress diffusion process of Fig. 8, 
is shown plotted by the dashed line on this plot. It is evident that there 
is a larger induction period in the case of the 14 gauge (0.065 inch wire) 
than in the case of the smaller 24 gauge (0.020 inch) wire and the ratio 
of the times is proportional to the square of the wire diameter ratio. 
This is an indication that we are dealing with a nucleation process. 
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The previous results show that the effect of temperature above room 
temperature on the solderless wrapped connection is to promote fusion 
which brings an additional set of forces into the picture with regard to 
the mechanical stability of the connection. In order to determine the 
effect of fusion at room temperature, it is necessary to evaluate the 
activation energy of diffusion. This is done by calculating the value H 
in the simple rate equation 


tT = to exp (H/kT) 


so that a given change of temperature corresponds to a given change of 
time on the fusion force curve. If we compare the 175°C curve with the 
100°C curve of Fig. 9, an activation energy curve as a function of residual 
hoop stress is obtained as shown by Fig. 10. The value of stress is deter- 
mined by the average time values used in the calculation of the activa- 
tion energy. It is shown that the activation energy for diffusion of tinned 
copper on nickel silver in the presence of stress is in the order of 20 to 24 
kilocalories per mole while that with no stress is in the order of 41 kilo- 
calories per mole. It appears that higher stresses lower the activation 
energy of diffusion just as they lower the activation energy for creep 
and stress relaxation. Very high stresses can reduce the activation energy 
to a value approaching zero which is the case of cold welding. The fusion 
occurring in the case of the solderless wrapped connection is different 
from cold welding in that it takes a finite time to develop any fusion 
forces. 

From the activation energy values of Fig. 10, it is possible to predict 
the rate at which the wrapped wire connection increases its strength at 
room temperature. Such a calculation applied to the case of tinned 
copper wire on a nickel silver terminal is given in Fig. 5 where it is shown 
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. Fig. 9 — Shearing strength of copper-tin-nickel silver connections as a function 
of time and temperature. 
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Fig. 10 — Activation energy for diffusion as a function of hoop stress. 


that the strength stops decreasing at the end of six months and even- 
tually increases to 20 per cent over the initial value at the end of forty 
years. Measurements on connections made one year and ten months 
ago, taken by V. F. Bohman, confirm this calculation. 

In all probability, the stress causing the diffusion is the shearing stress 
applied to the tin layer by the hoop stress of the wire acting on the fixed 
terminal. For a given initial winding stress, this shear stress will decrease 
in proportion to the relaxed hoop stress as confirmed by the data of 
Fig. 10. If we start with a lower winding stress, a smaller indentation 
is made in the wire and the area of contact is smaller about in proportion 
to the ratio of the hoop stresses. Hence, although the compressional 
force on the corners is smaller in proportion to the hoop stress, it is 
supported by a smaller area and hence the stress on the supporting 
area is approximately independent of the initial value of the hoop stress. 
This compressional stress produces shear stresses on the layer of tin 
since the material of the terminal and the wire have to slide with respect 
to each other in producing the indentation. Hence, we should expect 
that the diffusion forces would start at a time independent of the value 
of the initial hoop stress. 

This supposition is confirmed by the data of Fig. 11, obtained by 
winding tin plated copper wire on a nickel silver terminal with a hoop 
force half of the usual winding force of 1,300 grams for a 24 gauge wire. 
The initial strip off force was 60 per cent of that for a 1,300-gram winding 
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force confirming the data of Fig. 16 of Mallina’s paper, referred to in 
Footnote 2, which shows that the hoop stress decreases less rapidly 
than the applied tension. Fig. 11 shows the strip off force as a function 
of time and temperature. If we subtract the relaxed frictional force from 
the stripping force and divide by the area of 0.0027 square inches (60 
per cent of the area of the usual connection) the shear strength curves 
in pounds per square inch are shown by the lower solid curves. These 
start sooner than the values on Fig. 9 but require a somewhat longer 
time to reach their final value. The activation energy obtained from 
these curves is plotted on Fig. 10 (4800 psi) and shows that the activa- 
tion energy for the lowest stress is about the same as that for the higher 
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Fig. 11 — Stripping force and shear strength of copper-tin-nickel silver connec- 
tions as a function of time and temperature. Initial stress is 4,800 pounds per 
square inch. 


stress curve, showing that the shear stress is independent of the wrapping 
tension. It takes longer to complete the process and as a result, the strip- 
ping force at the end of forty years is equal to the initial value rather 
than 20 per cent higher as in the case for Fig. 5. Corrosion tests for these 
lower stripping force values’ show, however, that there is no indication 
of a loss of electrical contact. 

The method for analyzing fusion forces by subtracting relaxation 
data from stripping force data is verified by measurements on other 
metal connections. Fig. 12 shows the ratio of stripping force to intial 
stripping force for aluminum on nickel silver terminals as a function of 
time and temperature. Subtracting the relaxed force (the original strip- 
ping force multiplied by the relaxed stress ratio of Fig. 4) from the strip- 
ping force and dividing this difference by the contact area, the shearing 


5R. H. Van Horn, Solderless Wrapped Connections, Part III — Evaluation 
and Performance Tests, B.S. T. J., 32, May 1953. See Table II. 
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_ Fig. 12 — Stripping force curves for aluminum on nickel silver as a function of 
time for four temperatures. 


strength is obtained and is shown in Fig. 18. The shearing strength thus 
obtained approaches 6000 psi for long times at high temperatures. This 
value agrees approximately with the accepted value for bulk aluminum. 
Calculating the activation energy of diffusion as previously outlined, 
one obtains the curve for aluminum in Fig. 10. 

The advantage of a diffusion layer such as tin or aluminum as an 
element in the solderless wrapped connection is shown by the data of 
Fig. 14 which shows the rate of increase of the stripping force of bare 
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Fig. 13 — Shearing strength of aluminum on nickel silver connections as a 
function of time. 
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copper wire on nickel silver terminals at 150°C and 100°C. The increase 
in strength does not occur as fast as that for tin plated copper, and hence 
the standard connection with a tin plated copper wire is better than 
one formed from bare copper wire. 


PROPERTIES OF SOLDERLESS WRAPPED CONNECTIONS FOR OTHER TYPES 
OF PLATING 


A study of the factors governing the diffusion strength of the solder- 
less wrapped connection suggests methods for increasing the rate of 
diffusion and the — of the diffusing layer. These methods result 
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Fig. 14 — Stripping force and shear strength of bare copper on nickel silver 
connections as a function of time and temperature. 


in greater mechanical strength and a faster fusion of the wire and ter- 
minal. Experiments show that there are at least four metals which diffuse 
faster than tin, but due to the economic problems and the brittleness of 
some of the alloys formed, none of them are being considered for solder- 
less wrapped connection. 

The reason the diffusion forces in the tin plated solderless wrapped 
connection do not cause an increase in strength of more than twenty 
per cent during the life of the connection is that the limiting shearing 
stress in tin is only 3,000 pounds per square inch. If now we substitute 
for tin plating a plating with a larger limiting shear stress, the strength 
of the connection should increase by a larger factor. To be of use, how- 
ever, it is necessary that the diffusion forces shall develop rapidly. 

The data of Fig. 10 show that if we can produce a higher shearing 
stress on a layer of aluminum, the activation energy can be lowered and 
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Fig. 15 — Stripping force and shear strength of copper-zinc-brass connections 
as a function of time and temperature. 


be made to approach the cold welding condition. This requires that 
aluminum be placed on a stronger material such as brass, nickel silver 
or copper in order that the area of indentation for a given hoop stress 
will decrease and the shear stress will increase. Hence, the connection 
should form in a very short time. Furthermore, since the limiting shear 
stress for aluminum is near 6,000 pounds per square inch, one should 
except that the strength of the connection will nearly double. Since 
aluminum is not easily electroplated this suggestion probably is not 
practical. 

Of all the other metals examined, the next most promising are silver, 
cadmium and zinc. The activation energy for diffusion of zinc into copper 
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Fig. 16 — Activation energy for a copper-zinc-brass connection as a function 
of hoop Stress. 
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at no stress is given® as 38 kilocalories per mole which is 4 kilocalories 
higher than that for aluminum and 3 kilocalories under that for tin. 
The activation energy of silver and cadmium into copper also have 
low values. Hence, if the effect of stress parallels that for tinned copper, 
zinc and cadmium should diffuse into copper and nickel silver faster than 
tin. Furthermore, the limiting shearing 5000 pounds per square inch. 
Zinc, silver and cadmium arc readily plated on the terminals or the wire. 

Zinc plating has been tested experimentally by constructing a solder- 
less wrapped connection of bare copper wire on zinc plated nickel silver 
or brass terminals. The stripping force for a bare copper wire wrapped on 
a terminal plated with 0.001 inch thickness of zinc has been measured at 
175°C and at 100°C as a function of time with the results shown by Fig. 
15. The strength increases to 60 per cent over that found initially in a 
time of less than two hours at 175°C. At 100°C, the time required is a little 
over a day. If we subtract the relaxed force from the stripping force 
and divide by the area of the connection, the shear strength of the con- 
nection is as shown by Fig. 15, lower curves, for the two temperatures. 
T’'rom these two curves, an activation energy versus hoop stress can be 
obtained with the results shown by Fig. 16. These values allow one to 
extend the time variation of the shear stress down to room temperature 
with the result shown by Fig. 15. Adding these values multiplied by the 
area of the connection to the relaxed hoop stress, the indicated stripping 
force at room temperature is shown by Fig. 15, room temperature curve. 
This force increases at such a fast rate that the strength can be observed 
to increase at room temperature and corresponding measurements are 
shown by the circles. Since corrosion cannot occur in a region of fusion, 
a criterion for the corrodability of a connection is the time required to 
complete half of the total fusion at room temperature. Cn this basis, 
the zinc plated connection has the lowest half fusion time of any of the 
materials tested. 

Although zine diffuses more readily than other materials examined, it 
tends to form more brittle alloys and hence its use has not been seriously 
considered for solderless wrapped connections. 


6R. M. Barrer, Diffusion In and Through Solids, Cambridge University 
Press, 1941, Table 67, p. 275. 


A New Multicontact Relay for Telephone 
Switching Systems 


By I. §. RAFUSE 
(Manuscript received April 8, 1954) 


The trend in new telephone switching systems toward faster operation, 
longer life and lower cost, indicates a need for faster and more capable control 
circuits. This paper describes a new high speed wire spring multicontact 
relay designed primarily for these applications. The basic unit contains 80 
make-contact pairs. Two variations of the new design provide relays of 
60-contact capacity. They are mechanically and electrically interchangeable 
with all crossbar system multicontact relays. . 


INTRODUCTION 


In a modern dial telephone central office, many thousands of mo- 
mentary intraoffice control connections are made daily between the 
various parts of the switching equipment. For example, in the No.5 
crossbar system, seven major types of connectors’ are used to associate 
markers’ with other common control circuits, and with the switching 
frames, for brief intervals, to assist in setting up the talking connection. 
Connectors are required to simultaneously close a large number of circuit 
paths, as many as 240 in the trunk link connector. The earlier flat spring 
type multicontact relays’ used for this purpose provide large blocks of 
contacts per relay and provide an economical means for common or’ 
multiple wiring. 

The trend in new improved switching systems toward longer life, 
faster operation, lower cost and reduced maintenance, indicates a need 
for faster and more reliable connector circuits. This paper describes 
a new multicontact relay, designed primarily for these applications. 
It is a wire spring relay incorporating the improved manufacturing 
processes and many of the design features of the new general purpose 
wire spring relay.* The basic unit, Fig. 1, for use in new equipments, is a 
high speed 30-contact pair relay, with wiring terminals arranged for 
horizontal multiple connections. 


1111 
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Fig. 1 — New 30-contact relay with contact cover detached. 


This development also includes two modifications of the new wire 
spring relay for replacement use and for additions in existing crossbar 
equipments. They are 60-contact pair relays, each consisting of two 
30-contact units attached to a common mounting bracket. They are 
completely interchangeable with existing crossbar system multicontact 
relays. When a new improved relay is developed, it is almost invariably 
necessary to continue in manufacture and carry in merchandise stock, 
small demand codes of the old relays for many years. In this case, how- 
ever, manufacture of the old multicontact relays will be discontinued 
‘and all future needs will be supplied by the new product. 


OBJECTIVES AND REQUIREMENTS 


At the start of the project, all requirements from the standpoint of 
operating performance, circuit design and equipment use, were prepared 
in detailed form. The principal design objectives are summarized as 
follows: 


Electrical and Mechanical 


1. Operate and release times as fast as economically possible. 
2. Forty year life, or 200 million operations, with no adjustment 
necessary during the first 100 million operations. 
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3. No contact chatter. 

4. No false actuation due to armature rebound. 
5. No magnetic or vibrational interference. . 
6. 120-ohm and 275-ohm coils, to work with equipment already in use. 


Equipment 


1. Lower manufacturing costs. 

2. Reduced mounting space. 

3. Terminal arrangement for multiple wiring same as at present, or 
equivalent from a wiring standpoint. 


Mazntenance 


1. Contact failures due to dirt or insulating films should be sub- 
stantially equal to and preferably less likely than in the present relay. 

2. No contact locking due to contact erosion. 

3. Contacts should be replaceable in the field. 

4, Coil winding should be replaceable in the field. 

5. Field adjustment should be reduced to a minimum. 


Replacement Relay 


The design objectives also included modification of the new relay, if 
possible, to replace multicontact relays in existing crossbar equipments. 


Design History 


During the early stages of this development, considerable effort was 
directed toward improving the present flat spring multicontact relay. 
Later, many experimental models were constructed, to investigate other 
flat spring, and wire spring designs, and several contact actuating meth- 
ods. The most favorable designs of flat and wire spring multicontact re- 
lays were compared, and their differences were resolved by an analysis 
of manufacturing tolerances and their effect on performance and cost. 
Preliminary estimates of initial cost were only slightly in favor of the 
wire spring design. However, the wire spring relay offered significant 
advantages in (1) higher speed, longer life, less chatter, (2) better manu- 
facturing control of tolerances, (8) less maintenance, and (4) possibilities 
of future cost reductions as further improvements are made in mech- 
anized methods of manufacture. 


DESCRIPTION OF THE NEW RELAY 


General 


Stationary single wire and moving twin wire spring subassemblies 
are arranged in alternate layers attached to the core and mounting 
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Tig. 2 — Top view of relay showing location of parts. 


plate, as shown in Jig. 2. The wire spring assemblies form two rows, 
each containing 15 make contacts. Each contact pair consists of moving 
twin contacts on separate twin wires, associated with a single stationary 
contact. The stationary springs are supported close to the contacts by 
arms extending from the bracket. A detailed view of all parts and sub- 
assemblies for a 30-make contact relay is shown in Fig. 3. Moving con- 
tacts are pretensioned by relatively large pre-deflections as shown in 
Fig. 4. The method used in flat spring multicontact relays to obtain 
contact force is illustrated for comparison. It is apparent that contact 
force obtained by the “buckle” method depends on operating stud 
length and therefore is subject to change due to wear. The new preten- 
sioned wire springs are supported and actuated by a single molded 
phenolic card by the “card release’? method as illustrated in Fig. 5. 
In the unoperated position, the card is held against the core by a re- 
storing spring, which also supplies the force to open all contacts. In 
the operated position, the armature supplies the force to move the card, 
releasing the twin wires and closing all contacts. This method of actua- 
tion has some important advantages: 

1. Contact force is essentially independent of gauging and wear. 

2. The effects of wear at points in the relay which affect gauging are 
compensating to some degree, and therefore tend to minimize changes 
In gauging. 

3. Dimensional variations controlling contact separation and armature 
travel are reduced, making possible shorter armature travel, faster 
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Fig. 3 — Parts of the new 30-contact relay. 
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operate and release times, less contact chatter and increased mechanical 
life. 

Normally, the armature is held against the card by the lightly pre- 
tensioned armature hinge spring. However, when the relay is released, 
the armature motion becomes quite complex. Overtravel at the front is 
limited by the core plate backstop, and motion at the back, or heel sec- 
tions, is limited by the heel stop studs. This freedom of movement is 
intentional, its purpose being to dissipate armature energy into the core 
plate and core rather than back into the card. 


Magnetic Circuit and Armature 


Analytical and experimental studies’ show that one per cent silicon 
iron, with its higher resistivity, relative freedom from aging, and lower 
eddy-current losses compared to ordinary magnetic iron, provides opti- 
mum speed in the new fast relay. The contact load, about twelve grams 
per twin contact pair, is about one-half the load required in the present 
relay, and this together with winding space and heating,° largely deter- 
mines the size of the magnet. The magnetic structure is shown in Fig. 6. 
The core is a one piece “‘E”’-shaped section. The armature is a flat mem- 
ber made of low carbon steel having specific magnetic characteristics. 
This material simplifies manufacture, resulting in a cost saving with 
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no appreciable penalty in performance. The armature has adequate 
sections to carry the flux, optimum poleface area and lowest possible 
mass. Two small rectangular holes in the armature locate the base of the 
card in the horizontal direction only. The card is located vertically by 
the restoring spring as illustrated in Fig. 1. Fast release is obtained by 
a nonmagnetic separator strip, welded to the face of the armature. 
This strip also provides a smooth supporting surface for the molded 
card. Negligible wear at this critical point contributes materially to 
long life and stable adjustment of this relay. - 

A cellulose acetate filled coil’ is assembled to the center leg of the core. 
A nonmagnetic core plate, illustrated in Fig. 7, is then forced over the 
three core legs to hold them in alignment. The center hole in the core 
plate also functions as an armature backstop and permits a certain 
amount of overtravel of the armature when the relay is released. 


Coil Assemblies 


For circuit reasons, 120-ohm and 275-ohm coil resistances used in the 
old relays are required in the new relays. Nominal power savings which 
ordinarily would result due to an improved magnet and reduced load, 
are therefore sacrificed in the new relays in favor of increased speed. 
More than half of the new relays are expected to be used in circuits 
requiring maximum speed of operation and will, therefore, have 120-ohm 
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coils. In normal operation, these relays are not energized continuously, 
but operate only for short intervals while a talking circuit is established. 
The design provides for replacement of a coil winding in the field, but 
this is not expected to be necessary except in rare cases. 


M olded Wire Spring Subassemblies 


Wire spring subassemblies for a 30-contact relay are shown in Fig. 8. 
The two single wire subassemblies differ only in their terminal arrange- 
ment. The twin wire subassemblies are identical. Therefore only three 
basic molded parts are required which supply all needs for all new produc- 
tion relays. The wire spring sections are molded in continuous ladders® 
as in the general purpose relay. Spring bending, contact welding and 
coining? and terminal forming for solderless wrapped connections, ”° 
are performed in automatic tools developed by Western Electric Com- 
pany engineers. A comparison of the wire spring parts used in two new 
30 make relays and the corresponding parts of a 60-make flat. spring 
relay, is shown in Fig. 9. This illustrates the reduction in parts and 
simplicity of wire springs compared with flat springs. Seven types of 
subassemblies are used in the twelve layers required for an equivalent 
flat, spring relay. 

Terminals of the single wire subassemblies are formed for multiple 
wiring, and therefore differ in length and configuration. An improved 
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form of open wire multiple strapping has been developed for use with 
this terminal arrangement. It consists of bare wires held together in 
ladder form by means of phenolic plastic blocks molded successively in 
a continuous process. Fig. 10 shows relays with ladder type horizontal 
strapping soldered to the single wire terminals. The usual cable with 
new solderless wrapped connections is used on the twin wire terminals. 
Fig. 10 also illustrates the accessibility of multiple connections provided 
by locating them off to one side in the clear area between cable groups. 


Contacts 


All contacts in the new wire spring relays are palladium having a 
volume of contact material suitable for forty years life. This is equivalent 
to about 200 million operations for relays in high usage circuits. Con- 
tacts are easily visible, readily cleaned and may be insulated for test 
purposes. 

Contacts may be replaced in the field if necessary using Bell System 
field welding equipment.” Suitable tools and electrodes have been de- 
veloped to permit use of this equipment on all wire spring relays. 


Assembly and Adjustment 


The relay pile-up is securely fastened by two high tensile screws and 
a spring compliance member. Laboratory tests show that a pile-up of 





CORE PLATE 


lig. 7 — Core legs are held in alignment by the core plate, which is forced on 
the ends after the coil is assembled. 
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Fig. 8 — Molded wire assemblies for a 30-contact relay. 


this design will remain tight under widely varying atmospheric con- 
ditions. 

In the design of the new relay considerable attention was given to 
manufacturing control of tolerances, with reduced assembly and ad- 
justment costs as an objective. The molding process provides dowels 
for aligning the four layers of contact springs, trunnion supports to 
locate the single wires, and control grooves in the single wire subas- 
semblies to align the twin wires. These features provide accurate registra- 
tion for mating contacts, for the location of wire subassemblies relative 
to other relay parts and for the pretensioned restoring spring which in 
turn locates the actuating card in relation to the operating contacts. 
The card determines contact separation of all twin moving contacts 
in relation to their respective single stationary contacts. Due to these 
controls, and because the pre-deflected twin wire springs require no 
adjustment of contact force, no factory adjustments are anticipated 
except on relays which fall outside acceptable limits for back tension 
or contact gauging (or follow) as assembled. If necessary, therefore, the 
restoring spring may be adjusted to control the armature back tension. 
Contact gauging may be controlled if required by independent mass 
adjustment of the single wire contact rows. Fig. 11 shows this operation 
being performed by bending the bracket arm using a tool designed for 
this purpose. The mass adjustment feature is expected to simplify field 
maintenance practice. 
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Fig. 9 — A comparison of molded wire assemblies for a 60-contact replacement 
relay with the corresponding flat spring relay parts. 
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Fig. 10 — New 380-contact relays showing multiple wiring by the new ladder 
type strapping method, and the new solderless wrapped cable connections. 





Fig. 11 — New 30-contact relay showing method of mass adjusting stationary 
contacts. 
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Replacement Type Relays 


Future replacements of flat spring multicontact relays in existing 
crossbar equipments for maintenance reasons require a slightly modified 
form of the new relay, capable of complete interchangeability. As these 
relays are commonly used in connector circuits, operate and release times 
must be comparable with older type relays in order to avoid circuit 
interference. This was accomplished in the new wire spring relay with 
only minor changes in design and in the manufacturing process. 

As illustrated by Fig. 12, the interchangeable or replacement relay 
consists essentially of two 30-contact units assembled on a common 
mounting bracket having the same vertical mounting centers as the 60- 
contact flat spring relay. Since horizontal mounting space required by 
the new relay has been reduced, the 60-contact wire spring unit may 
also be used to replace 30-, 40- and 50-contact flat spring multicontact 
relays. The model shown in Fig. 12 has terminals arranged for horizontal 
multiple wiring. In another variation of the replacement relay all termi- 
nals are arranged for cable wiring. 

Modifications necessary to slow down the new relays for replacement 
use are (1) longer armature travel, requiring a different card; (2) an 
armature of larger mass — although this is a different piece part it 
has the same contour as the armature for the fast relay and may be 
punched in the same tool setup; (8) a core of low carbon steel in place 
of the one per cent silicon iron used in the high speed relay; (4) a leak- 
age reluctance shunt element shown in Fig. 13 to by-pass a portion of 
the total magnetic flux; and (5) coil windings having resistances of 120 
and 275 ohms as required for circuit reasons, but having the number of 
turns calculated for slower speed consistent with reliable operating 
capability.” 


Relay Performance 


Measurements have been made on laboratory-built models carefully 
prepared and adjusted to simulate extreme ranges of manufacturing 
tolerances, and more recently, on representative samples of pre-produc- 
tion relays. Although some of the performance characteristics studied 
will be determined accurately only after long-term use in the field, it 
has been possible by designed experiments and comparative tests, to 
obtain a fair appraisal of relay capability. These tests and measure- 
ments indicate that design objectives stated earlier in this paper have 
been substantially achieved. 
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Fig. 12 — A 60-contact replacement relay with one contact cover detached. 
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Load and Pull Characteristics 


Load and pull curves are measured under essentially static conditions. 
Spring load and armature motion are both observed at the center line 
of the card. They are measured and simultaneously recorded in chart 
form in a modification of a tensile testing machine.” A typical load and 
pull chart is shown in Fig. 14. The abscissa shows armature motion, as 
the armature moves the card and the spring load, through a distance 
of 0.030 inch to the operated position, and back again. The ordinate 
shows spring load on the armature on both operate and release, and 
also the magnetic pull which is developed in the armature for various 
numbers of ampere turns in the winding. 

The load and pull chart provides a comprehensive picture of over-all 
relay performance. For example, starting from the released position, 
the force or back tension, holding the card against the core is about 140 
grams. Following the upper curve, the spring load increases slowly as 
the armature moves toward the core, until the first contacts make at a 
load of about 200 grams. The load increases rapidly as the remaining 
contacts are closed until the last contacts are closed at about 650 grams. 
Further travel of the armature to the operated position increases the 
spring load to a final value of about 700 grams. As the armature is 
allowed to return to the original position, the lower curve is traced. 
The area between the two curves is a measure of mechanical hysteresis, 
or friction, in the relay. This energy loss is a very small fraction of the 
spring load at all values of armature travel. 

The pull curves show ampere turns necessary to assure operation of 
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Fig. 13 — Core assembly for replacement type relays showing leakage reluc- 
tance shunt. 
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Fig. 14 — Typical load and pull characteristics of a 30-contact fast relay. 


the relay. The maximum ampere turns required are determined by the 
“critical load point.” This occurs at 0.010 inch armature travel and 
about 650 grams. Under static conditions, therefore, 160 ampere turns 
would be required for complete operation. Circuit uses for these relays 
do not include nonoperate, hold, or release requirements. This informa- 
tion could however be obtained from the pull curves in a similar manner. 


Operate and Release Speed 


The new high speed multicontact relay operates two to three times as 
fast as its predecessor, the flat spring type relay. Operate and release 
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Fig. 15 — Comparison of operate and release times of wire spring versus flat 
spring multicontact relays. 


times are shown in Fig. 15. The improved performance of the new relays 
is shown by nominal operate and release values, and also by greatly 
reduced spread, or difference between minimum and maximum values, 
as compared with corresponding data for flat spring relays. 

Minimum operate times of the replacement and existing flat spring 
relays are comparable. It will be noted however, that operate and re- 
lease time spreads are much less in the new relay. This generally should 
improve the operation of existing crossbar circuits as new relays are 
used for replacements or additions. 


Contact Performance 


As speed increases, relays and other switching mechanisms become 
more susceptible to false operation of contacts. It is obvious also that 
faster operation adds to life requirements and therefore extends the 
period or increases the number of operations during which trouble-free 
contact performance must be provided. For these reasons, extensive 
studies were made of chatter, unprotected erosion, and locked contacts 
as applied to the new multicontact relay. Additional longer range tests 
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are planned to study protected contact erosion and susceptibility to 
open contact failures. 

When contacts of the new high speed relay are operated, initial chatter 
does not normally exceed 0.1 millisecond. There is no shock chatter 
caused by spring vibration due to impact of the armature with the core 
or backstop. There is no chatter caused by hesitation of the armature 
in its travel at the point where it picks up the contact load. This is due 
largely to the low mass of 0.020 inch diameter twin wires, the low mass 
and short travel of the actuating mechanism, the type of card operation, 
and the rigid mounting of the relay structure. 

Electrical erosion of contact material is reduced in the new relay, 
because of the reduction in chatter. For this reason, the contact size 
provided is expected to be adequate for all normal use for the life of the 
relay, and contact maintenance should be greatly reduced. 

Locked contacts are substantially eliminated in the new relay by the 
single card release method of operation. Static and dynamic forces 
associated with the restoring spring and card system are powerful enough 
to break loose any random pair of locked contacts. 

Contact failures due to dirt or the formation of insulating films on 
the contacts are difficult to check in laboratory tests. Long-term ac- 
celerated tests are necessary, with a large test sample, under carefully 
controlled dust conditions. Many precautions were taken in the design 
to minimize failures from this source, as follows: (1) a dust cover, Fig. 
16, encloses the contacts, but does not enclose the coil; (2) the cover 
partially segregates the contacts in groups of three pairs of contacts, 
reducing air movement in the vicinity of the contacts; (8) palladium 
contact material is used on all relays; (4) twin contacts are coined — 
the rounded surface reduces the area in contact, effectively restricts the 
area which may trap lint or other foreign matter, and increases contact 
pressure; (5) card release actuation and wire springs with large pre- 
deflections insure that no appreciable loss of contact force will occur due 
to age or erosion; and (6) twin contacts are attached to completely 
independent wires. 

Rebound chatter is another form of false operation which occurs in 
the form of contact reclosures caused by rebound of the armature after 
striking the backstop when the relay is released. Fundamental studies 
were made of rebound behavior in relay structures” and various models 
were constructed and measured. As a result there is no rebound chatter 
in the new high-speed wire spring relay within the range of normal 
adjustment. A comprehensive survey was made to determine the prob- 
ability of reclosures due to rebound, in relays having limiting adjust- 
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Fig. 16 — Contact cover for the new relays showing compartments in which 
contacts are grouped. 


ments. The probability of reclosures exceeding one millisecond duration 
is estimated as one in 28,000 relays. This grade of performance is due 
primarily to: (1) an armature having low travel and the lowest possible 
mass; (2) an armature suspension designed to dissipate rebound energy 
into the core plate and core rather than into the actuating card; and (8) 
a stiff mounting bracket to reduce the natural amplitude of core vibra- 
tion due to armature impact on operate and on release. 


Life 


Less than 5 per cent of the new multicontact relays will be required to 
operate more than 100 million times in crossbar systems during an esti- 
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mated life term of forty years. Life tests show that no readjustment 
should be necessary during the first 100 million operations. The tests 
also indicate that not more than a small percentage of relays will re- 
quire readjustment prior to an estimated maximum life of 200 million 
operations. In extreme cases, still greater life may be obtained if re- 
quired, by replacing the molded card. This is an inexpensive part and 
replacement is easily accomplished. 


Stability 


When the new relays are exposed to extreme temperature and hu- 
midity cycles, the greatest change in contact separation is, in general, 
about 0.004 inch, and only a small percentage of relays are likely to be 
used in this manner. Tests indicate that changes of this magnitude leave 
adequate margin for 100 million operations before readjustment is 
necessary. 

For economy, most equipment is shop assembled and wired on a 
frame basis, and shipped complete, ready for installation as equipment 
units. It is important, therefore, that apparatus units should be capable 
of withstanding physical shock far in excess of normal usage. Design 
features in the new relay which provide an adequate margin of safety 
in this respect are: (1) a rigid mounting bracket; (2) the wire spring 
pile-up is attached securely to the bracket with two specially heat 
treated steel screws; (8) the cover is held in place by the bending moment 
of an embossed section of a spring clip with a force many times greater 
than the compressive force of a single spring; and (4) guard surfaces 
molded in the cover prevent twin wires from leaving their respective 
guide notches in the single wire combs. 

Excessive shocks during shipment have, at times, damaged flat spring 
relays by bending their brackets. The new relays have been subjected 
to shocks of similar magnitude without damage. 


Magnetic Interference 


Under certain marginal conditions, a relay may be affected by leak- 
age flux from adjacent relays entering its magnetic circuit, and changing 
its operate and release values. Tests show that interaction is negligible 
between the new relays and also between new and old type relays when 
they are used in adjacent positions. 


CONCLUSIONS 


The new 30-contact relays provide faster operate and release times, 
longer life, improved contact performance, reduced maintenance, and 
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greater adaptability in new circuit and equipment units than previous 
multicontact relays. The new relays also require less vertical and hori- 
zontal space in new equipments. As a result of these improvements, 
substantial savings are expected when these relays are used. The design 
includes many features which permit the use of mechanized manufac- 
turing processes, which, in turn provide better control of tolerances. 
Yor these reasons, lower initial costs are expected as manufacturing 
and assembly methods continue to improve. 

The new 60-contact relays are completely interchangeable with all 
codes of flat spring multicontact relays in existing crossbar equipments, 
and in addition, they provide superior performance, longer life and 
reduced maintenance. Therefore, manufacture of the flat spring multi- 
contact relays will be discontinued as soon as new relay production 
becomes adequate for all uses. 
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Topics in Guided Wave Propagation 
Through Gyromagnetic Media 


Part III — Perturbation Theory and Miscellaneous Results 


By H. SUHL and L. R. WALKER 


Some problems, complete discussion of which would be extremely drfficult, 
are treated approximately by means of perturbation theory. Among these 
are the partially filled cylindrical waveguide, and the problem of multiple 
internal reflections in a sample of finite length filling the cross section of a 
cylindrical guide. Propagation in a ferrite between parallel planes, mag- 
netized along the propagation direction is discussed by the methods described 
in Part I. The paper concludes with an addendum to Part I — a numer- 
ical study of field patterns of the TEy,-limit and TM ,-limit mode for 
various de magnetic fields. 


INTRODUCTION 


Parts I and II of this paper were devoted toa number of specific propa- 
gation problems, whose solutions, though frequently quite complicated, 
could be discussed with a reasonably modest investment of effort. Un- 
fortunately, not all of these problems pertain to situations met with in 
actual gyromagnetic devices. Actual devices frequently employ struc- 
tures whose performance could be predicted only as the result of lengthy 
computing programs. For example, the microwave gyrator using Faraday 
rotation usually employs a ferrite sample whose cross-section only partly 
fills that of the cylindrical waveguide. Although it is easy to formulate 
the corresponding equation for the propagation constant, the classifica- 
tion and survey, let alone the computation of solutions, would be very 
difficult to carry out. 

Thus, one must often be content with approximate results, and the 
bulk of the present paper is devoted to perturbation methods. These 
take as starting point a situation whose propagation problem is essen- 
tially solved. The small change in propagation constant due to a slight 
change in the original state of the system is then calculated. The small 
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change of state may be the weak magnetization of an originally un- 
magnetized specimen occupying a substantial part of the structure, or 
the introduction of a very small specimen with arbitrarily large mag- 
netization into the originally empty structure. Under the heading, 
“Small Magnetization — Arbitrary Sample-Size,” we shall discuss the 
propagation constant for a pencil of ferrite of any radius, coaxial with a 
cylindrical waveguide, the space between guide wall and pencil being 
filled with an isotropic medium whose dielectric constant equals that of 
the ferrite. This is discussed in preparation for the practically more 
important case of a ferrite pencil of any radius in an air-filled guide. Here 
the unperturbed state of the system, when the pencil is unmagnetized 
and therefore isotropic, is already rather complicated and require some 
preliminary calculations. Under the heading “Small Sample-Size — Ar- 
bitrary Magnetization,” we consider the case of a thin pencil of ferrite 
in an originally empty guide. 

Another topic, not easily treated except by perturbation methods, is 
that concerning end effects in samples of finite length. After a prelim- 
inary discussion of internal reflections in an extended slab of ferrite (a 
problem which can be treated rigorously), two cases are considered: a 
ferrite slug of arbitrary length, closely fitting a cylindrical guide, and a 
thin disc normal to the guide axis, of arbitrary size. In these cases in- 
terest centers around the effect of sample length on Faraday rotation, 
though for the ferrite slug a subsidiary effect, that of mode conversion, 
is also mentioned briefly. 

It should be emphasized that the perturbation methods employed here 
are not in themselves novel. They are standard to most linear eigenvalue 
problems of physics, and have been used in connection with electromag- 
netic problems by many authors.’ 

The remainder of the paper is devoted to a discussion of a ferrite-filled 
“cable” in plane parallel form, using the methods of Part I. The treat- 
ment is kept in terms of saturation magnetization and magnetizing field, 
and is based on Polder’s equations. The paper concludes with an adden- 
dum to Part I, which reports some calculations and graphs of field pat- 
terns in a cylindrical waveguide completely filled with ferrite. 


1. PERTURBATION METHODS 
1.1 General Method 


A number of authors have made applications of perturbation theory 
to the problems of propagation in gyromagnetic media and the exposi- 
tion which follows is included mainly for completeness. We shall develop 
the subject in the following fashion: it will be supposed that the unper- 
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turbed system is a wave guide containing a medium whose permeability 
and dielectric tensors are diagonal and isotropic, but may vary over the 
cross section of the guide, although not in the z-direction along the guide. 
For this system it will be assumed that a complete set of normal modes 
exists for which appropriate orthogonality relations are known. The 
perturbation of the system will then consist of changes in the permea- 
bility and dielectric tensors of the medium, including the addition of 
non-diagonal terms. If these changes are to be genuine perturbations, 
they must be of one of two kinds. Either, the variation in the properties 
of the medium is confined to a limited region, small in volume in some 
appropriate sense, in which case its magnitude may be large, or, we may 
have a small fractional change in the material properties extending over 
a considerable volume. The fields in the guide may be expanded in the 
normal modes and a system of equations is developed for the 2-de- 
pendence of the amplitudes of these modes. These equations are then 
solved approximately, making use of the smallness of the perturbing 
terms. The results may then be specialized to the various situations of 
interest. 

Let us suppose that the unperturbed permeability and dielectric con- 
stant are w(x, y) and e(x, y) respectively and that the system is now 
altered so that it possesses a permeability tensor 


u(t, y)  — —Jjx(a, y) 0 
jx(a, y) uaa, y) 0 
0 0 us(a, Y) 
and a, dielectric tensor 
e(t,y) —Jjn@,y) 0 
jn, y) ex(z, y) 0 
0 0 é3(a, ¥) 


Maxwell’s equations for the perturbed system may be written, using 
the notation of Parts I and IIJ,f in the form: 





i 
V*H, — ae — jweli; — onl ;* = 0, 
G1 Es i 
2 02 JOM2L t Kit~ = VU, (1) 


V-H;* — joesE, = 0, 
V-E# + jousH, = 0. 


We omit the vector signs from all transverse vectors, which are sufficiently 
labelled by the subscript “t.”’ 
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These may be rearranged as 








: oH. . 
V H, —— D2 —= quel; = jules = eli; + anki; => At, 
* ae . * 
V*h, — + jou, = = —jJw (pe eae mi), — wk * = B,, (2) 


V — joel, = joes — a1), = 
VE + jomH, = —jo(us — w)H. = B., 


where A;, B;, A, and B, are introduced as abbreviations for the terms 
on the right hand side of the equations. L, and H, may be eliminated by 
substituting in the first two equations the expressions 


ee V-H* — A, 
Juey 
and 
Ge SVE Be 
Jou 


The two equations so obtained are 





- HH * 
ve (~ at \+ % eet = joel, = —A; anes = Vee Bs 
jw M1 M1 
and 
.U* 
Lye (YH) - ee + jomlls = ye a ea (3) 
jo € jw €1 


We now suppose that /; and H,; can be expanded in the form 
Ey = 2) an(z)Ein(2, Y) 


and 


H; = dX bn(2)H tn(2, Y), 


where 
E — iB nz d H —IBnz 
tn€ an tn€ 


satisfy the unperturbed form of equations (3) and the boundary condi- 
tion that tangential # vanishes at the guide wall. These equations have 
solutions for certain values of 8, only, but, if Ez, , Him are solutions for 
Bn = c > 0, then E:, , — Hi, are solutions for 8, = —c. We shall assume 
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hereafter that H:, and H:, pertain to positive 8, values. For a given 
bn (2) 
Gn (2) 
of propagation reverses. Substituting these series for H, and H; in equa- 
tions (3), one finds 


unperturbed mode it follows that reverses sign when the direction 





n | Gz qw Ta 
and 
n Zz Jw €y 


The orthogonality relation between functions of different n has been 
given by Adler.” It is 


if En*-Him dS = 0, 


where n ~ m; the tilde denotes the complex conjugate and the integra- 
tion goes over the cross section of the guide. We consider the fields to be 
un-normalized and write 


| But Ha a5 = Ln 
Clearly 
[ Bat Bon dS = —An . (5) 


Multiplying equation (4b) by Hi» and equation (4a) by #,,- and inte- 
grating each expression over the cross section we have 


—A, E ae ia = — [| AcE dS + Ea | Ba V* as, 
dz Jo My 





An EB oF Babs | = — [BeBe dS — = f BaD A; ds. 
dz jw €1 


Now we use the identity 
[ @-v+Pas = —[ PG@eds) + | PV-G*as. 
surface boundary surface 


This yields 
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Biggel. Pegg ca: [= V -Ein* dS, 


i. Ein * 7 
My b’dry My 


= jo | Bf dS, 


since L;,-ds = 0. H,, is the H, field appropriate to the nt» mode. Again, 
we have 
/ Hus ° we A: a0 


Bie ic -ds) + [4 V -Hin* d8, 
€1 b’dry € 
= —jo | Abn ds. 
since A, = 0 on the boundary. Z,, is now the Z, field of the n‘” mode. 


Making use of these results we obtain the equations for the amplitudes 
in the form: 


. 1 ji; a 


(6) 
dan x 2 i 
T+ Babe = —| — | BeBe aS + { AB dS |. 
Restoring the full iin: for A,, B,, Az and B,, we have 
on - Brn = x i (€2 = edly: Bim dS =o q | nl B in ds 
| (7a) 
+ i (us = ui) HHen as | ’ 
an + jBnbn = —— a (uz — ws) H,- Hin dS — jf ett Hin dS 
(7b) 


+ | (€3 —_ «Be fhan as . 


Equations (7a) and (7b) are, so far, exact, but they involve, on the right 
hand side, the functions Z,, H,, EH, and H, which are still unknown. 
We are interested in those cases where the integral terms are small, 
either as a consequence of the terms (e, — e), 7 and so forth being small, 
or of their being finite only over a small region. In the first case the fields 
E,, H,, #, and H, may be replaced in the integrals by the values which 
they would have before the perturbation was made. In the second case 
this is not possible since a large change in the material constants of a 
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region alters the field substantially within that region. Here, then, we 
have a preliminary problem to solve, namely that of determining the ficld 
in the perturbed region in a zero order approximation. 

Perturbation problems may be divided into two classes by another 
distinction. The changes in material properties may be independent of 
the z2-coordinate, so that the new problem is to consider propagation in a 
uniform guide differing slightly from the original one. Typical of such 
problems is that of a waveguide containing a ferrite rod of infinite length 
parallel to the z-axis; the perturbation consisting here of the change in 
the properties of the rod when it is magnetized. Clearly, in such eases, 
solutionsfor which all field components vary as exp — jG@z are still possible 
and the perturbation equations (7) become equations to determine 8. 
On the other hand, there is a class of problems for which the perturbation 
is confined to a limited region in the z-direction, and we are interested, 
perhaps, in the reflection and transmission coefficients for a wave inci- 
dent upon the obstacle. Here, for example, we might think of the case 
of a disc of ferrite across the guide. If we remain in the range for which 
perturbation theory is valid the changes in the amplitude of reflected 
and transmitted waves will be small, but the changes in phase may not 
be, if the perturbed region is sufficiently long. In the latter case, it would 
be possible, if the perturbation were uniform in z over the region in which 
it exists, to find solutions going as exp j@z, as described above, and to use 
these to fit the boundary conditions at the ends. It is also possible if the 
perturbed region is long, with slowly varying properties, to obtain suitable 
approximate solutions by the WIKB method. Some of these cases will 
arise in the examples which we treat below. 


1.2. Perturbations Uniform in z 


We consider first the general case in which the perturbation is uniform 
in z. In the absence of the perturbation the mt» mode is to be present. For 
the fields EL; and H, , in the perturbed region we write @n(Z)Eimo(x, y) 
and Om(z)Hzmo(x, y), respectively, where a,,(z) 1s the amplitude function 
for the m*» mode. If the perturbed region is one in which, for no magneti- 
zation, the material properties differ only slightly from their unper- 
turbed values, we may justifiably identify Lino with Li, and Heno with 
Hm. Tf the material properties are appreciably changed even in the 
absence of a magnetic field, Himo and Hzmo have to be calculated by 
an independent method. For a,, we put Ane’, where B = Bn + 68 
and 66 is small. Similarly for H,; and E,, we write b,, Himo and bn Emo , 
with b, = Bne”’. With such assumptions, the m* set of equations (7) 
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gives an equation for 8, while any other set, with n # m, gives the 
excitation of the n*™ mode. Substituting in (7) we have 


I 


IBmAm ca JBBm a ( [Cee = €1) E imo E im an jnE tmo*+ Etim 


Am 
(8a) 
+ (us = 111) Hemold em| as) Ay = jLAm 
and 
Joubm — IBAm = — (/ (ue = 1) tmo* H tm — jeH tmo** Him 
. (8b) 


i (es re €1)EzmoL em is) Bi = IMB. 


Ignoring squares and products of small quantities, one then has 
68 = 16 (L + M). (9) 


The first example to be considered is the effect on the propagation in a 
circularly cylindrical waveguide, when a coaxial, magnetized pencil of 
ferrite of very small radius is introduced. The guide radius is ry and that 
of the pencil is 71 . Before the ferrite is introduced, w1 = yo and a = 6, 
where po and eo are the free space values. The unperturbed fields are 
those of the usual TE and TM modes in round guide. It is necessary to 
calculate first the zero order electric and magnetic fields within the 
magnetized pencil; it will be sufficient to work out the magnetic case and 
deduce the electric one by analogy. Since the cross-section of the pencil 
is very small and transverse propagation effects consequently negligible, 
the internal field may be calculated by solving a static problem. The 
transverse magnetic field before the pencil is inserted is H;,, and it is 
assumed that the pencil is so small that over a circle with a few times 
the radius of the latter, Him is essentially uniform. We must now solve 
Laplace’s equation for a cylindrical rod immersed in a magnetic field 
which is to be uniform at large distances. Within the rod, B, = pH, 
— jxH;*, and at its surface the usual boundary conditions prevail. 
Hereafter we write p for pe. 

The fields are derivable from potentials ®..4, Pin , which are of the 
form 


Pin = (H imo); 


(a-7) 
a? 


Pout = (H tm*?) = 
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where r = (a, y), a is a constant vector and the coordinate system has 


its origin at the centre of the rod. Continuity of tangential H at the sur- 
face of the rod requires 


and then 
2 
as r 
Pout = Hum? + a (Htmo ~ Him): Tt. 


The normal derivative at the surface of the rod is 


1 
= [#5 tae ty ye 


1 - - 
a im _ (A tmo = Him) -1). 
1 


or, externally, 


Matching the normal B’s at the surface then gives 
o[2 Hum -— Hino = LH tmo —_ KA tmo*, 
or 


> * 
H imo es! Quol (u a pie = Jem ] ; (10a) 


In a similar manner, one would find if the dielectric constant of the rod 
were €, 





2eolt tm 
lim = —— - 10b 
Hang = 20 (106) 
The longitudinal fields H., and Hz, are unchanged within the rod. 
Turning now to the expression (9) for 68, we have in the present-case, 





Sai, ae Qeole — €0) ; _ sinh 
oe 7 2Am is a | ae 8 | Etim | a (e €0) | Loam | 

tg eas ; HitBln | 

Ho (ene re ioe ie tm 4Lim |, 


where we have anticipated that A,» is real, which we verify below. 
Since the integrand is constant the integral may be replaced by ary’ 
times the value of the integrand at r = 0. We consider now a TE-mode 
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with variation ¢”"’, n ~ 0. We shall have, 
Lim = —jJopoV*¥m , 
Him = —JBmV Vm , 
where 


Vn = Jan (i ry ei”? J i(ham) = 0 


0 


and 





Unm 
Bae = © €ofto ee es 
To 
For the fields on the axis, one finds forn = +1, 
: . Batten in 
A, = H = a ee. ies 
n 4H jn or, € 
and 


ing 





If |n| # 1 there is no first order perturbation. We now have 





2 2. 2 
Ww 2 €E€— €& 2 2Unm Hob Unm - + NK — Ho 

68 = —-— an E ofp? tm, 4 HB ton BE a = 
2LNGe e+ & To" 10° w+ nk + bo 


But we have 


ro 
Pc = an [ Bim™® +H tm? dr, 
0 


= —2n- enim [ E (a)? + we) | a ds, 


° 2 
Jn (ae 


= —2repoBm = 1). 


and then, 


2 2 
1 rn Unm 


Bn re Fulton)? (tan? = 1) (11) 


2b + NK — po 2€1 — €0 
E ae eee “| 
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Tor TM modes we have 
j 2 
0 
Him = Joe V *Xm 


where 
. %& jny : 
Xn = Jn Jnm —)é ? es © = 0 
To 
and 


Y 2 Jnm 
Bm == & €0 — —. 





To" 
Proceeding as before, we find, 
1 nv 1 [ 4.2 ecb mK = Mo | 
i6B=s>-°s° ry aren pi aE ea m . (12 
28m To” TGae para et & ey) 


A problem which is of some interest, although not of immediate prac- 
tical significance, is that of a ferrite pencil of arbitrary radius and infinite 
length in a round wave guide, with the remainder of the waveguide filled 
with a non-magnetic dielectric, whose dielectric constant, «1 , is equal to 
that of the ferrite. The ferrite is supposed to be only weakly magnetized. 
For such a problem, we have, 


66 _ oe i [Cu a 0) Him * Him = jtH im® + Hin] ds. 


Hm is the field of an unperturbed TE or TM mode in the dielectric-filled 
guide. » — wo and « are supposed small, but 7, , the radius of the pencil 
need no longer be small. 

For TE modes, we have as before (again excluding the case n = 0), 


Lim = —jopoV*¥m , 
Hum = —jBaV Yn » 


yh - 
hs, Unm ea eo" 
To 


2 2 
2. 2 nm 2 Unm 
Bn = @ eo — —> = Br — oe 
To ro 


Vn 


I 


and 
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Hence, assuming the pencil coaxial with the guide. 


OKs a aie 
=" 5g = Oe [w= w)Br’ E a (vm ") 
@: 0 To rT 





0 


: 
+ a 5 ie (eam “| rdr 
r To 
+ QrKBm” + 2 iia | Jia (on ") Jas (om x | dr, 
To 0 To To 
Unm (71/179) 
nm 1/709 j ga 
= ani | u = wo) f  2a'e(e) + 1] a 


Unm (71/79) 
+ 2«n f LOI ae], 
0 





2am (a z= po) Saad!) $e) 


r 
) + xnJ?, (om ")] : 
z=(r1/ToQ)Unm 70 


Making use of the value of A,, found in the preceeding paragraphs, we 
have 


= Bm Bo ’ x 12 


K 72 Ty 
Z=uUnm (71/7 9) Mo To 


Lim = =m VX ) 
Him => goa V Yin ’ 


é r j 
= Jn (jum) et 
To 


2 
2 2 dnm _ 
Ba = Wey ——> = Bi —-—— - 
To To 





(13) 





For TM modes, 


In this case, 


74m 68 = Qa (wer)” (« He — po) [ (ts 2G or ") 
pin y 12 (Jam 2) | dr + 200 [Fu (Jom =) Je! (om) ar) 
To 0 To To To 
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The value of A, for this case is 
j 2 
12 .. 
Am = —27weBm-: as Soe CFasads 


The value of 68 now becomes 


— Be BL _ ( ; ; x 5 
68 ~~ Bindnmed 2 (Jum) (e 1) tI nla) Jn (a) + a Fie (x) 


2 
i ae F@)) oe ae pe (im a) 
2 2=jnm (71/79) Ho To 


We note that for a ferrite filled guide with r, = 7» , the nonreciprocal part 
of 68 vanishes which confirms a result found in Part I of this paper for 
weak magnetization. 

The very high dielectric constant of the ferrites (about 10) puts rather 
severe restrictions on the size of the pencils to which perturbation theory 
is applicable, even for weak magnetization. This limitation would be 
- substantially relaxed if we possessed exact solutions for rods of high 
dielectric constant inserted into round guide, which could be used as the 
basis for magnetic perturbation calculations. Unfortunately, the only 
extensive published calculations of this kind are for dielectric constants 
less than 3. However, at the suggestion of M. T. Weiss, a calculation of 
the propagation constant of the lowest’ mode varying as e~’’ in a wave 
guide containing a coaxial dielectric rod (e. = 10) has recently been made 
in the Mathematics Department, for varying rod diameter, but for a 
single value of guide radius equal to 0.4 times the free space wavelength. 
With the aid of this information, which was made available to us, the 
magnetic perturbation calculation has been carried out. 

As before, the radius of the guide is 7) and that of the rod is r,. The 
dielectric constant of the rod is a . We consider first the propagation in 
the unmagnetized case. Since we are considering only one mode, namely, 
that with an angular variation, e”, and of the lowest order radially, we 
need not identify the £’s and H’s by a label. We use a subscript “1” for 
fields in the dielectric and ‘‘0” for fields in the empty part of the guide. 
In general, we have 


aE, = —jlonV*H, + BVEL, 
aH, = —j[@VH, — wV*E,], 





(14) 





where 
2 2 2 
a=we-— Pp, 


WH. = —oH,, 
VE, = —oE,. 
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and ¢, uw refer, for the time being, to the dielectric constant and perme- 
ability of the region considered. It will be convenient to put 


€0 77 € x BK = 
A. 4 = — H } oS 46; — = BB, 
Ho €0 Ho 
B 7 a 


SSS B; =a 
WV €oko OV E00 : 


1 
and to measure lengths in terms of ==. We shall continue to use 
wV oe 








V and V* with the understanding that they refer to the scaled units. 
We now have 


ja E, = pV*H., + BVE., 
jxH, = BVH, — eV*E., 





Oe = 8, 
VH. = —vH., 
WE, = —WwE, 
At the surface of the rod E,, H,, E, and H, are continuous. We must 
have 
iB oH, 1 oH, 
3, |B Barn -( \|-a[# Bin) — (22) | 
Qo? LTo or /o a? Or /1 
and 


| (2 ‘) +P ney] ge 2/a (4 ), +2 ae]. 
Qo or /o a? or Ti 
where 


€& = 


; ao =]—- B and ay = €& — B’. 


pa 
an | Gy 
o le 


fi = 1, everywhere, if the unmagnetized ferrite has the permeability of 
free space. 
These equations may be rearranged in the form: 


B( a ga) BOD = {ra aL" men (Gen 
+ gal ates (Ge), ] eo, 
(a a) a.) = {-2 [ners (B),| 








(15) 
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Within the dielectric, since all fields are bounded at r = 0, both /, and 
H, ave proportional to J,(ar) and, evidently, 


= Ty () = ry (2) = ar 1 (ear) = Flory) 
H.(r1) or /i _ Er) Or, 1 J 1(a171) wee 


in the notation of Part I. (#,)) and (H,), similarly, will be those two 
linear combinations of J;(@r) and Y1(aor), which, respectively, vanish 
and have zero normal derivative at 7 = 79, in order to ensure the van- 
ishing of the tangential fields there. The functions 

















CED ; r O(E.)o 
(H.)o or ie (Jo or 


will be called H(ar) and G(ar) respectively. 
Eliminating 1,(7;) and H,(r;) from equations (15) we obtain the char- 
acteristic equation of the problem in the form 


B( ] = 1 2 (Mex) = ye F(a) = Het) 


Qo” a,” a,” A” a,” OX? 


The perturbation in the present problem is that due to a mild mag- 
netization of the rod and referring again to equation (9) we have (in 
unscaled units) 


68 = re i [Cu = uo) H,: H, == jx *+ i] ds, 
21Na rod 
with 


An = E,*- Al, ds. 


guide 


Thus, in the scaled system, 


(4-1) lap Get | re 
i. 0 Ko Ho 
5 s 


TO ~ 
if EY*-Hyr dr 
0 





6B = (16) 


The evaluation of the normalizing integral in the denominator is an ex- 
ceedingly tedious business and it seems advisable to avoid it. This may 
be done in the following manner. The characteristic equation has been 
solved for numerous values of r; and 8 may be considered to be a reliably 


dp 


known function of r,. In particular the slope or is known. But we may 
ry 


also deduce a , by a perturbation calculation in which we start with a 
1 
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rod of radius 7, and increase the latter to 7, + dr, by changing ¢e to « 
in the shell 7, < r < 7; + dr, . For such a perturbation, since £, and FE, 
are continuous at the boundary, they suffer no change when the boundary 


moves; //, however is discontinuous and (H,), = e (E.); where (/,), and 
€0 

(Z,); are the normal fields just outside and just inside the rod. The per- 

turbation formula (unscaled), thus gives 


68 = 5 (a — ab( | P+ [BoP + |B): '). ae 


or, scaled, with r and 7; representing scaled radii, 
1@- DCF. [+ |B + «|: oe 
2 [ E*- Hr dr 





B= — 


The formula (16) for the magnetic perturbation may now be written 


1 dp (4-1) mek 55 Ghetto |r a 
rdry | 0 Mo Ko 


5p OD ee NE fe ee rear 
oo a ACA AC ee 


Inside the red, we may write 














rF H.(n) . 
H,=-—= LE, = jcH, 
and then 
jor H, = jcB ans —ak., 
or 
_ OH, 
jeu A or” 


The integrals are readily carried out and are as follows: 
x Ty Cs 
[| Bebrar 
0 


= FO] at + do Fen) + 2S) 4 BT) - oeba], 








| (A,*-H,)r dr 
0 


iB? 22. 
= 2 G + cB") — 2czB (Faun oe er) 5 “ +) I. 
ay 9 
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The term in the denominator may be evaluated by using 


go, 





jor E, = 


hae, 


ju Ey 
We obtain 
| Bry) P + | Bom) P + & | (A) P 
= E(n) E — cF(ar)) + ale — BF (aar)] a 





ry a4 


The perturbation may now be written: 








oh 
ae 7 
(17) 
(4 = 1 AG + cB’) —_ 2cBe| + a la + cB” = 2A cép] 
* 0 0 
r?ayt + [8 — cF(%rn)P + ale — BF(@n)P 
where 
2 
A= EOD) seg) ee 





and c may be obtained from equation 15, and the definition of c 


G(Gor1) _ z, F (ar) 


bot Oo? a? = ar G(aor1) = & a0 F’ (@ir1) 


1 1 B(@ — 1) 
B (2: - za) 


Ni(aoro) J 1 (&or1) 3 Ji(Goro) Ni (aor) 
Ni (Goro) J 1(Gor1) = J (Goro) Ni (Gor) 


Fig. 1 shows the propagation constant as a function of the relative 
diameter, 71/70, of the dielectric rod "8 the unmagnetized case, with 














with 


G(aor1) = Bort 


= 10. Fig. 2 shows the derivative 77—7— as a function of r,/r). The 


d(r - 0) 
guide radius is 0.4 times the free space wavelength. From equation (17) 
58 may be written as 


68 = P(u — mo) + Qk. 
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3.2 
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Vig. 1 — Propagation constant of a circular guide containing an unmagnetized 
coaxial, dielectric rod (€/eo = 10). r, is the radius of the rod and ro the radius of 
the guide. ro = 0.49 , where Xo is the free space wavelength. 


The computed values of Q and P — Q are shown in Fig. 3 as a function 
of relative rod radius. P — Q is plotted since P and Q are very nearly 
equal.* 


1.3 Perturbations Non-Uniform in z 


So far we have been concerned only with structures indefinitely ex- 
tended in the direction of wave propagation. In practice the non-recipro- 
cal element is, of course, finite, but end effects can frequently be ne- 
glected, since the element is matched at its ends (by tapering of the finite 


* H. Seidel and Miss M. J. Brannon, at the suggestion of M. T. Weiss, have re- 
cently calculated the dielectric loss for the guides containing a dielectric rod de- 
scribed above. By combining such information with that obtained here it is possi- 
ble to discuss figures of merit (degrees of rotation loss in db) for various pencil 
radii. Such an analysis is being made by M. T. Weiss and will appear in an ar- 
ticle, by S. E. Miller, A. G. Fox and M. T. Weiss, in a forthcoming issue of the 
JOURNAL. 
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sample, for instance). The matching could be accomplished in such a 
way that the transition region, whose characteristics would be very diffi- | 
cult to compute, should contribute little to the overall non-reciprocal 
behavior. Therefore, in many cases, the theory for the indefinitely ex- 
tended sample is adequate. For some special purposes, however, it is 
desirable to mismatch the sample deliberately. For instance, Rowen’ 
has suggested that the change in Faraday rotation, due to internal re- 
flections in an unmatched specimen, can offset to some extent the fre- 
quency dependence of the rotation which is implied by the Polder rela- 
tions, broadening thereby the useful bandwidth of the device. 
Consider an infinite slab of ferrite, magnetized in a direction normal 
to its two parallel plane bounding surfaces. A circularly polarized wave, 
normally incident on the slab, will be partially transmitted, and, since 
for such a wave, the medium behaves as though it had an ordinary scalar 
permeability, the phase and amplitude of the transmitted portion are 
readily calculated. It is clear that, as the result of multiple internal re- ~ 
flections, the phase of the emerging wave will differ from the value ap- 
propriate to a single trip through the slab (such as would be obtained 
were the slab perfectly matched). Both amplitude and phase of the 
transmitted wave will depend on the electrical thickness of the slab and 
its dielectric constant, and on the effective permeability. The latter 
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Fig. 3 — Q and P — Q versus r2/ro . 


differs for right and left circular waves, so that a plane polarized signal 
(which is the sum of equal right and left circular components) will 
emerge, in general, elliptically polarized, with the major axis of the 
ellipse tilted from the polarization at incidence by an angle which differs 
from the single trip value as the result of internal reflections. It is clear 
that this change in rotation can be calculated in a very elementary way. 

When the sample is confined to a waveguide a similar effect occurs, 
but its calculation becomes extremely involved, at least for arbitrarily 
large magnetizations. The reason, which should be clear from Part I of 
this paper, is that the circularly polarized modes no longer have the same 
field configurations inside and outside the sample.* Therefore any inci- 
dent mode excites all of the normal modes of the ferrite; these, in turn, 
give rise to all the mode patterns of the air-filled portion of the guide. 
Even if all but one of these are cut off, the excitation modifies the phase 
and amplitude of the reflected and transmitted portions of the propagat- 
ing mode. 


* This is due to the fact that there is now no ordinary effective scalar perme- 
ability as for infinite geometry. 
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Thus all modes have to be included in the problem, which conse- 
quently takes the form of an infinite system of linear equations for the 
mode amplitudes. This can be solved only to some approximation whose 
general validity it would be hard to establish. The problem could also 
be stated as an integral equation involving a complicated Green’s func- 
tion, with no greater chance of complete solution. 

We are therefore forced to restrict the problem to the ranges of mag- 
netization, or of sample size, in which perturbation theory is applicable. 
However, we begin with a discussion of the infinite, plane, loss-free slab, 
a problem which can be solved completely, and which has some bearing 
on the perturbation problem for a slug of ferrite whose cross section 
completely fills the waveguide. — 

Let the plane boundaries of the slab be normal to the z-axis, which is 
also the direction of magnetization and the propagation direction of a 
circularly polarized wave incident on the boundary z = 0 (see Fig. 4). 
In terms of the parameters p and o of Section 2.1, Part I, the effective 
permeability for a circularly polarized wave is 





p 
= —j 1 
BEEK Le 10 ( £7 .), 


the upper sign referring to right circular polarization. The correspond- 
ing propagation constants in the slab are then 


Bs = oan 4/1 & 2 


lag 





—_ 
INCIDENT Cc 
FIELD FIELD 


Fig. 4 — Normally magnetized ferrite slab. 
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where ¢ is the dielectric constant of the ferrite. If dis the thickness of the 
slab, the electrical thickness is 


2a d 9 
x= ped = 4 aa y/1 a P, 


where % = w+/ued is the electrical thickness of the unmagnetized 
sample. Let us now confine the discussion to the right circular wave. If 
the incident electric field is taken to be e 7 where Bo is the free space 
propagation constant, w ~/uoeo , the incident magnetic field will be 

ae Bo ge 0? 

Oo , 

and if the reflected electric field is pe”, the reflected magnetic field will 
be 


since {> reverses sign. Inside the slab, the electric field consists of forward 
and backward travelling parts 71e iB and 7.e+?, and the corresponding 
magnetic fields are 
= Bx T eo +? 
M4. 
and 


Bs. ret? : 


U4 


Finally the transmitted electric and magnetic fields will be denoted by 


~j8g(e—d 
736 58 y(z—d) 


and 

— Bo, gated) 

* Oo 
respectively. In general p, as well as the 7’s, will be complex. To obtain 
73, we write down the equations of continuity of all tangential fields 


across the boundaries z = 0 and z = d. Since the fields are confined to a 
plane normal to the z-direction, these equations are: 


l+p=n+r, 


Po (1 — p) = Ft (ry — ») 
Od LL9 Gio 
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and 
—j8.d ~jBsd 
73 = Tye tt A pe + : 


Bo 


—_— T3 = Bs (rie P+? = roe? +4), 
Who Go 


These four equations in four unknowns are easily solved for 7; . Writing 


ps toe Pp 
ty = (+k 
Bs, Bn 4/t /oy, 
w/ bo € 

a 





and noting that 


= Ly. ’ Say, 
where 
a= /* 
€0 : 
one finds the solution to be 
1 
i i ee 
Cos O% + 5 (+ ah *) Sin O24. 

= | 73 yews * 
where 

tan ®, = : (2 + ts) tan v4. (18) 

Ba : 

and 


[73 [4 = | eos Bors. + i(2 + =) sin” box | - (19) 
nO AND Sil 
Similar results apply to a left circular wave; it is necessary only to reverse 
the signs of o, p in the expression for 7, . Equations (18) and (19) show 
that equal right and left circular incident waves emerge with different 
amplitudes and phases; hence an incident plane polarized wave emerges 
elliptically polarized with its major axis inclined to the polarization direc- 
tion at incidence. The inclination and the ratio of minor to major axis 
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are determined as follows: The right and left circular fields, upon emer- 
gence, may be written in terms of rectangular components (with the 
polarization of the total incident field along the x-axis) 


J 


I 


+ ‘Mm + jut 
EH, —_ jy T+€ 3 

= “Th jut 
Ee + jE, = r-e”; 


from which the resultant field in the x-direction is seen to be 


Ler = Bo Be => | v4. | cos (wt — &,) + | 7_| cos (wt — P_) 
and in the y direction (20) 
E,r = £,° + EF = |r | sin (wt — 6_) — | 74] sin (wt — 4). 


The amplitude at time t, (Her + E,7’)"”, is thus given by 
Ep + Eyr = | 74 ? 4 Te 
+ 2 74.|-| 7] cos ot — @. + 4,)] (21) 


The major axis of the ellipse is the maximum of (Z,7° + E,7’)"” with 
respect to wt. It equals | 74 | + | 7_| and is attained at 


wt = lg@_ + ,). 
Similarly the minor axis is the minimum and equals | 74] — | r_ |. 
The ratio of minor to major axis is therefore 
[ret =e 
[ee | fre 


The angle between the z-axis (the incident polarization direction) 
and the major axis is found by substituting wt = }4(@_ + #,) in (20). 
This gives 


Pp, — P_ 
Lip = (|74| + | 7) cos ~{—, 


. B — € 
Eyr = (|r| + | 7-) a 
which shows that the angle is 


P, = P_ 
aaa; aa ° 


|z| and ® are plotted versus x in Fig. 5. 7, , ®; and 7_, &_ at given 
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1.6 


NEAR X2 = & » THE ELUCTUATIONS IN 
¢ AND |%| SUBSIDE, SINCE THE SLAB 
IS THEN MATCHED. BEYOND x2=£; 


THE FLUCTUATIONS RISE AGAIN 
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Fig. 5 — Phase and sara of transmitted circularly polarized wave as a 





function of xz = 1+ 


A 


, With =~ £10: (a), top, ® for 0 = and 37; (b), 
1 " Co €0 


bottom, @ for 6) = 77; (c), | 73 | for 60 = and 3z; and (d), | 73 | for 6 = 77. 
| o |, |p| are found by choosing positive and negative signs for « and p 


in 
= p 
= 4/472. 


Note that x can be imaginary corresponding to cut-off in the range 
—-l<oao< -l1—p(p <0). 
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Fig. 5(c), top and (d), bottom — See Fig. 5. 


Near « = 1 (corresponding either to small p or to sufficiently large o), 
® differs from its value © for the isotropic case by a small amount 6. 
The rotation, to first order, is then just one half the difference between 
the 6 for positive and negative p, o. Writing 





- = 1 p 
eae See eat ee eae 
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and expanding equation (18), we obtain 
cosh W sec” 0 — ban fo sinh Y 


SG “Pai, Ge a ee 
- ote 1 + cosh? V tan? 0 : 


€ Vv 
a= —~=€e. 
€0 


[This result holds even near 0) = (n + 14)a where tan 0 = ©, as can 
be seen by expanding the reciprocal of equation (18).] The quantity 
166(6x4. — 6x_) is the rotation corresponding to a single trip through 
the sample. The actual rotation is }4(66, — 66_). Hence we may define 
a rotation gain as the ratio 





where W is defined by 


3 (6b, — 6B_) 


9 ag, Se 
a me a) 500(624. = 5x_) ; 
an Mo (22) 


cosh © sec” 0 — : sinh VY 





0 


1 + cosh? ¥ tan? 6 


In many cases, % >> 1, that is, the thickness of the sample is much 
greater than a reduced wavelength in the specimen. Then the second term 
in the numerator of g is always negligible compared with the first. g then 
simplifies to 


cosh 


a cos? 0) + cosh? © sin? 6° 


This expression is plotted in Fig. 6 as a function of 6 for various a = e”. 
: : 5s i 1 1 
For given ¥ it has minima equal to ana at % = (n + 4 7,and max- 


ima equal to cosh V at 0 = n(n = 0,1, 2, +--+). Whena >1 (a~3 
for many ferrites), cosh Y is replaced by 4e” = a, and then 


=e 
91 max = 9 


—_— 


91 min 


Q 1b 


It is to be noted that when 6 = nz, the condition for maximum g; , 
the unperturbed reflected amplitude is zero, and the ellipticity vanishes 
to first order in éz. 
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Fig. 6 — Rotation gain, g:(@0), versus electrical thickness, 6> , for small mag- 
netization. 


Perturbation theory enables us to solve approximately, for weak mag- 
netization, the problem just considered for the case in which a right 
circular cylinder of ferrite snugly fits a cylindrical waveguide. We will 
show that, with a suitable reinterpretation of the constants, equation 
(22) for the rotation gain of the extended slab continues to apply here. 

Before magnetization, a particular right circularly polarized mode, 
say the m", is present in the sample. The small magnetization distorts 
the pattern of this mode and changes the propagation constant slightly. 
The distorted pattern can be expanded in the series of normal modes of 
the unperturbed material. In these expansions only the coefficient of the 
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originally present m‘" mode will be large; all others will at most be of the 
order of the perturbation. Denoting perturbed quantities by the superfix 
+, we have for the distorted fields 


+ _ —IBmtz 

Ein =€ id wel Dmnlitn ry 
+ —Ij8mtz 

Hin el as al mnt tn ) 


where Pm and Gm are large compared with the other coefficients. By, 
and the p’s and q’s are determined from equations (7a) and (7b). a, in 


these equations is identified with ¢ 7" nny On With e * mvs Ymn - Since 


all perturbation integrals involving H; and H,, /, vanish in the present 
case, we obtain 


BnQmn a BEDmn => £| fo ==, po) HtnH tn dS = 7 i KH t, agin as | 


and 


In the first of these equations, Hj, , the perturbed magnetic field, con- 
sists Of QmmlZ tm , plusan admixture of other modes with coefficients them- 
selves of order u» — po, x. Therefore, to first order, it suffices to write 


Htn = Qmmll tm 
in the integrand, with the result: 


BnQmn Ts BiDmn = Oahies 
+ 


Pmn = 


ae Me a 


where 
Ih, = = / (22 SG as: 


Elimination of Pmn gives 
(Ba = Bm dma = Bal mad mas 
The case n = m determines Bi: 
Bm = Ba(l — Timi) (23) 
All other cases give, to first order, 


Bales 
Quam = B — Qmm 3 
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which, again to first order, equals 


Bad sin 
B2 — B?, 


One of the quantities Pinm, Umm , is still arbitrary. If it is required that 
the perturbed field is, to zero order in I%,,, normalized to the same value 
as the unperturbed field, it is readily found that one can take 


_— Bm 
Ba. 

We are now in a position to consider the problem of a right circular 
mode, say the r’™ incident upon the end plane, z = 0, of the ferrite cylin- 
der extending to z = d. One simplifying feature of this problem is that 
the unperturbed modes inside, and the modes outside, the sample have 
the same dependence on radius since the sample fills the whole guide- 
cross-scction. However the modes inside and outside may have different . 
numerical coefficients. Thus, if we distinguish quantities outside the 
sample by primes, the TE, mode can be represented outside and inside 
the sample by 





Dmm « (24) 


Pam = 1; Qmm 


/ 
i, = i, = V *wWen . 
“sr va 8h 5 


“sT 


, , 
/ Tr sgn. By 
FA ,,, = S phig = Vise ? 
0 


ih, = Br 8Bn Br oy 
Qo 


Here sgn @ means: sign of the propagation direction, +1 and —1 for 
forward and reverse respectively. The function y,, is given by 


r ; 
War = Ss Ur — ers 
To 


where ¢ is the azimuthal angle, ro the guide radius, and u, the r" zero of 
J,/(x). Tor TM modes we have similarly 


Ei.. = Ei., = V vers 


I 


H,,, aT (A*Per) sgn B; , 


H:,, = a (V*Ysr) sgn Br 
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where in the definition of the y, wu, is replaced by 7, , the r’® zero of J,(2). 
Since the perturbations considered here do not couple modes with dif- 
ferent azimuthal number, s will be considered fixed hereafter, and only 
the suffix r will be retained. 

The field at z = —O0 (just outside the ferrite) will consist, not only of 
the incident field and its reflection p,, but also of other modes excited 
by the perturbation. Taking the incident FH,’ to have unit amplitude, we 
have, from the continuity of tangential electric fields, 


(1 + pr E tr + DS palin = DS (tm + tra) Do PaL ie, 
nfr n £ 
(25) 
= 2d, (Tar + Tn2)Pntlh te , 


where the tn, tng are respectively the forward and backward traveling 
amplitudes of the n™ perturbed mode. In the same way, we have for the 
tangential magnetic field: 


(I — rH tr —S 2 Pall tn = > (Ti — Tn2)QneH tt . (26) 


The changes in sign of the coefficients of the backward waves have al- 
ready been explained in Section 1.1. Let us now suppose that the incident, 
mode is the TE, mode. Multiplying both sides of (25) by V*¥, and both 
sides of (26) by Vy, and integrating over the cross-section we have, 
from the orthogonality relations of these functions, 


1+ pr = > (Tn1 + Tra) Dae ’ 


6,(1 _ pr) = > Brn = Tn2)Qne » 


But in these series, all 7 except 7, will be small, as will be all p, g except 
Drr} Grr Hence all terms, except those for which n = 1, will be small of 
second order, and can be neglected. Further, 


Dr = 1; Grr = 7 e 
Thus we obtain finally, 
1+ pr = Tr + Tr2 4 
2 
| Le = Pr = (Tr pe ra) gar 


In view of equation (23) for 8,*, these equations can be written, to 
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first order in I7, , 


I+ pp = tn + Tr2 


+ (27) 
1—p= Brito (Tr — T;2) 
Br Mr+ 





where 


Similarly at the output plane, z2 = d, the transmitted amplitude 7,3 
satisfies 


I 


—j6 70 
Trié + Tre ’ 


ny ’ . (28) 
a a (tae a Tr20"'), 
rber+ 


where 0 = 8,' d. Comparison of equations (27) and (28) with the cor- 
responding equations for the infinite slab shows that they have the 
same form; however, equations (27) and (28) hold only for small mag- 
netization. Hence only those results for the slab that relate to small 
magnetization can be generalized to the present casc. One such result is 
that for the rotation gain. If we re-identify variables by 


Tr3 





Tr = 





dat —> da,* = VT _ er 
2 B, 2 Br 
00 — Or = By d, 
a—> dy = B,/Br, 
we obtain, for the rotation gain, 
cosh y, sec” 6, — a us sinh y, (29) 


ue 1 + cosh? y, tan? 6, : 


where 
? 
e'=4,= Be/ Bes 


The formulae and conditions for maximum and minimum rotation gain 
derived in the previous section, apply here also, in terms of the re-defined 
variables. 

The conversion of part of the incident mode to the s™ mode in the 
transmitted wave can be examined by multiplying the matching equa- 
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tions by V*#, or Vv. and integrating over the cross-section. It is then 
found that the transmitted amplitude of the n™ mode is 


(Z, + Zn)(cos On — cos 6,) + 7(1 + Z,Z,) (sin On, — sin 6,) 





n => 2 ayn 7aW...n | Fo | F\ > NSAP... nN 1 “*/7F79 ) 4X25. nN YD 
inh SVE" OT, 608 9n +5 (Ze 1) sin bn) (22, 6080, + 5(Z2 + 1) sin8,) 
where 

6, = Br d, Len = 
and 
Raken 
Ym Bo BS 


The perturbation theory just outlined assumed a guide closely fitted 
with a slug of ferrite of finite length, and slightly magnetized. The per- 
turbation consisted in the small changes in permeability. Here we treat 
another kind of perturbation in which the ferrite does not fill the guide 
completely, and in which*'the magnetization can be arbitrarily large, 
but the sample is a thin lamina, mounted normally to the guide axis 
which only slightly perturbes the field pattern. Its shape can then be 
considered arbitrary; its thickness will be assumed! uniform and very 
much smaller than a wavelength in the material. Psd 

Under these conditions, the equations for the perturbed amplitudes 
of the right-circularly polarized n™ radial mode are 


an, . : . , 
a oF IBnOn = 2 | f (u as po) HH in dS — J [ tHe dS 
Be [& oe 6) LE, Ben as|, (30) 


Fe + itate = 5-5] f @- Be Baas | 


in the usual notation. The terms on the right are assumed different from 
zero only in the range 2, z -++ 6z occupied by the sample, and the integrals 
are extended over the cross section of the sample only. /:, , Hin are the 
mode functions for the empty guide. The fields H; , ZH; , Hz are the actual 
fields inside the sample. For the first-order calculation H; and H; are 
identical with the incident field and EH, (by continuity of the normal 


component of D,) is *° times the incident E, . Now the incident fields may 
€ 
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be taken simply as 
Anol’ tn ; AnolT tn Y Anon 


(since, in vacuo, bro = Gno Sgn By and sgn B, = +1 for the incident field), 
and the amplitude dn) may be taken to be unity. Therefore equation (80) 
may be written 


dan _ «ff 
= -i2|f (u — ps) Hinlen AS — 5 f edeg* Bley dS 
i ae = *) eoenLan as | = jBnDn ? 
dbn = gee | (e —_— €0) Lf; E, ds —- j6 An « 
dz An sae - 


We solve these equations by integrating through the small thickness of 
the sample. Since we are interested only in results of first order in éz, 
a, and b, on the right hand side of these equations may be replaced by 
ano, Ono ; that is by 1, 1. Writing 


I, = oe if Gi = po) HinH in dS — 4 ae Se: ds 
L\n LY J 


+ | (1 - = elnLen as | , 
In= 2 | ( - @)Bubin dS, 
we have, upon integration from the incidence plane (1) to the exit plane 
(2), 
Anz — An = J(Ia — Bn)é2, 
On2 — ba = JT — Bn)éz. 


Gm, the amplitude at the entrance plane consists of dno = 1, the incident 
amplitude, plus a small reflected amplitude p. Thus a,1 = 1 + p. Simi- 
larly b, consists of dno sgn 8, = 1, and the small reflected amplitude 
psgn B, = —p (corresponding to a backward wave). Dr2 is equal to Gre 
(since the transmitted wave is a forward wave). Hence we have 


Ang — (1 + p) = ja — Bn) 82, 
ang. — (1 — p) = js —_~ Bn) 62, 


GUIDED WAVE PROPAGATION THROUGH GYROMAGNETIC MEDIA. 11 1167 


whose solutions for the transmitted and reflected amplitude are 


Qng = 1 -i(6, ~ ttt) 5, 


and 


p= genes La 62. 


2 


Thus, to first order, the transmitted wave only undergoes a phase change, 


of amount 
50, = (6. ca iat ts) éz. 


A similar result applies to a left circular wave; the only difference arises 
in the contribution 


ae / cHin® Hm AS 





to “f which will merely change sign. The remainder of J, , and also Iz 


are unchanged. Thus the Faraday rotation of a plane polarized wave, 
which is one half the difference of the phase change for right and left 
circular waves, is equal to 





. we 2 “ ~ cya) = disc KH en* Hen oe 

The integral on the right is real, and, when the disc is circular of radius 
71, 18 of exactly the type that has already been encountered in the case 
of the ferrite embedded in a material with the same dielectric constant 
(Section 1.2). Here, however, x need not be small, so long as 6z is suffi- 
ciently small. In using the results of Section 1.2 it must be remembered 
that the mode functions there related to a dielectric, while here they 
relate to air-filled guide. When this is taken into account, one obtains 
for the specific rotation of a TE-mode: 


Ye) = Bn NK 2 "1 
ie Clan? — 1)d2 (tan) | ied (an 2) I, oe 


a TM-mode: 
iY) GF kes E 2 (; a) 
— = eee — I um — ; 3lb 
4 BmJnm?d i (dan) Ko To ( ) 
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results which show that the rotation in a thin disc is independent of the 
dielectric constant of the disc. 


2. THE PARALLEL PLANE CABLE 


Measurements on ferrites are sometimes made by a coaxial cable tech- 
nique. The ferrite fills a section of the cable, and is magnetized axially. 
In this section we consider a simplified problem of the same kind: that of 
wave motion between parallel conducting planes bounding a ferrite slab 
magnetized in the propagation direction. The characteristic equation for 
the propagation constant has been established by a number of authors, 
particularly Van Trier, but as in the case of the cylindrical waveguide, 
no analysis of the equation in terms of the laboratory variables appears 
to have been made. We shall make such an analysis by means of the 
methods already used in Part I for the cylindrical waveguide. The dis- 
cussion will be limited to the dominant (TEM-limit) mode, and to the 
principal features of the incipient modes (shape resonances). 

The distance between the two conducting boundaries will be denoted 
by 2a. The system of axes is as in Jig. 7; the y-axis is perpendicular 
to the walls (located at y = -ta), the z-axis is the magnetization and 
propagation direction. The microwave fields are assumed not to vary 
along the z-axis. 

The solution of Maxwell’s equation proceeds just as in the case of the 
cylindrical waveguide (Part I, Section 3), with the simplifying feature 
that one of the coordinates (x) drops out of the problem. We shall use 
the reduced notation of Part I. 6, the propagation constant will be meas- 
ured in units of w+/e, the propagation constant of the unmagnetized, 
unbounded medium. 4 will denote aw+/ye, the half-spacing measured 


in terms of the reduced wavelength on of the unmagnetized, extended 


medium. 






7 | 
CONDUCTING 
BOUNDARIES : 

\ 2a | 
\ 


\ 
\ 


| 
SY 


_ Fig. 7 — Parallel plane cable filled with ferrite magnetized along propagation 
direction, 
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It is found that (see Part I, Section 3) the longitudinal electric and 
magnetic fields are then given by 








Aw — Ars Wi — Ye 
UP = ————— 5 HH: = 
eG ss ae Oe 
where W1,2(y) are solutions of 
0 
Vas + x12 ti. = 0, (32) 
oy’ 


where A;,2 are the roots of the quadratic (16), Part I, and where, finally, 
the x’s are defined in terms of the A by equations (17a) and (17b), Part I. 
The x-component of FE is given by [see equation (22a), Part J] 


2 
0 
(Ai — A,)QE, = | esos = pura + va (2 — a] a 
Vi Oy 


or the same expression with suffixes 1, 2 interchanged. (‘oie that the 


x component of Vy is now zero, that of V*y is A If the boundary con- 
ditions are to be satisfied at both y = -++a and y = —a, the solutions of 
equation (32) must be either even or odd functions of y. Since #, must 
vanish at y = -ka, we have for the symmetric (or even) case 








COS x12 COS x22 
Wi _ 1 x y. yo" = Ae Y , 
COS X10 COS X20 
and for the antisymmetric (odd) case 
ie sin x11 Sin x2? 
Ww = Ay a Yo = Ag —— Xey 
sin x1@ sin x20 


The characteristic equation for these two sets now follows from the 
fact that LZ, = 0 at y = ta. Expressing the A’s in terms of the x’s by 
means of equation (17b), Part I, and writing BAe,1 = 1,2, We obtain 


: a& tan xia ! Xd tan xX 
X1 1 5X2 2! 
Mx Nox2” 


for the symmetric case, and 





2 tan X10 2 tan Xoo 
| ere reenicces AoX2 = 


Ax = = 
X1a Xo 


for the antisymmetric case. ‘lo bring these equations into a form suit- 
able for graphical discussion, we express the x in terms of A, o by means 
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of the quadratic for \, and by means of the relations between the tensor 
components of the permeability (see beginning of Section 4.11, Part I). 
We then obtain 


L 1 oi a 1— A? 
WV Tle Ves ex 
1 1 — are . 1 — dr 
= t 
Wis. ™=Vi-a, & 
for the symmetric modes, and 
1— A? - 1 — A? 1 — rd £ 1 — d2? 
t = t 
M1 1 — ov aa | 1 — od dof 1 — ode cam | 1 — ahe (34) 


for the antisymmetric modes. When )1, dz pairs satisfying one of these 
equations, and the Polder relation 


























Ar + Ae — ove = o + Pp (35) 
have been found, 6’ is given by 
B — —) Ao * 


Appearances to the contrary, equations (83) and (34) describe the 
ordinary TI and TM modes when p = 0, regardless of o (see the dis- 
cussion opposite Fig. 3. of Part I). When o = 0, the Polder relation trans- 

2 


2 
forms bean into eae so that either of equations (83) and (84) 


1 = ode 1 ct oN ; 
can be satisfied only by 


1-AY _1—-d oe ( Ve 











Lai. Looe ee or 


and the two alternatives respectively give © 


te 


9° 
- 


Apo = BP=1l—xn 


yee 


Since both equations (33) and (84) are satisfied by either the n-, or the 
(n + 14)- dependence, a problem of classification arises: “Which mode, 
TE or TM, is described, as p ~ 0, by the solutions of (83) and (84)?” 
Leaving aside the question of the TEM mode we note that the following 
degeneracy exists in the isotropic case: 


or 
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A TE mode with antisymmetric H, and a TM mode with symmetric 
i, have the same propagation constant 


2 1 2 
g= 1-5 (n+) (n = 0,1,2 °°). 
ita 2 
A TE mode with symmetric H, and a TM mode with antisymmetric 
E, have the same propagation constant 
2 2 
p=1-7t (n = 1,2, °°°). 
Since application of a very small magnetization will not destroy the 
symmetry properties, it follows that the solution of equation (84), and 
the solution of equation (83), which in the limit p = 0 reduce to 


2 © iL ; 
correspond to a TE-limit and a TM-limit mode respectively. Likewise, 


the solutions of equation (33) and the solution of equation (84) which 
in the limit p = 0 reduce to 


2 
2 27 


correspond to a TE-limit and a TM-limit mode respectively. When 
p ~ 0, the new #’ are thus given by different equations, so that the 
magnetization is seen to have removed the degeneracy of the isotropic 
case. 

In the special case of the TEM limit mode (8” = 1 when p = 0) only 
one of equations (33) and (34) has real roots. For, when p = 0, the 
Polder relation gives 


1—A = eee Ae 
1 — o\y 1—- OX” (36) 
= 0, 





if B’ is to equal unity. But equation (36) satisfies equation (34) only 
[equation (83) would demand \; = Az which is impossible for real £]. 
Thus the TEM mode exists as the limit of an antisymmetric mode only. 
to 


In fact it is easily shown that for a value of @ too small (tes than 5) 
admit any except the TEM mode in the isotropic medium, the only solu- 


tions of equation (83) for general c, p describe incipient modes. 
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Thus for d < 5 the analysis may be confined to equation (34) which 


will give the course of the TEM limit mode. The graphical analysis 
follows closely that of the cylindrical waveguide (Part I, Section 4.1); 
the part played by the G-function there is now played by the function 


_ =a ey eae 
b= r4/2=™ tana ee 


A contour map is drawn of the function Z in the upper half of the \-o 
plane. Only the upper half plane need be considered since the L-equation 
(34) is unchanged by reflection in the origin. This, incidentally, means 
that the present situation is reciprocal, so that co, p hereafter can be 
considered positive. Pairs of values \i(>0) and A2(<0) satisfying the 
L-equation can then immediately be read off, but while they will give 
8° = do, they do not necessarily, for given o, p, satisfy the Polder re- 
lation, equation (35). To ensure this, the quadrant, \ < 0, ¢ > 0, 
and the L-contours therein are transformed on to \ > 0, ¢ > 0 by the 
Polder relation for fixed p: 








¢o+tp-M 
l1— on 


The surfaces LZ; = L(\y, ¢) and Le = L(T(A,), oc) will intersect in a 
number of curves, along whose base curves in the A-o plane both the 
Polder equation and the L equation are satisfied, and along which #” = 
—)id2 is known as a function of o, p. The zero and infinity curves of 
L(Q, c).are denoted by 0, J when \ > 0, and by 0’, I’ whend < 0. The 
transforms of 0’, [’ onto \ > O are denoted by (0’)7, (1’)r. The suffix 
n denotes the infinity (zero) curves corresponding to 


1—N lV 9x 1—W © nt 

i =X =(n+$) Be? (; = X= ntZ) n integral. 
The lines \ = 0, +1, —1, are all zero curves denoted by 04, Oz, Oz’, 
respectively. The line Av = 1 is a conditional infinity curve, called J, . 
It is an J curve when viewed from od > 1 for \ < 1, and when viewed 
from od < 1 ford > 1. Otherwise (for oA < 1,4 < LlandoA > 1, > 1) 
it is a limit curve of all possible curves LZ = const., where the constant 
takes on any value indefinitely many times. (See Part I, Section 4, 
where the curve oA = | is a conditional zero curve, 0, .) Fig. 8, drawn 
for d ~ 1, p ~ 0.5, shows the part of the first quadrant allowed by the 
Polder relation divided into regions of like and unlike sign of L(, c) 
and L(T(A1, «)) by the various 0, J and (0’)7, (’)r curves. The un- 


de = = T(A1). 
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shaded regions are regions of like sign, and all these carry solution 
curves (dotted lines), by the same reasoning as was employed in Part I, 
Section 4.11. Two branches of the TEM limit mode exist; in the area 
bounded by (04) 7, Oz , (0s’)7 and in the area Oz, [1 , (O2’)r. The branch 
in the first region begins at ¢ = 0 with \1,A2 given by 


MV 1 — AZ tan G1 — AZ = Av 1 — AZ tan GV/1 — 2.2, 
Ai + Ae = Pp 


and ends at the intersection of (04)r, Os with ¢ = 1 — p and p’ = 0 
(since Az = 0 on (04)7). The branch in the region 03 , I, , (O2’) r begins 
with 8. = ©,o = 1 and proceeds towards ¢ = «, 6” = 1. 

The region bounded by J; , (04)7, J. contains an infinity of solution 
curves, the incipient modes. The n“™ of these begins at o = 1,» = 1, 
the intersection of J, and J, (the line I, is also the transform of \z = — ©, 


Og Oa. In ale Op (O's) + 





eae ne 
eee 


Sos 
— re 


oe 


Sass 


~— 
oe. 


oR 


(Ca); (0's); 
=| ie) 1 2 
<—A2 Ay 


Fig. 8 — Zero and infinity contours of L(A, ,¢) and L(7'(\1), c) in the first 
quadrant of the d - o plane. Cross-hatched regions are excluded by the Polder 


relation; shaded areas are regions of unlike sign of the two functions. Dotted 
curves are solution curves. 
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along which Z, = ©); crosses from one allowed region to the next 
through the intersection of 0, with (0,’)7 and proceeds to its cut-off 
point (0, — (O4)r). 

Formulae giving #’ in the neighborhood of special points on the solu- 
tion curves, or for special values of the parameters, are obtained exactly 
as in Part I, by expansion of the L-equation. The results are tabulated 
below. They relate to antisymmetric modes only. 


Cut-off of the TEM-limit mode (« = 1 — p, p” = 0) 
Near cut-off 
dee 1 
dc tan @)\ (37) 
P (1 ~ 8G 


For p > 1, the TEM mode has no lower branch. 





Behavior near resonance (¢ = 1): 
2 





vis 
TEM mode: p” = i lk Ns 
—o 
7 i, (38) 

n+ ;) 
n incipient mode: 6° = ((r+ 3)" Si js 
24? l—o 

(n = I, 2, be +) 


Cut-off of the n" incipient mode (parametric representation): 
6 8 a rig 
M1 = €; c= e (2 — £) + na é 


T 


p=ao—-%X 8 =0 (39) 


2 2 2 2 2 2 2 = 
. = ma( me — 1) coth 0 /| (™ — ets ota — ‘|. 
(The reader may note the similarities of these formulas to those 


for the cylindrical waveguide). 


Spot-point for the n™ incipient mode: 


« 2 « 2 
le 7 aa) / 6 |+ nit (40) 


p= (t+ oe — 1). 


Q 
ll 
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Approximate formulae for all antisymmetric modes at small p, 











oo Ts 
TM modes: 6” = Brit i se ee (41) 
TE modes: 6” = Bn,2° (1 i = (42) 
where 
: vee aie 
Pra = 1——S m= 1,2, ++ 
and 
is 
Bug = 1 — “—-— n= 0,1, 2,--- 
TEM mode: 
g=1- ?.. (43) 
— 6 


TEM mode for large o: 
g=1+2+0/(4). (44) 
o o 


TEM mode for very small a: 


for sufficiently small a, the 6” for the TEM mode is independent 
of a. For general o, p, except ¢ = 1, 


_1l-(+ Pp) 


~ lop —e 


6 (45) 


Reference to the expression for u and « in terms of o, p will show that 
the wave progresses as though the medium had a scalar permeability 


Mere given by 
1 1 1 1 
Licht) 
por 2\ee oe Kp 


3. FIELD PATTERNS IN THE FILLED CIRCULAR GUIDE 





In Part I we discussed propagation in a circular guide filled with 
ferrite and longitudinally magnetized. Formulae or the field com- 
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ponents were given there, and it is of some interest to examine numeri- 
cally or graphically the distribution of Z and H. Since all field compo- 
nents vary (in complex notation) as e%”’t®?, with n = 1 in the present 
case, it is clear that the field patterns are stationary in a system rotating 
with angular velocity, —w. Writing ® for the angular variable, @ + at, 
in this system, one finds the following formulas: 








E, = —E_ sin 24, E, = (L, — E_) sin ®, 
fy, = H_cos28+ L,, Ky = (£4, + E_) cos , 
E, = Ey cos ®, (46) 
H, = H_cos2®+ H,, H, = (H, + H_) cos 9, 
H, = H_sin 24, H, = (H_ — H,) sin 9, 
A, = Hysin 9, 
where 
1 
od +1 | 
Eh — OS _ 
(r) = OxWh Gare J 2(xar" ) 
{8s | 
yD} = — 
mG ) QxaJi(xaro) o(x17) 
HAG) = My 
(7) ore 7 Baoan) o(xir) — (47) 
H.(r) = My - 
+(7) ds are Cree o(xi7") 
Ji(xar) | | 
rH = _ 
o(r) J (xan 0) 2? 
1 Ji(xar) -| | 
H = : 
o(r) as J (xr 0) 2 


The terms in square brackets are in each case the same as the cor- 
responding unbracketed terms, A: and x2 replacing \; and x1. The quan- 
tities £4, and E_ are the amplitudes of the left-handed and right-handed 
components of circular polarization into which the transverse E-field 
may be resolved at each point. H, and H_ have a similar significance 
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for the transverse H-field. This may be readily verified by examining 
the vectors (f, + jE,)e"*' and (H, + jH,e*, which represent the 
transverse field vectors in the laboratory system. The transverse fields 
are elliptically polarized at any point and the ratio of minor to major 
axes of the ellipses are 


[eI= (Hl) og Hel = [HEI 
[| #4] + | £_|I || H.)+ | H_|| 
The fields so far are normalized only by the choice of a simple form for 
the function, Ko(r); all components may, of course, be multiplied by the 


same arbitrary constant. There is some virtue to a normalization based 
upon power flow. The power flow is given in unreduced units by 


i: (EX H). ds. 
guide 


Using the scaled units of part I with 


r 


Tactual = 
WV Moe 


and H replaced by 0 FT to give it the same dimensions as F (e is 
€ 





the dielectric constant of the ferrite), the power flow becomes in the 
present variables 


p= 4/2 | (iH, + ELH) rar, 


Mo Woe 


Tv 
= -4/* ; IT 
Uo W' Loe 


We shall normalize the / and H fields by dividing the values given by 
equations (47) by J’”. This makes the power per (isotropic wavelength 
in the ferrite)’ a constant. In Fig. 9 we show the normalized fields for 
the TEn-limit and TMy-limit modes as a function of r for the case 
ro = 5.75, | p | = 0.6 and several values of c. The behavior of the modes 
as a function of o and p for this radius is shown in Fig. 14 of Part I. 
We also show the amplitudes for the isotropic cases, > = p = 0. It may 
be recalled from the discussion of cut-off points in Part I that, for the 
TM mode at this radius, when cut-off is reached at c = —0.4, the am- 
plitude of the H™ field is overwhelmingly greater than that of the others. 
Even when normalized to the same power flow, the field amplitudes 
for a given o are undetermined to a factor of +1. This factor has been 




















RADIUS | 


Sanus 


n-Hmit and TMy-limit 
= 0.6, except for Figs. 9(f) and (p) 


= 0. The fields are normalized to a constant 
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Fig. 9 — Amp 
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Fig. 9(c) — See Fig. 9. 
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Fig. 9(d) — See Fig. 9. 
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Fig. 9(e), top, and (f) 
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Fig. 9(g), top, and (h), bottom — See Fig. 9. 
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Fig. 9(i), top, and (j), bottom — See Fig. 9. 
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Fig. 9(k), top, and (1), bottom — See Fig. 9. 
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Fig. 9(m) — See Fig. 9. 
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Fig. 9(p) — See Fig. 9. 
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Fig. 9(q) — See Fig. 9. 
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Fig. 9(r) — See Fig. 9. 
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Fig. 9(s) — See Fig. 9. 
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Fig. 9(t), top, and (u), bottom — See Fig. 9. 
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chosen for each pattern according to the following considerations. For 
large | o | the amplitudes should tend to those of the isotropic case; for 
| «| < oo the amplitudes should develop with increasing | ¢ | away from 
the isotropic amplitudes. Thus, if the phase of the isotropic amplitudes 
is first fixed, that of all the others may be fixed by the conditions that 


TRANSVERSE H-FIELD TRANSVERSE E-FIELD 


TEy—-LIMIT MODE 
oT = 06 





2 3 4 
RADIUS l 






TEy,—LIMIT MODE 
o = -0.6 
= -0.6 





4 3 2 3 4 5 5.75 
fs RADIUS RADIUS lo 


Fig. 10 — Transverse field patterns in a rotating system for four particular 
cases of Fig. 9. For any pair of E and H patterns the length of the arrows is propor- 
tional to field strength, but the patterns for different o-values have not been 
given a common normalization. (a), top, TEy-limit mode, a = 0.6, p = 0.6; (b), 
bottom, TE,;-limit mode, ¢ = —0.6, p = —0.6; (c), TMy-limit mode, « = 0.6, 
p = 0.6; and (d), o = —0.3, p = —0.6. 
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they shall fall into an orderly sequence. In Fig. (10) the actual transverse 
j- and H- field patterns are shown for some representative cases. 

It will be noted from Fig. 9 that for both TE- and TM-limit modes, 
the amplitude patterns for 0 < « < oo resemble those for o < —1, 
while those for —o) < o« < 0 resemble those for 1 < o. Further, the 
patterns in the first two ranges of o are quite similar to the isotropic 
patterns; those of the latter ranges depart markedly from the isotropic. 
All of these similarities are more clearly seen for the TE modes than for 
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Fig. 10(c), top, and (d), bottom — See Fig. 10. 
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the TM modes. In Part I where the propagation in guides of very large 
radius was briefly examined, it may be recalled that in both of the ranges 
0 <0 < oando < —1 the field patterns approached those of a plane 
wave rotating in the same sense as the pattern as a whole. Again, it was 
found that for —o) < « < 0 orc > 1, the field at most points in the 
guide was locally rotating in the opposite direction to that of the whole 
pattern. The points of similarity to the present case are clear and it is 
also evident that the TE mode more nearly approaches the large guide 
situation because 79 = 5.75 is much further above the cut-off radius for 
this mode than it is beyond the cut-off for the TM-mode. 
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